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Foreword 



One of the advantages of microwave techniques for diagnostics and monitor- 
ing applications is that microwave signals penetrate within dielectric struc- 
tures and they are sensitive to the presence of interior flaws and interfaces. 
Broadband microwave techniques provide additional information either 
through incorporating finite range resolution or multi-frequency material 
characterization. 

Microwave reflectometry is commonly implemented in a one-sided man- 
ner, which in turn makes it more attractive from practical point-of-view. 
The interest in broadband microwave reflectometry for materials diagnostics 
and for monitoring physical parameters of materials covers a broad realm of 
applications including: civil engineering and infrastructure, agriculture and 
medicine. Broadband microwave reflectometry is an area of engineering and 
science from which many publications have resulted over the years. 

The authors of this monograph have expertly brought together informa- 
tion from many of such papers and by many investigators as well as their 
own. Of course, this monograph does not reflect all works in this field, nor 
does it answer all questions with respect to diagnosis and monitoring appli- 
cations. However, it serves as an excellent summary of important broadband 
reflectometry approaches including the time domain reflectometry (TDR), 
the frequency domain reflectometry (FDR) and the TDR/FDR combined 
approaches. It is also important that their specific applications for the charac- 
terization of liquid materials, for monitoring of water content and for antenna 
measurements are considered in detail. 

They include simultaneous measurement of the levels and the dielectric 
characteristics of liquid materials in layered media with consideration of 
measurement accuracy improvement using appropriate probe design, custom- 
made fixtures for calibration and a targeted optimization routine. 
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Another application I would like to mention focuses on moisture measure- 
ments and includes estimation of moisture content directly from TDR and 
through TDR/FDR combined approach. Though these methods and tech- 
niques are developed for soil measurements, they can also be applied for 
varieties of materials. 

I believe this monograph will be useful to scientists, researchers and prac- 
titioners as well as students for future comprehensive studies, investigations 
and applications. 

Rolla (MO), Prof. Sergey Kharkovsky 

February 2011 Missouri University of Science and Technology 

Rolla (MO) 



Preface 



Monitoring and diagnostics are essential in many application fields: for the 
industry, for laboratory applications, as well as for countless other areas. 
Therefore, over the years, considerable research effort has been devoted to ex- 
plore innovative technologies and methods that could guarantee increasingly 
reliable and accurate monitoring solutions. In this regard, electromagnetic 
methods have attracted great interest, also thanks to their vast potential 
for nondestructive testing. In particular, broadband microwave reflectomctry 
(BMR) has established as a powerful tool for monitoring purposes; in fact, 
this technique can balance several contrasting requirements, such as the ver- 
satility of the system, low implementation cost, real-time response, possibility 
of remote control, reliability, and adequate measurement accuracy. 

On such bases, the central topic of this book is the investigation of in- 
novative BMR-based methods for monitoring applications. More specifically, 
throughout the book, the different approaches of this technique will be consid- 
ered (i.e., time domain reflectomctry - TDR, frequency domain rcflectometry 
- FDR, and the TDR/FDR combined approach) and several applications will 
be thoroughly investigated. For each considered application, particular at- 
tention will be focused on innovative strategies and procedures that can be 
adopted to enhance the overall measurement accuracy. There are many ap- 
plication areas where TDR and FDR (or the combination of the two) have 
proved useful and promising. Therefore, it comes as no surprise that the 
applications considered herein are very diverse from each other and cover 
different fields. The present book is structured as follows. 

In Chapt. [H a brief overview of the contexts in which monitoring has 
assumed a paramount importance is given. 

Chapt. [^introduces the theoretical principles that are at the basis of BMR, 
and describes the parameters of interest in this technique. Furthermore, this 
chapter introduces to dielectric spectroscopy, which is one of the pivotal ap- 
plications of BMR (in fact, as well known, many of the applications of BMR 
often descend from dielectric spectroscopy measurements). 
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In Chapt. El the TDR, FDR, and TDR/FDR combined approaches are 
thoroughly discussed, and the related advantages and shortcomings are ad- 
dressed. In particular, some effective strategies for enhancing the accuracy of 
BMR measurements; all these strategies basically aim at compensating for 
the effect of systematic errors. 

Chapt. |31 presents some innovative BMR-based solutions for the simultane- 
ous monitoring of qualitative and quantitative characteristics of liquids. This 
is quite a broad application area, and the investigated applications cover di- 
verse sectors: from the monitoring of stratified liquids (particularly useful for 
the industry of petrochemicals) to the analysis on edible liquids (useful, for 
example, for anti-adulteration control on vegetable oils). 

Chapt. [5] is focused on the qualitative analysis of granular materials and 
on moisture measurements. This last application, which is particularly useful 
in soil science and agriculture, can be extended also to different areas. In 
fact, the qualitative status in the production of agrofoods, in the production 
of inert materials, and in the fertilizer industry (just to name a few of the 
investigated application fields) is strongly related to the amount of water 
content of these materials. 

Finally, Chapt. [B] deals with the characterization of antennas through 
BMR. In particular, some useful guidelines for antenna characterization are 
presented. More specifically, the proposed procedures can help obtaining a 
higher measurement accuracy even when measurements are not performed in 
dedicated facilities (i.e., anechoic chambers). 

In all the described procedures and methods, the ultimate goal is to endow 
them with a significant performance enhancement in terms of measurement 
accuracy, low cost, versatility, and practical implementation possibility, so as 
to unlock the strong potential of BMR. 

Lecce (Italy), Andrea Cataldo 

February 2011 Egidio De Benedetto 

Giuseppe Cannazza 
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Chapter 1 
Introduction 



'Common sense is the collection of prejudices acquired by 
age eighteen '. 
Albert Einstein 



Abstract. This chapter introduces to broadband microwave reflectometry (BMR), 
which is the expression used in this book to indicate an electromagnetic technique 
used for monitoring and diagnostics purposes. The adjective broadband emphasizes 
the fact that the analysis can be carried out over a wide frequency range. A brief 
survey of the typical applications of BMR is also provided. 



1.1 The Concept Behind This Book 

Under the term monitoring, the Cambridge dictionary states the following definition 

to watch and check a situation carefully for a period of time in order to discover 
something about it. 

Similarly, the word diagnosis indicates 

a judgment about what a particular illness or problem is, made after examining it. 

Such broad, and yet unambiguous definitions seem to emphasize the urgency asso- 
ciated with the terms themselves |1|]. As a matter of fact, now more than ever, 'keep- 
ing the situation under control', for cognitive or predictive purposes, is a transversal 
need in a number of interdisciplinary sectors, especially in industrial, environmen- 
tal, civil, and agricultural areas. As a direct consequence, strict regulations have 
been enforced and guidelines have been issued that dictate specific standards (either 
qualitative or quantitative) that should be complied with, in order to ensure the suit- 
ability of the 'product' or of the 'process' . In this scenario, monitoring and diagnosis 
procedures are used for the optimization of process and for triggering possible cor- 
rective actions, or for checking the compliance of the product with pre-established 
quality standards. 

In general, a monitoring program gathers data for several purposes, such as 

• to draw comparisons against standard or target status; 

• to make comparisons between different conditions or to conduct analyses; 
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2 1 Introduction 

• to establish long-term trends in systems under observation; 

• to estimate changes with respect to a reference status. 

There are four major features that can summarize the desirable requirements that 
monitoring systems should satisfy: i) time saving; ii) energy saving; iii) economi- 
cal feasibility (or profitability); and iv) accuracy and reliability. As a matter of fact, 
the Scientific Community is in constant struggle to individuate innovative monitor- 
ing solutions that can help to simultaneously satisfy such requirements or, at least, 
achieve an optimal trade-off among them. 

Monitoring implies measuring, either directly or indirectly, a quantity that repre- 
sents the value of the parameter of interest. Most measurement techniques rely on 
applying a stimulus to the system under test (SUT) and on analyzing how the system 
responds; one of the most popular approaches consists in pulsing the SUT. A broad 
classification of measurement techniques based on pulsing includes transmission 
measurements and reflectometric measurements: the former rely on the analysis of 
the signal that is transmitted through the SUT; whereas the latter are based on the 
analysis of the signal reflected by the SUT. Reflectometric measurements typically 
require a simpler measurement setup; hence, they are generally more adaptable to a 
wider range of conditions. Depending on the specific application field, the stimulus 
signal may be acoustic, electromagnetic, or optic. 

The present book focuses on the use of electromagnetic (EM) signals as stimula: 
this approach is herein referred to as broadband microwave reflectometry (BMR). 
This expression emphasizes the fact that the analysis is performed over a wide fre- 
quency range, that theoretically goes from Hz up to the microwave region of the 
frequency spectrum (this is contrast, for example, to resonance methods which rely 
on analysis at a specific value of frequency). 

In such a context, this book is intended to provide a comprehensive overview 
on the capabilities offered by BMR: in fact, this technique arguably encompasses 
the aforementioned needs, thus allowing to achieve, simultaneously, low implemen- 
tation costs, accuracy of results, and in situ implementation in numerous practical 
applications. 

As will be detailed in the following chapters, BMR-based measurements can be 
performed either in the time domain (time domain reflectometry - TDR) or, equiva- 
lently, in the frequency domain (frequency domain reflectometry - FDR). Depending 
on the specific application, one approach may be more suitable than the other. Typi- 
cally, instrumentation operating in time domain (TD) is usually less expensive than 
instruments operating in frequency domain (FD). In fact, portable low-cost units are 
readily available on the market, thus making the TDR technique appealing for in situ 
applications. Commonly, step-like voltage function and impulse are used as stimula 
for TD measurements. 

On the other hand, reflectometric measurements performed directly in FD are 
typically carried out through vector network analyzers (VNAs). These instruments 
use sinusoidal signals as stimula, and they usually provide high measurement accu- 
racy, thanks to the possibility of directly performing calibration procedures, which 
are crucial for reducing the effect of systematic error sources. 
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An optimal trade-off (in terms of low cost, in situ applications, and measurement 
accuracy), can be achieved through a TD/FD combined approach. This particular 
approach involves standard measurements in one domain and, through a suitable 
data-processing, the corresponding information is retrieved in the other domain. 

On such bases, the aim of the book is not only to provide a clear picture of the 
most interesting state-of-the-art applications of BMR, but also to give helpful hints 
that can allow to fully exploit the potential of BMR. In particular, throughout the 
book several strategies are presented (such as calibration techniques and innovative 
processing strategies), all aimed at enhancing measurement accuracy, while still pre- 
serving low cost and possibility of practical implementation. Particular emphasis is 
given to the possibility of implementing a TD/FD combined approach that, start- 
ing from measurements involving low-cost TDR apparata, could provide accurate 
results in FD, useful for the intended monitoring purposes. The book also reports 
some significant experimental results that corroborate the proposed methods. 

It is important to point out that the explored applications all have strong poten- 
tial also for automation, which makes BMR particularly appealing for industrial 
applications. 

Finally, it is worth mentioning that, although the present book focuses on mi- 
crowave reflectometry, many of the core ideas can be readily extended to applica- 
tions of reflectometry that employ different stimulus signals. 

1.2 Survey of Typical Applications of BMR 

Thanks to the intrinsic versatility and accurate performance, applications of mi- 
crowave reflectometry are numerous and cover a wide range of fields. Just to give a 
rough idea of the broad spectrum of applications in which BMR can be successfully 
adopted as a monitoring tool, here follows an overview (which is far from being ex- 
haustive) of some of the most well-established and/or promising application fields. 

1.2.1 Electrical Components Characterization: Testing and 
Localization of Faults 

This application strictly relies on measurements of the electric impedance of the 
SUT: in fact, faults along wires or failure of electrical components cause specific 
impedance changes that can thus be taken as indicators of malfunction. Several 
reflectometry-based approaches have been investigated over the years, a clear de- 
scription of such methods can be found in |8|]. However, TDR still remains one of 
the most advantageous techniques to fulfill this task. In fact, TDR allows to directly 
measure the impedance profile of the SUT as a function of the electrical distance. 
As a result, also the spatial location of the fault becomes very straightforward. 

Locating faults along cables and wirings (in buildings, aircrafts, and so on) is the 
foremost application for which TDR was originally developed |2l,|2l|]. Recently, an 
approach to locate wire faults using reflectometry without physical contact with the 
wire conductor was presented in 13811 . In 113211 . a method based on TDR response and 
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genetic algorithm was proposed. The benefits of the TDR/FDR combined approach 
have been considered for diagnostics on electric power cables (also to monitor in- 



cipient defects 13711 *): in B33H . faults are localized by monitoring the cross-correlation 
characteristics of the observed signal in both TD and FD, simultaneously. BMR is 
also used for the characterization of power electronics systems [39]. Based on the 
same principle, TDR has also been used for testing of instrumentation loops in gas 
turbines during scheduled maintenance H20h . 

1.2.2 Measurements on Soil for Agricultural and Geotechnical 
Purposes 

Soil measurement is another pivotal application area of BMR. In particular, soil 
water content measurements represent one of the major practical applications of 
this technique, with uses that focus on environmental monitoring 13411 and on the 
optimization of irrigation cycles and water resource management [26]. Moisture 
measurements through BMR start from the estimation of the relative dielectric per- 
mittivity of the sample under test: in fact, even a small amount of water (which 
has a relative dielectric permittivity of approximately 80) increases significantly the 
overall dielectric permittivity of soils (which, in dry conditions, typically have much 
lower relative dielectric permittivity). 

TDR is also used in geotechnical engineering and environmental control for mon- 
itoring liquefaction of soils 128(1 and for the detection of organic pollutants in sandy 
soils Irk . 

BMR is particularly useful also for monitoring the static electrical conductivity of 
moistened soils. In fact, BMR measurements are performed in the whole frequency 
range between Hz and the maximum frequency of the used instrument; as a direct 
result, the problems of direct current (DC) measurements of traditional methods is 
overcome. 

Finally, resorting to TDR, by monitoring the deformation of coaxial cables it is 
possible to monitor landslides Ifi Il2n . and to assess distributed pressure profiles 



(useful for evaluating the strength of the ground for constructions B27I1 ). 



1.2.3 Civil Engineering and Infrastructural Monitoring 

Reflectometric techniques are providing promising results also for monitoring ap- 
plications in civil engineering. For example, an FDR-based resonant method for the 
accurate evaluation of the depth of long and shallow surface damages was proposed 
in II 1411 , Other nondestructive reflectometry-based monitoring solutions include the 
detection of chloride ingress and the evaluation of its distribution in reinforced con- 
crete structures (useful for monitoring corrosion) 12311 : the determination of the 
material content of concrete (useful to infer the compressive strength of concrete) 
yQ; and the determination of water/cement content in fresh Portland cement-based 



materials 11811 . 
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Additionally, in [36] a TDR-based system for sensing crack/strain in reinforced 
concrete structures was presented; fault location on concrete anchors was investi- 
gated in 19D; a TDR-based solution for monitoring cement hydration was described 
in [11], and so on. 



1.2.4 Dielectric Spectroscopy Measurements 

BMR has also been extensively used for dielectric spectroscopy of materials, which 
involves measuring the frequency-dependent dielectric permittivity of the consid- 
ered material. Typically, this kind of measurements relies on the combination of 
reflectometric measurements and mathematical models that describe the expected 
dielectric behavior of the considered material. Dielectric properties are of paramount 
importance because they can be related to other non-electric properties of the SUT, 
thus obtaining useful information. 

It would be impossible to list all the specific materials for whose characteriza- 
tion BMR has proven useful; however, some of the materials investigated through 



this technique include alcohol mixture s |2|| , mixtures of ester with alcohol I30LI31I1 . 



poly(propylene glycol)water mixtures B29|]. chlorobenzene with n-methylformamide 



1 2211 . Nylon-11 with phenol derivatives 113H , powders (e.g., zeolite powders [35]). 
Dielectric characterization through BMR has also been extensively used for food 
analysis Biol and for application on biological materials (e.g., for breast tissue anal- 
ysis B25I1 . for diagnosing skin cancer II 1511 . for cerebral edema-related studies 11011 . 
etc.). 



1.2.5 Applications in Industrial Monitoring 

Typical applications of BMR for industrial use include measurement of liquid lev- 
els in tanks or containers. Basically, the level of the liquid is inferred from the 
travel time of the EM signal along the probe immersed in the monitored liquid [19]. 
Thanks to the capability of discriminating interfaces between materials with dif- 
ferent permittivity, TDR can be used for the localization of multiple interfaces in 
layered media lfll24ll. 

Other applications of BMR include monitoring of levels in dykes and rivers; 
moisture evaluation of agrofoods (such as cereals and coffee |5[]); monitoring of 
the osmotic dehydration process for fruit and vegetables; and other quality-control 
applications that resort to measurement of dielectric characteristics. 



1.3 Organization and Content of the Book 

The book comprises two main parts: one is more theory-oriented, whereas the other 
is more practice-oriented. The first part (from chapter two through chapter three) re- 
calls the theoretical background of BMR measurements and presents the strategies 



6 1 Introduction 

that are used for enhancing measurement accuracy. The second part (from chap- 
ter four through chapter six) is a representative collection of experimental results 
that demonstrate the potential of BMR for various monitoring tasks. Let us see the 
content of each chapter more in detail. 

In the second chapter, the basic theoretical principles behind BMR are recalled. 
First, a preliminary introduction on the transmission line theory is given. Second, 
the most important electrical parameters related to BMR are introduced. Particular 
attention is given to the quantities that are directly measured through BMR, namely 
the reflection coefficient in time domain (p(t)) and the reflection scattering parame- 
ter (S n(f)). Finally, the basic concepts behind dielectric spectroscopy are recalled. 

The third chapter describes the 'actors' involved in typical BMR measurements. 
First, TDR and FDR are presented and the related instrumentation is fully described. 
Successively, the TDR/FDR combined approach is described in detail: as aforemen- 
tioned, this approach can help exploit the benefits of both TDR and FDR, without 
necessarily employing two different measurement setups. It is worth mentioning 
that great emphasis is given on this issue, since it represents the solution that can 
ultimately lead to the optimization of the cost-benefit ratio in the implementation 
of the systems. Particular attention is also dedicated to the design of the probes, 
which is crucial for obtaining accurate results. In the third chapter, also the methods 
for measuring the dielectric characteristics of materials starting directly from the 
analysis of TDR waveforms are presented. Additionally, the major error sources in 
BMR-based measurement systems are described in detail: specific focus is given to 
the strategies for enhancing measurement accuracy. 

The fourth chapter deals with the quantitative and qualitative characterization of 
liquid materials. In this chapter, BMR-based methods for measuring simultaneously 
the levels and the dielectric characteristics of liquid materials are presented. Results 
show that all the proposed methods have strong potential for practical implementa- 
tion in the field of industrial monitoring. 

First, an approach based solely on the analysis of TDR measurements is described 
and validated. Secondly, a further enhancement of this method is accomplished by 
resorting to the combination of TDR with FDR. In this second measurement method, 
the accuracy is enhanced also through the adoption of a transmission line modeling 
of the measurement cell and through the realization of a custom-made fixture that 
allows performing an ad hoc short-open-load (SOL) calibration. 

Finally, BMR-based measurements of the dielectric characteristic are extended 
to edible liquids: in particular, vegetable oils are considered. This was done in view 
of the possibility of performing quality and anti-adulteration control. 

The fifth chapter focuses on BMR applications for monitoring water content and 
static electrical conductivity of granular materials, focusing on applications on soils 
and on agrofood materials. 

First, a TDR-based method for inferring water content from measurements of the 
apparent dielectric permittivity is presented. This approach, which relies on the in- 
dividuation of so-called calibration curves, is discussed in detail and a metrological 
assessment is provided. 
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Successively, a method that takes into account the frequency-dependence of the 
dielectric permittivity of the moistened granular material (considering the permit- 
tivity of each single constituent) is presented. This application is also used as a 
test-case for validating an innovative calibration procedure that becomes especially 
useful when the traditional SOL calibration cannot be performed. Furthermore, the 
adoption of antennas in place of the traditional probes is discussed, thus assessing 
the possibility of achieving a noninvasive approach. 

At the end of chapter fifth, two innovative strategies for enhancing and simpli- 
fying TDR-based measurements of electrical conductivity (typically used in soil 
science) are also addressed. 

Finally, the sixth chapter provides some significant examples on the use of 
TD/FD combined approach for the characterization of antennas. In particular, some 
guidelines are given for obtaining accurate results even without resorting to expen- 
sive facilities such as anechoic chambers. The presented procedure can be extended 
to other electronic devices or systems as well. 
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Chapter 2 

Basic Physical Principles 



'If we knew what it was we were doing, it would not be 
called research, would it?'. 
Albert Einstein 



Abstract. In this chapter, the basic concepts behind broadband microwave reflec- 
tometry (BMR) are recalled. First, a brief overview of the transmission line theory 
is provided, and the most common electromagnetic structures are introduced. Sec- 
ondly, the major parameters that are used to characterize electrical networks are 
introduced, and the related theoretical background is briefly discussed. Finally, a 
short overview on dielectric spectroscopy is provided, thus anticipating its connec- 
tion with refiectometric measurements. 



2.1 Transmission Line Basics 

In electric circuits, when the wavelength of the propagating signal is large compared 
to the physical dimensions of the system, the electrical characteristics of the system, 
at a given time, can be assumed to be the same at all points {lumped-element model). 
On the other hand, when the the physical dimensions of the system are comparable 
to the wavelength of the propagating signal, the dimensions of the cables, connectors 
and other components cannot be ignored: in such cases, it is particularly useful to 
model the system through transmission line (TL) segments. 

TLs are typically used to transfer a signal from the generator to the load, by 
guiding the electromagnetic (EM) signal between two conductors |5[1. A TL is a 
distributed-parameter network and must be described by circuit parameters that are 
distributed throughout its length. For the purposes of analysis, a TL can be mod- 
eled as a two-port network. The model represents the transmission line as an infinite 
series of two-port elementary components, each representing an infinitesimal seg- 
ment of the transmission line. The elementary section of a TL can be modeled as 
shown in Fig. 12.11 This model includes four parameters, referred to as primary line 
constants, which are generally defined 'per unit length' . The primary line constants 
are the series resistance (R), the series inductance (L), the shunt conductance (G), 
and the shunt capacitance (C). For a uniform transmission line, the primary line 
constants do not change with distance along the line. These constants are used to 
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Fig. 2.1 Equivalent circuit 
modi 

line 



model for a transmission m ^L J^^ 



define the secondary line constants, namely the propagation coefficient (y) and the 
characteristic impedance (Zo). 

For an infinitely long line, the characteristic impedance Zq (which is defined as 
the ratio of voltage V to current / in any position) can be written as 






where (0 = 2nf is the angular frequency and i 2 = — 1 . 
The propagation coefficient is given by 



y=y/(R + i(»L)(G + i(oC). (2.2) 

It is useful to separate the imaginary part (/3), which gives the phase-shift coefficient, 
from the real part (a), which gives the attenuation coefficient: 

_R_ GZo 

2Z 2 



P^coVLC. (2.4) 

For a lossless TL (i.e., when R = and G = 0), Zq can be written simply as 

Z = ^L/C. (2.5) 

From ( 12.5b . it can be seen that for a lossless TL, the characteristic impedance is 
purely resistive, although given by reactive elements (C and L). It is important to 
point out that this does not mean that the line is a resistance. 

In the following subsections, the most common types of TLs are considered, 
namely coaxial, two-wire, and microstrip. 

2.1.1 Coaxial Transmission Line 

Coaxial lines are made of a central conductor with diameter a and a hollow outer 
conductor with inner diameter b. The space between the conductors is usually filled 
with a dielectric material: the electric and magnetic fields are confined within the 
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Fig. 2.2 Schematic of the 
cross section of a coaxial 
line 



dielectric. Fig. 12. 21 shows the schematization of the cross section of a typical coaxial 
line. The capacitance and the inductance per unit length are 

^ _2ne 

(-■coax — u V-^-OJ 

icoax = z— In- (2.7) 

In a 

where jjl is the permeability of the transmission line medium; and e is the dielectric 
permittivity of the dielectric material. Therefore, the impedance per unit length is 



^E * fe„*. ( 2.8) 

Qoax 2n\j £ a 

For many materials, the permeability is equal to that of free space (i.e., jx = }Xq = 

4n x 10 _7 Hm _1 ). Additionally, considering that e = £n£ r , where e% = 8.854 x 

10 _12 Fm _1 is the dielectric permittivity of free space, for practical purposes, l !2.8b 

becomes 

60 b 

Z coax = —=ln- . (2.9) 

y/Sr a 

Typically, coaxial line can support transverse electromagnetic (TEM), transverse 
electric (TE), and transverse magnetic (TM) modes. In radio-frequency applications 
up to a few GHz, the wave propagates in the TEM mode only. However, at frequen- 
cies for which the wavelength (in the dielectric) is significantly shorter than the cir- 
cumference of the cable, TE and TM modes can also propagate: when more than one 
mode can exist, bends and other irregularities in the cable geometry can cause power 
to be transferred from one mode to another. For many microwave applications, most 
of coaxial lines are designed to work at the TEM mode [4]. Higher-order modes 
are prevented from propagating, when the frequency is below the cutoff frequency, 

/c— o- 

/c-o= ,„J; t= (2.io) 
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where c = 3 x 10 8 ms _1 is the velocity of light in free space; /i r is the relative 
permeability of the transmission line medium; and e r is the relative permittivity 
of the dielectric material. 



2.1.2 Bifilar Transmission Line 

The configuration of this TL comprises a pair of parallel conducting wires separated 



by a constant distance (Fig. 12.31 ). For a bifilar transmission line [25] 



ne 
Cbif= ln(2D/d) (2J1) 

L w =-ln— -. (2.12) 

% a 

Hence, the characteristic impedance of a bifilar TL is given by 

120 2D 
Z bif =— Z«— (2.13) 

y/F T d 

where e r is the relative permittivity of the dielectric between the wires/rods; D is 
distance between the rods; and d is the diameter of each rod. 



Fig. 2.3 Schematic of a 
two-wire line 



2.1.3 Microstrip Line 

The geometry of microstrip lines is very simple: a thin conductor and a ground 
plane separated by a low-loss dielectric material. The fabrication typically relies 
on printed circuit techniques. The use of a substrate with a high dielectric constant 
reduces the fringing field in the air region above the conductor. Because the electric 
field lines remain partially in the air and partially in the substrate, microstrip lines 
do not support pure TEM mode, but a quasi-TEM mode |3fl. For 0.05 < w/h < 1, 
the characteristic impedance is given by 
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Fig. 2.4 Schematic of a 
microstrip line 



conductor strip 




ground plane 



cross section t « w 
t«h 
a »w 
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whereas for 1 < w/h < 20 
120^: 



'£eff 



— + 1 .393 + 0.667/m — ) + 1 .49 

h \h 



-l 



(2.14) 



(2.15) 



where e e ff is the effective permittivity; w is the width of the microstrip; and h is the 
thickness of the substrate [126]. 



2.2 Reflected Waves 

Let us consider a Zo-impedance TL, terminated with a generic load impedance, Zl 
(Fig. 12.5b . When Zl does not match Zq, part of the incident wave is absorbed by the 
load, and the rest is reflected along the line. 



Fig. 2.5 TL connected be- 
tween the generator and the 
load 



I- Zii 



*.(\> 



generator ■ transmission line ' load 



At any point on the line, the reflected voltage is given by B25I1 

v R = rvi 



(2.16) 
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where F is the reflection coefficient and VJ is the incident voltage. F is given by 

^ Z L -Z 



Zl + Zo 
with \r\ < 1. By rearranging ( I2.17l i 



(2.17) 



i+r 

Zl=Z 0t -^ (2.18) 

Generally, F is a complex quantity and it is often defined through magnitude (p) 
and angle phase (<p): 

r = pl(j). (2.19) 

It is worth noting that usually both F and p are referred to as reflection coefficient; 
however, the former is a complex quantity, whereas the latter is scalar. As long 
as there is no change in the impedance of the TL, the magnitudes of the incident 
and reflected voltages do not change with position; therefore, p does not change 
with position. On the other hand, the phase of the reflection coefficient changes as 
position changes 12911 . 

When a sinusoidal signal propagates down a transmission line and the termina- 
tion of the line is not matched, the incident voltage and the reflected voltage create 
a typical interference pattern. The envelope of the sinusoidal voltage will maintain 
a constant shape (standing wave). The magnitude varies with distance, but the volt- 
age at each point on the line varies sinusoidally. The voltage standing wave ratio 
(VSWR) is the ratio of the maximum V max and minimum V m j n of the envelope: 

VSWR = ^- . (2.20) 

* min 

The VSWR is related to the magnitude of the complex reflection coefficient by the 
following equation: 

rara = I±H. (221) 

VSWR is greater than or equal to one: VSWR = 1 when no mismatch occurs; 
VSWR — ► °° for an open or short circuit. 

2.3 Characteristic Parameters of Electrical Networks 

An electrical network is typically used to indicate a circuit (either single or multi- 
component) or a generic TL. Generally, an electrical network is represented as a 
black box with input and output terminals, and the analysis of the network is done 
at the terminals. A network can be characterized through different parameters that 
relate current to voltage; the most common parameters are summarized in Table l2.3l 
When dealing with multi-port networks, it is particularly useful to represent network 
parameters through matrices. The choice of the most suitable representation depends 
not only on which the characteristics of interest are, but also on the frequency range 
of analysis and, obviously, on the available measurement instrumentation. 



2.3 Characteristic Parameters of Electrical Networks 



17 



Table 2.1 Typical parameters used for describing electrical networks 



symbol 


matrix 


ABCD* 


transmission 


Z 


impedance 


Y 


admittance 


H* 


hybrid 


G* 


inverse hybrid 


S 


scattering 



* Cannot be used for networks with more than 2 ports. 



At high frequencies, the influence of parasitics (i.e., capacitances, cable induc- 
tance, undesired coupling effects) becomes more significant; hence, it becomes in- 
creasingly difficult to measure voltages and currents. In this case, it is preferred to 
measure power, and if necessary, to successively extrapolate the other parameters of 
interest. In this regard, scattering parameters (S- parameters) are particularly conve- 
nient; in fact, S-parameters relate to signal flow rather than to voltages and currents 
directly. In particular, the 'measured quantities' are traveling waves. 



Fig. 2.6 Schematization of 
an Alport network 
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Considering an Alport network (Fig. 12. 61 ). scattering variables at the generic port 
n are defined in terms of the port voltage V„, port current /„, and a normalization 
impedance Zq. The voltage and the current can be synthesized by an incident scat- 
tering variable, a„, and a reflected scattering variable b n , given by 



a n 



Vn 



2V2o 



InVZa 



(2.22) 



Vn 



2VZn 



I/iVZq 



(2.23) 



The scattering parameters are measured in terms of a„ and b„. For a Alport network, 
there are A^ 2 scattering parameters, Sn, with i,j= l,...,N. 

The S-parameters where i — j are referred to as reflection scattering parameters. 
They are practically measured as the ratio between the outgoing and incoming wave 
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at the i' h port, when a matched generator is connected to the i rl 
ports are terminated with a matched load (i.e., aj = 0,V/ 7^ j) 



ort and all the other 



por t 



Sit — — L=0,V#i- 



(2.24) 



The terms S, ; - are referred to as transmission scattering parameters. They are mea- 
sured as the ratio of the outgoing wave at the i th port (b,) and the incoming wave at 
the /* port (aj) when a matched generator is connected at port j and the other ports 
are connected to a matched load: 



Oij — — | fl ,=o,v;/;- 



Or 



(2.25) 



Often, the matrix notation is used 



b N 






Sm 



Sin 


X 


fl2 


Snn_ 




. aN . 



(2.26) 



For a reciprocal network, the scattering matrix is symmetrical, i.e. S mn = S nm 
(m,n = l,2,...,N). As can be seen from the previous equations, to measure the 
generic scattering parameter Stj, which involves ports i and j, the remaining ports 
have to be terminated with matched load^l (rather than with short or open circuits, as 
must be done with other network parameters). This is great advantage of character- 
izing electrical networks through S-parameters; in fact, at microwave frequencies, 
open and short circuits are more difficult to implement and may lead to instability 
(severe distortion or oscillation may occur) 112911 . 



Fig. 2.7 Two-port network 
scheme 
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Fig. 12.71 shows the schematic of a two-port network: Iy and I2 are the currents 
at the terminals; V\ and V2 are the voltages at the terminals; and a and b are the 
traveling waves. The square of the coefficient represents the power: \a\ \ and \b\ \ 
can be considered as the power entering port 1 and exiting port 7, respectively. For 
two-port networks, the system can be modeled by two equations: 



Clearly, when performing reflection measurements (i.e., 5„), all ports except the i port 
should be terminated with a matched impedance. 
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b\ =S\\a\ +5i2«2, (2.27) 

h = S2iai +S22C12 ■ (2.28) 

The scattering parameters are therefore measured as 

Su = — k=o, (2.29) 

a\ 

S 2 i = — k=o, (2.30) 

Si2=— k=o, (2-31) 

fl2 

5 22 = -k=o- (2.32) 

<22 

The potential of S-parameters for monitoring and diagnostics purposes will be clearer 
in the following chapters. In fact, it will be shown that, when the system to be moni- 
tored can be considered as an electrical network, the S-parameters of the system can 
disclose useful information on the system itself. In particular, the scattering parame- 
ters can be used to extrapolate the dielectric characteristics of the SUT The retrieved 
information, although related to the dielectric behavior, may be appropriately asso- 
ciated with several quantitative and/or qualitative properties of the sample. 



2.4 Reflectometry Measurements for Dielectric Spectroscopy 

Dielectric spectroscopy, also called impedance spectroscopy, measures the dielectric 
properties of a medium as a function of frequency. It is based on the interaction of 
an external field with the electric dipole moment of the sample, often expressed by 
permittivity. This technique measures the impedance of a system over a range of fre- 
quencies, and, therefore, the frequency response of the system, including the energy 
storage and dissipation properties. Almost any physico-chemical system possesses 
energy storage and dissipation properties that can be analyzed through dielectric 
spectroscopy. 

There are several approaches that can be used to measure the dielectric proper- 
ties of materials, such as open-ended probe technique, free-space transmission tech- 
nique, cavity resonant systems, time domain reflectometry, and so on. An overview 
of some of the most common methods for measuring dielectric properties of mate- 
rials can be found in B27I1 . 

In particular, dielectric spectroscopy performed through microwave reflectom- 
etry, has gained enormous popularity over the past few years and is now being 
widely employed in a wide variety of scientific fields that go from food analysis 
12011 to biomedical applications J24lll9L ll2ll . As will be detailed in the following 
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chapters, for microwave reflectometry-based dielectric spectroscopy of materials, 
either the direct frequency domain (FD) approach or the time domain/frequency 
domain (TD/FD) combined approach can be adopted 12811 . 



2.4.1 Dielectric Relaxation Models 



Materials can be classified according to their dielectric permittivity. The dielectric 
behavior of a material is suitably described by the frequency-dependent complex 
relative permittivity, £,.* (/), which can be written as follows |4J]: 



<?(/) = «ft/)-ik(/) 



Ob 



2nf£o 



(2.33) 



where i 2 = — 1, £o = 8.854 x 10 _12 Fm _1 is the dielectric permittivity of free space, 
/ is the frequency, £ r '(/) describes energy storage, ef (/) accounts for the dielec- 
tric losses, and On is the static electrical conductivity (which is related to the ionic 
losses). Generally, when (ao/2nfe' r ) »1, the material can be considered a good 
conductor. On the other hand, when (ao/2nfe' r ) <<1, the material is considered a 
dielectric (i.e., lossless or low-loss material). Materials with a large amount of loss 
inhibit the propagation of electromagnetic waves. Those that do not fall under either 
limit are considered to be general media. A perfect dielectric is a material that has 
no conductivity, thus exhibiting only a displacement current. 

The dielectric relaxation refers to the relaxation response of a dielectric medium 
to an external electric field of microwave frequencies. This relaxation is often de- 
scribed, in terms of dielectric permittivity, through the so-called relaxation models. 
Fig- EU shows a schematization of typical relaxation phenomena. 

There are several mathematical models (obtained either theoretically or empiri- 
cally) that are used to describe the dielectric properties of materials. 

For example, for pure polar materials, the equation known as Debye model de- 
scribes the dielectric permittivity 190 : 



Fig. 2.8 Generic dielectric 
permittivity spectrum over a 
wide range of frequencies. 
The figure shows the main 
processes that occur: ionic 
and dipolar relaxation, and 
atomic and electronic res- 
onances at higher energies 
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<(/) 



■u£ 



(2.34) 



where £<*, is the dielectric constant at frequency so high that molecular orientation 
does not have time to contribute to polarization (i.e., at infinite frequency) J2D; £ s 
is the static dielectric permittivity; f Te \ is the relaxation frequency (defined as the 
frequency at which the permittivity is (e s + £»)/2) 11611 : and / is the frequency. 

As a matter of fact, only few materials of interest behave as pure polar materials 
with a single relaxation frequency; therefore, other models have been developed to 
describe more accurately the dielectric behavior of materials \T\. The Cole-Cole 
model, with symmetric distribution of relaxation times, is described through the 
following equation: 

e r * (/) = &• + £s ~ £ °;_^ (235) 
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where/3 is a parameter that describes the spread in relaxation frequency (0 < /3 < 1) 
Ild.ll7n . When /3 = 0, the Cole-Cole model reduces to the Debye model. 

A further generalization of the Debye model is given by the Havriliak-Negami 
formula SQ: 

e* (/) = 6. + T £s ~ £ : nxia (2.36) 



1 



('£) 



(1-/3)' 



where a describes the broadness of the permittivity spectra. 

Finally, the Cole-Davidson model with asymmetric distribution of relaxation 
times follows for /3 = and for < a < 1 |fl]: 
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(2.37) 



When the static electrical conductivity is not negligible, an additional term must 
be included in the equations above; therefore, for example, the Cole-Cole model 
described by (12.35b becomes 
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(2.38) 



The five parameters are referred to as Cole-Cole parameters: they are different for 
each material and they define its spectral signature. The Cole-Cole relaxation model 
has proved suitable in describing relaxation phenomena of various materials, espe- 
cially for lossy media talSD. an d f° r a wide number of materials (e.g., concrete, 
paper, petrochemicals, and even blood) H 1 3L 1 1 5f1 . 
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In particular, the estimation of the Cole-Cole parameters of foods (e.g., egg white, 
butter, bacon, etc.) has received considerable attention, thanks to the possibility of 
associating the dielectric parameters to the quality status of the products/materials 

As a matter of fact, the evaluation of such parameters is not always a straight- 
forward task, especially when dealing with complex materials, such as mixtures, 
powders, etc.. 

For the sake of example, Table 12.21 shows the Cole-Cole parameters for some 
common liquids, as reported in 01111 . The Cole-Cole parameters of these and of other 
liquids are typically used as reference for calibrating and validating experimental 
setups for dielectric measurements. Fig. |2.9l shows the complex dielectric spectrum 
of water, as reported in 11811 . 



Table 2.2 Values of the Cole-Cole parameters for some common liquids, as reported in II ill 



material 
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(ps) 
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air 


1.0005 








chloroform 


4.82 


2.28 


6.25 




ethyl-acetate 


6.04 


2.48 


4.4 


0.06 


toluene 


2.40 


2.25 


7.34 
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Fig. 2.9 Complex dielec- 
tric spectrum of water at 
20°C. The dashed line indi- 
cates the extrapolated high 
frequency permittivity, £«,, 
of the principal relaxation 
of water. The dotted line 
shows the extrapolated high 
frequency permittivity of 
an additional relaxation, 
with relaxation frequency 
1/2MT 1 , around 1000 GHz 
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Chapter 3 

Broadband Reflectometry: Theoretical 

Background 



'A scientist 's aim in a discussion with his colleagues is not 
to persuade, but to clarify'. 
Leo Szilard 



Abstract. In this chapter, the basic approaches of broadband microwave reflectome- 
try are described in detail. More specifically, time domain reflectometry (TDR) and 
frequency domain reflectometry (FDR) are presented and the involved instrumenta- 
tion is fully described. Successively, the FDR/TDR combined approach is described 
in detail: this approach can help exploit the benefits of both TDR and FDR, without 
necessarily employing two different measurement setups. Additionally, since the 
sensing element (or probe) plays a major role in all the aforementioned approaches, 
a comprehensive description of its design and of the corresponding performance is 
given. Finally, the basic principles leading to the possibilities of enhancing accuracy 
in BMR measurements are presented. 



3.1 Broadband Microwave Reflectometry: Theoretical 
Background 

Broadband microwave reflectometry (BMR) is a powerful technique that can be ef- 
fectively employed for a number of practical applications; in particular, the high ver- 
satility, the real-time response and the potential for practical implementation have 
contributed to the success of microwave reflectometry for monitoring purposes. 

Generally, in BMR, a low-power electromagnetic signal is propagated into the 
system under test (SUT): the analysis of the reflected signal along with specific 
data-processing are used to retrieve the desired information on the SUT J32I1 . 

Two main elements are involved in BMR measurements: 

1 . the instrument for generating/receiving the electromagnetic (EM) signal, and 

2. the measurement cell, which includes the sensing element (or probe) and the 
SUT. 

Microwave reflectometry-based measurements can be performed either in time do- 
main (time domain reflectometry - TDR) or in frequency domain (frequency domain 
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reflectometry - FDR). Depending on the specific application, one approach may be 
more suitable than the other. 

TDR instrumentation is usually less expensive than instruments operating in fre- 
quency domain. Nevertheless, instruments operating directly in frequency domain, 
such as vector network analyzers (VNAs), despite being more expensive, usually 
provide higher measurement accuracy. An optimal trade-off (in terms of low cost 
and measurement accuracy), can be achieved through a time domain/frequency do- 
main (TD/FD) combined approach: in this case, measurements are carried out, for 
example, in TD and the corresponding FD-information is retrieved through ap- 
propriate processing of the acquired data, thus discovering additional information 

lira. 



It is worth pointing out that a substantial advantage in operating in FD relates to 
the possibility of performing calibration procedures (as will be detailed in Sec. 13.41 ). 
thus enhancing the final measurement accuracy. This procedure can be implemented 
in the TD/FD approach; consequently, the combination of the two approaches is 
highly regarded as a powerful tool for guaranteeing at the same time low cost of 
the setup and adequate measurement accuracy. This aspect, along with the detailed 
description of the specific implemented procedure, will be addressed in Sect. 13.41 In 
the following subsections, the basic principles behind all the three approaches are 
recalled, and pros and cons are discussed. 



3.2 Time Domain Reflectometry (TDR) 



In TDR measurements, the EM stimulus is usually a step-like signal that propagates 
along the probe, through the SUT. The basic block diagram of a typical TDR system 
is shown in Fig. 13.11 It consists of a step-pulse generator, a sampling scope, and a 
signal analyzer. The step generator produces a positive-going incident wave that is 
applied to the SUT. The step travels down the transmission line at the velocity of 
propagation of the line. Any impedance variation causes the partial reflection of the 
propagating signal. If the load impedance is equal to the characteristic impedance 
of the line, no wave is reflected, and the oscilloscope displays only the incident 
voltage step. If a mismatch exists at the load, part of the incident wave is reflected 
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Fig. 3.1 Functional block diagram for a time domain reflectometer | 
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14|]. The reflected signal is acquired by the oscilloscope, and its voltage amplitude is 
displayed as a function of time (or as a function of the traveled electric distance). 

The ratio between the amplitude of the reflected signal, v re a(f), and the amplitude 
of the generated signal, v; nc (/), gives the value of the reflection coefficient in time 
domain, p(f) (jjl: 

P(0 = ^f (3-D 

Vinc(0 

where -1 < p(t) < +1. 

It is important to point out that, most often, the oscilloscope functionality is in- 
tegrated with the signal generator within one single instrument (which, as will be 
detailed later in this chapter, can be either portable or benchtop). 

Clearly, the behavior of p(t) is strictly associated with the impedance variations 
along the electrical path traveled by the electromagnetic (EM) signal. As an exam- 
ple, Fig. 13.21 shows the schematization of the TDR waveforms observed when the 
SUT exhibits purely-resistive behavior. 

If an increase of impedance is encountered (i.e., Zl > Zq), then a positive step 
is observed and the reflection coefficient will be positive. Conversely, if a decrease 
of impedance is encountered (i.e., Zl < Zq), then a negative step is observed and 
the reflection coefficient will be negative. In particular, if the impedance of the load 
equals the characteristic impedance of the line (i.e., Zl = Zq), no wave will be re- 
flected (p = 0); hence, the TDR display on the scope is a flat line. If an open circuit 
is encountered, then the reflected voltage will equal the generated voltage and the 
reflection coefficient will be +1. If a short circuit is encountered, the reflection co- 
efficient will be -1. It goes without saying that the actual value of an 'unknown' Zl 
may be inferred from the reflection coefficient displayed by the TDR unit. 

Also of interest are the reflections produced by complex load impedances. Four 
basic examples of these reflections are shown in Fig. 13.31 where the incident step 
voltage is indicated as £;. A direct analysis (in TD) of the reported waveforms in- 
volves evaluating the reflected voltage at r = and at t = °°, and assuming any 
transition between these two values to be exponential. 



v(r) 



z, > z„ 



Vin/4 



Z,= (p = -l) 



Fig. 3.2 Schematization of TDR waveforms for purely-resistive terminations 
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Where k = —!— 
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Fig. 3.3 TDR waveforms when reactive components are present |4]. The incident step signal 
is indicated as E\ 
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For the sake of example, let us consider the series R-L case |4|]. Also, for simplic- 
ity, let us assume that the time is zero when the reflected wave arrives back at the 
monitoring point. 

At t = 0, the reflected voltage is +E\. This is because the inductor will not accept 
a sudden change in current: it initially behaves like an infinite impedance; hence, at 
t = 0, p will be unit. 

Successively, the current in L builds up exponentially and its impedance drops 
toward zero. Therefore, at t = °°, E r (t) is determined only by the value of R: 

The exponential transition of E T (t) has a time constant (t) determined by the effec- 
tive resistance seen by the inductor. Since the output impedance of the transmission 
line is Zq, the inductor sees Zq in series with R: 

(3.3) 



R + Za 



A similar analysis can be carried out for the other three cases reported in Fig. 13.31 

It must be pointed out that an in-depth analysis of the impedance characteristics 
of the SUT should rely also on a FD approach (in particular, through the evaluation 
of the reflection scattering parameter). However, as seen from Fig. 13.31 the TDR 
waveforms can provide a concrete idea on the impedance-behavior of the SUT. 

This approach is particularly useful for TDR-based individuation of faults along 
cables. In fact, faults do show specific impedance behavior, whose meaning can be 
inferred directly from the TDR waveforms. 

3.2.1 Typical TDR Measurements 

When the TDR signal propagates along a transmission line, impedance mismatches 
cause the reflection of a portion of the signal. The TDR signal's propagation veloc- 
ity, v, is related to the relative dielectric permittivity (e r ) of the medium, which is 
assumed to be lossless or at least with negligible conductivity, and to the relative 
magnetic permeability (|i r ), which is equal to 1 for most materials jj]: 

" = -£= (3.4) 



where c is the speed of light in free space (c = 3 x 10 8 m s _1 ). 

For the sake of example, let us consider a typical three-rod probe, whose elec- 
trodes have length L. A typical TDR waveform for such a probe, immersed in a 
generic homogeneous dielectric, is represented in Fig. 13.41 The analysis of TDR 
waveform directly leads to the evaluation of L app (also called electric distance), 
which can be directly associated to the dielectric characteristic of the medium. In 
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fact, L ap p can be considered as the distance that would be traveled by the EM signal, 
in the same interval of time, if the signal were propagating at a speed c. 

Typically, in TDR analysis, the apparent distance is evaluated through the indi- 
viduation of the beginning and end of the probe, thus it can be associated to the 
physical length, L, through the following equation: 



L. 



app 



(3.5) 



where t t is the travel time (round-trip time taken by the signal to travel between the 
beginning and end points), and e app is referred to as apparent dielectric permittivity 
of the material in which the probe is inserted 14J]. 

A more rigorous expression for £ app is given by the following equation 13711 : 
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where £o = 8.854 x 10 _12 Fm _1 is the dielectric permittivity of free space, / is the 
frequency, e' Y {f) describes energy storage, e"(f) accounts for the dielectric losses, 
and (Jo is the static electrical conductivity. 

However, when dealing with low-loss materials, the imaginary part of the com- 
plex permittivity can be neglected. Additionally, considering low dispersive materi- 
als, the dependence of £^(f) on frequency can be considered negligible. Under these 
circumstances, £ app (/) can be considered approximately constant: 



£a PP (/) =ei(f) = const. 



(3.7) 



On the basis of the above described theory, TDR-based dielectric measurements 
simply rely on the individuation of the points corresponding to the beginning and 
to the end of the probe. To individuate these two points in the TDR waveform, 
different approaches based on the so-called tangent method can be used I24L I41I1 . 
However, as can be intuitively deduced from Fig. 13.41 L app can be easily evaluated 
through the derivative of the TDR waveform, which typically exhibits prominent 
peaks in correspondence of the probe-beginning and probe-end sections. Indeed, 
not only is this 'derivative method' simple and quick, but it also proves particularly 
advantageous when significant impedance changes are masked by other impedance 
variations (for example, by the transition coaxial cable/rods in Fig. 13.4b . as will be 
clarified later in this book. 



3.2.2 Typical TDR Instrumentation 

As aforementioned, the pivotal components of a TDR instrument are a high- 
bandwidth oscilloscope and a high-speed pulse generator: the characteristics of 
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Fig. 3.4 Typical TDR waveform (black curve) and corresponding first derivative (red curve) 
for a 30 cm-long three-rod probe in air, obtained through a Tektronix TDR80E04. The am- 
plitude of the incident step-like signal is 0.25 V. The first derivative of the TDR waveform 
emphasizes the impedance variations. The section k = 1 indicates where the sensing portion 
of the considered probe begins 



these two components have the major influence on the performance of the instru- 
ment. In particular, the frequency content of the signal, /b, is approximately given 
by the following equation: 

0.35 
h = — (3-8) 

h 

where f r is the rise-time of the generated step signal. 

A higher bandwidth of the signal makes the TDR measurements significant 
in a wider frequency range; this concept will be detailed in Sect. 13.41 when the 
TDR/FDR combined approach is described. The lower the rise-time of the gener- 
ated signal, the more sophisticated (and costly) is the TDR instrument. 

It is worth mentioning that the minimum spatial resolution, AL m \ n , is related to 
the rise-time of the TDR step-like signal, ? r llfJl : 



AL n 



trC 



(3.9) 



Equation d3.8t is a simplification of /b = — 2 J, ■ which is an expression typically used 
in electrical engineering to describe the frequency characteristics of low-pass filters. This 
expression remains accurate when the energy contained in the voltage pulse is equally 
distributed across the frequency bandwidth and when there is no significant dispersion 
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Indeed, there is a wide availability of TDR units on the market: ranging from 
low-cost portable units to high-performance benchtop instrument. 

Most of the experimental results reported in this book were obtained using 
two TDR instruments with different characteristics, namely the Hyperlabs 1500 
(HL1500) and/or the Tektronix TDR80E04. 

The Hyperlabs 1500 is a low-cost portable TDR unit (see Fig. 13.5b |7]. The rise 
time of the generated signal is approximately 200 ps, which corresponds to a fre- 
quency bandwidth of around 1 .7 GHz. The incident pulse amplitude is 250 mV and 
the pulse duration is 14 (as. The maximum number of acquisition points, for the se- 
lected time window, is 2048. The HL1500 is a single-port instrument; the output has 
a BNC connector and the output impedance is 50 Q.. The HL1500 can be controlled 
through a computer. This TDR instrument supports multiplexer, thus allowing the 
simultaneous connection of several probes to a single TDR unit. This feature is 
particularly useful for practical purposes (e.g., for monitoring along the production 
lines, for geotechnical analysis of soils, etc.), as it would lower the implementation 
costs. 




Fig. 3.5 Picture of the 
HL1500 portable TDR in- 
strument |7J] 



As for the Tektronix TDR80E04, it is a highly sophisticated TDR module that 
is mounted on the Tektronix Digital Serial Analyzer DSA8200 (a benchtop instru- 
ment) shown in Fig. 13.61 The TDR incident step signal has a rise-time of approxi- 
mately 23.5 ps, thus the frequency content is about 15 GHz. The incident pulse am- 
plitude is ±250 mV and its duration is approximately 2.65 (xs. The maximum num- 
ber of acquisition points, for the selected time window, is 4000. The TDR80E04 has 
two ports (as it also allows transmission measurements to be performed). The output 
ports have 3.5 mm connectors. 

These two specific instruments were chosen, because they are representative 
of two categories of TDR instrumentation. The HL1500 is a compact and rugged 
TDR unit, and this make it particularly useful for in situ applications. On the other 
hand, the TDR80E04 is a high-performance instrument, and it is typically used in 
laboratories. 



3.3 Frequency Domain Reflectometry (FDR) 
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Fig. 3.6 Picture of the Tek- 
tronix DSA8200 equipped 
with the TDR80E04 unit d 




3.3 Frequency Domain Reflectometry (FDR) 

In FD, the parameter that corresponds directly to the time-domain reflection coeffi- 
cient, p(t), is the reflection scattering parameter, Su(f). In this case, the frequency- 
domain response function, R(f), can be written as 



R(f) = V (f)Sn(f) 



(3.10) 



where Vq (f) is the input signal B26l 13811 . Traditionally, measurements performed 
directly in FD rely on VNAs: in this case, the excitation stimulus is a sinusoidal 
signal whose frequency is swept over the desired range of analysis. 

VNAs usually have more than a single port; therefore, they can also be used for 
transmission measurements. A generalized block diagram of a VNA is shown in Fig. 
13.71 showing the major signal-processing sections. In order to measure the incident, 
reflected and transmitted signal, four sections are required yfl: 

1 . source for stimulus: most network analyzers have integrated, synthesized sources, 
providing excellent frequency resolution and stability. 

2. signal-separation block: this block is generally referred to as the 'test set' (which 
can be either integrated within or separate from the VNA). There are two func- 
tions that the signal- separation section must provide. The first is to measure a 
portion of the incident signal to provide a reference for ratioing (this can be done 
with splitters or directional couplers). The second function of the signal-splitting 
hardware is to separate the incident (forward) and reflected (reverse) traveling 
waves at the input of the SUT (for this task, bridges are often used). 

3. receivers that downconvert and detect the signals: there are two basic ways of 
providing signal detection. One relies on diode detectors which convert the RF 
signal level to a proportional DC level. If the stimulus signal is amplitude mod- 
ulated, the diode strips the RF carrier from the modulation (this is called AC 
detection). Diode detection is inherently scalar, as phase information of the RF 
carrier is lost. The other solution for detecting the signal relies on the tuned re- 
ceiver, which allows preserving phase information. The tuned receiver uses a 
local oscillator (LO) to mix the RF down to a lower intermediate frequency (IF). 



34 



3 Broadband Reflectometry: Theoretical Background 



incident 
signal 



stimulus 

signal 



DDE 



F^>4 



transmitted 
signal 



DUT 



> 



Sy Reflected 
I signal 



incident 
(R) 




Fig. 3.7 Generalized block diagram of a VNA [lj] 



The LO is either locked to the RF or the IF signal so that the receivers in the 
network analyzer are always tuned to the RF signal present at the input. The IF 
signal is bandpass filtered, which narrows the receiver bandwidth and greatly im- 
proves sensitivity and dynamic range. This architecture is called superheterodyne 
receiver, and it is used not only in VNAs, but, as well known, also in spectrum 
analyzers. 

processor/display for calculating and reviewing the results: this is where the re- 
flection and transmission data are formatted in ways that make it easy to interpret 
the measurement results. 



The picture of a VNA (model Agilent E8361C) is shown in Fig. 1X81 

As will be detailed in the following sections, one of the advantages of FD mea- 
surements is that there are well-established error correction models that, through 
appropriate calibration procedures [such as the short-open-load (SOL) calibration], 
can be used to reduce the influence of systematic errors |2|]. 

Although VNAs are still costly instruments, in |9[], a low-cost single-port VNA 
for high frequency measurements was developed. The major feature of the VNA 
proposed in ]9[] is that it eliminates the need for expensive heterodyne detection 
schemes required for tuned receiver-based VNA designs; as a result, the realization 
costs can be reduced dramatically. 
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Fig. 3.8 Picture of an Agi- 
lent E8361C VNA @|. 
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The FDR approach is often used for the dielectric spectroscopy of materials and 
for impedance characterization of electronic devices and components. In particular, 
a VNA used in conjunction with a dielectric probe kit can be employed to charac- 
terize the dielectric behavior of materials. However, it is important to point out that, 
although these are well-established experimental setups, their application is far from 
being universal, and often specific solutions must be designed. 



3.4 TD/FD Combined Approach 

The adoption of the TD/FD combined approach allows taking advantage of the ben- 
efits of both the approaches. In particular, estimating the frequency-dependent scat- 
tering parameter starting from TDR measurements, can help disclosing useful infor- 
mation that is masked in time domain (e.g., single or multiple dielectric relaxation) 
1 2111 . This strategy is regarded as a powerful tool for guaranteeing simultaneously 
low cost and portability of the instruments, and measurement accuracy. The Si\(f) 
of the monitored system can be calculated as the ratio between the discrete Fourier 
transformation (DFT) of the reflected signal and that of the input signal B22I1 : 



Sn(f) = 



DFT[r(t)} 
DFT [v (r )] 



(3.11) 



where vn(*) is the input function, which should be the unaltered signal generated by 
the TDR instrument and propagating (along the cable) up to the section where the 
5n(/) is being measured [29]. 

This approach provides reliability and adequate measurement accuracy, although 
relying on relatively-inexpensive TDR instrumentation. 

In this regard, it is worth mentioning that, unlike VNAs, portable and affordable 
TDR instruments (with wide operating frequency range) are commonly available on 
the market. For example, for the time being, the cost of portable TDR modules with 
a frequency bandwidth of 10 GHz and of 2 GHz is approximately 7000 USD and 
3500 USD, respectively. 
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3.4.1 Preserving Measurement Accuracy in the TDR/FDR 
Transformation 

There are some crucial aspects that need to be considered to perform an accurate 
TD/FD transformation, thus reducing errors on the assessment of the spectral re- 
sponse of the system and in the extraction of S\ \ (/) 12911 . 

First, the time duration of the acquisition window, T w , should be sufficiently long, 
so as to include all multiple reflections due to the signal traveling back and forth 
from the probe-end toward the generator. However, an excessively-long T w would 
lower the sampling frequency (since the maximum number of acquisition points for 
a given instrument is fixed), thus limiting the possible frequency range of analysis. 

Additionally, to reduce the spectral leakage caused by the rectangular truncation 
window, two main strategies can be adopted. One is the application of the Nicolson 
algorithm: a linear ramp is subtracted from the original TDR waveform, and each 
of the acquired N points is scaled according to the N' h final point; in this way, the 
value of the TDR waveform at the truncation point is zero. The resulting scaled 



signal, w* (nT c ), can be written as IB3I1 



«, s , n w(NT c )n 

W *(nT c ) = W {nT c )~ ^ c; (3.12) 

where T c is the sampling period, w(nT c ) is the rectangular- windowed signal, and 
n = 1 , 2, ..., N refers to each sampling point. The important feature of this algorithm 
is that the FFT of the Nicolson-modified signal produces the same response as the 



FFT of the original TDR signal [29:]. Alternatively to the Nicolson algorithm, it 
is possible to use the derivatives of the time domain signals calculated with the 
backward difference method; however, generally, noise is slightl y en hanced when 



the derivative of two similar quantities is taken in digital systems 124 12911 . 

Moreover, the addition of a number of zeros to the Nicolson-modified waveform 
(or to the derived signal) artificially enhances the frequency resolution of the ob- 
tained FD data, thus corresponding to an ideal interpolation in FD B291I . 
The steps at the basis of the TD/FD approach are summarized in Fig. 13.91 this proce- 
dure can be implemented, for example, in a MATLAB -based algorithm that, starting 
from the TDR waveforms, extrapolates the S\\ (/). 

In practical applications, the specification of the input signal, Vn(f), is crucial 
for obtaining accurate results. As aforementioned, the input function should be the 
unaltered signal generated by the TDR instrument and propagating (along the cable) 
up to the section where the scattering parameter is being measured 12911 : in fact, the 
Sn (/) should be evaluated at this very section, vrj(f) may be acquired by removing 
one of the electrodes from the probe [12411 ; however, this can be unpractical for most 
commercial probes, like the probe shown in Fig. 13 .41 in which the central conductor 
is attached to the molded probe-head. 

A possible alternative is the use an artificial TDR waveform modeled as an ana- 



lytical function [23]: 



1+ erf [<*(! -t )] 

V0,art(0 = ~ (3.13) 
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Fig. 3.9 Schematic diagram 
illustrating the algorithm 
for the TD/FD combined 
approach and for extrapolat- 
ing the reflection scattering 
parameter from TDR mea- 
surements 
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where erf is the error function, a t is the inverse of the rise time of the step-like 
signal, and ?o is the position where the input signal starts the rise. The parameter a t 
is evaluated by fitting the voart(0 to the TDR waveform of the probe in air I29n . 

As will be seen later, the adoption of the SOL calibration on FD-transformed 
TDR measurements can help circumventing the need of evaluating vq(V ). 



3.5 The Sensing Element 

The sensing element (or probe) is responsible for the interaction of the stimulus sig- 
nal with the SUT, and it is the ultimate factor that influences the accuracy of results. 
Since microwave reflectometry senses the changes in impedance, it is extremely im- 
portant to employ a probe with a well-known impedance profile; in this way, it is 
easier to discriminate and interpret the impedance variations due to the SUT. 

The high versatility of microwave reflectometry is also related to the possibility 
of customizing the probe configuration, thus employing an ad hoc solution for the 
specific needs. 

Coaxial probes, which are widely used for monitoring and diagnostics on liquids, 
are the most simple to design. A coaxial probe is composed of an outer cylinder 
(acting as the outer conductor) and a rod along the center line of the cylinder (acting 
as a central conductor). The impedance profile (in the TEM propagation mode) can 
be easily determined from the transmission line theory, as reported in Sect. 12. 11 121. 
11: 

Z(/) = -^=ln{b/a) (3.14) 

where Z(f) is the frequency-dependent impedance of the probe filled with the con- 
sidered material; £ r (/) is the relative dielectric permittivity of the material filling the 
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probe (which is assumed lossleso; b is the inner diameter of the outer conductor; 
and a is the outer diameter of the inner conductor. 

On the other hand, for soil monitoring and for granular materials monitoring 
in general, the multi-rod configuration is a widely used solution. In fact, this con- 
figuration allows an easier insertion of the probe into the granular material. For a 
two-wire transmission line, the characteristic impedance can be derived from (12.13b 
and it may be written as 

where D is the distance between the center of the conductors, and d is the diameter 
of each conductor. 

However, three-rod probes have become a widespread solution since their electric 



behavior resembles coaxial cells 14211 . Unfortunately, for configurations different 
from a coaxial line, models are not always available 112511 . and the evaluation of the 
probe impedance is not straightforward, especially for unbalanced probes. However, 
for multi-rod probes, an analytical formulation relating the impedance of the probe 



in air, Zp, to the geometric characteristics can be found in [12]: 

1 



Zp 



2iz{n-\) 
where H is given by 



H+^H 2 -l 



(3.15) 



( 2 2\n-l 2(»"1) 

H= /' a °> - °' (3.16) 

a ';- l {(s + ao)"- l -(s-a Q )>'-i} 

In the previous equations, n is the number of rods of the probe, a, is the radius of 
the probe center conductor, uq is the radius of the outer conductors, s is the distance 
between the middle of the center conductor and the middle of the outer conductor, 
£o = 8.854 x 10 _12 Fm _1 is the dielectric permittivity of free space,and jj.q = An x 
10 _7 Hm _1 is the magnetic permeability of free space. The cross-section of a four- 
wire probe is shown in Fig. l3.10l 



A simplification of (13.161 ) for three-rod probes is reported in 12711 : 



^p. three— rod 




l_ism (3 17) 



2(g/s 



where s is the center-to-center rod spacing, and g is the rod radius. 

Recently, the adoption of a single-rod probe for monitoring soil moisture was 
investigated. This probe exploits transverse magnetic (TM) propagation, rather than 
transverse electromagnetic (TEM) propagation. This kind of probe is simple to 



It is worth pointing out that when dielectric losses of the material cannot be neglected, 
equation ( 13. 14b can be extended by considering the complex dielectric permittivity. 
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Fig. 3.10 Cross section of a 
four-wire TDR probe 11211 
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realize and easy to insert; nevertheless, the analysis of the TDR waveform is some- 
what more complicated and the region of influence is relatively high B3411 . 

Fig. l3.11l shows the geometry and corresponding electrical fields of some typical 
prob e configurations. A three-rod probe and a two-rod probe are shown in Fig. 13.121 

in. 

It must be pointed out that, when a noninvasive approach must be preserved, it is 
possible to use surface probes, like the one shown in Fig. l3.13l l401. Alternatively, it 
is also possible to use antennas as sensing elements: this possibility will be discussed 

in Sect. E6| Hi- 
lt is worth mentioning that the probe configuration introduces some limitations 
on the useful frequency range of analysis, mostly because of the geometric char- 
acteristics. In particular, the TEM mode propagation assumption is true when the 
operating wavelength is higher than the dimensions of the probe 12611 . For example, 
for a coaxial probe, the circumferential resonance frequency / up . c irc (below which 
the TEM mode propagation assumption is valid) is given by 



:k 



2 c 



up.circ 



{a + b)n^/e v 



(3.18) 



where c is the speed of light and e r is the dielectric permittivity of the material filling 
the probqij- Similarly, the longitudinal resonant frequency is given by 12611 



/up, 



up.long 



2L^F r 



(3.19) 



where L is the length of the probe. 

As a matter of fact, several probe configurations are commonly available on the 
market; nevertheless, their configurations often do not allow an easy calibration 
procedure to be performed. As will be detailed later in the book, the calibration 
procedure consists in measuring the reflected signal when some standard (known) 



Equation d3.18t is equivalent to !2. 101 
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Fig. 3.11 Various probe designs (parallel plate and multi-wire configurations), typically used 
for TDR-based soil moisture. Electrical field lines generated for different probe configura- 
tions are also shown: closer line spacing is associated with a more concentrated field (i.e., 
greater influence on permittivity) 1 30] 



Fig. 3.12 Three-rod 
TDR probe (0.15 m- 
long; 3.20mm-diam. rods; 
12.0 mm rod spacing), and 
two-rod probe (Campbell 
Scientific CS616) fil 




impedances are connected at the section where the sensing portion of the probe 
begins. Fig. l3.14l clarifies this aspect; however, the importance of calibration will be 
discussed in detail in the following chapters. 
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Fig. 3.13 Serpentine- 
shaped probe 14011 
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Fig. 3.14 Schematic of 
the longitudinal section 
of a three-rod probe. The 
section where the reference 
standard impedances should 
be connected is indicated 



section where the 

calibration standards 
should be applied 




Indeed, the calibration procedure plays a major role in terms of minimization of 
systematic errors, thus leading to a significant enhancement of the final measure- 
ment accuracy. A viable alternative is to specifically design and realize the probe 
that would suit the specific needs [39] and to realize the corresponding calibration 
structure (as described for the application discussed in Sect. 14.3b . 



3.6 Strategies for Enhancing Accuracy in BMR Measurements 



As well known, measurements are typically affected by both systematic and random 
errors. The effects of systematic errors can be compensated for through appropriate 
calibration procedures 12811 . On the other hand, the residual random measurement 
fluctuations (essentially due the environmental and operating conditions, noise, etc.) 
are taken into account through a statistical analysis performed on a series of repeated 
measurements. 



3.6.1 Calibration Procedure in Time Domain 



In TD measurements, since the information on the SUT is retrieved from direct mea- 
surement of the apparent length of the probe immersed in the SUT, it is necessary to 
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Fig. 3.15 Signal flow graph 
describing the error model 
for the SOL calibration 
procedure 11511 



signal in . 9" 




locate exactly where the beginning and the end of the probe are. This is not a trivial 
task especially when unmatched probes are used. For example, for the widely-used 
three-rod probe schematized in Fig. 13.41 the actual beginning of the sensing element 
is section k = 1 and this point has to be individuated in the corresponding TDR 
waveform. In fact, when the probe is inserted in a material with specific dielectric 
characteristics, only L app will change. The offset length (Lq) has to be individuated: 
Lq corresponds to the electric distance between the cable end and the section k = 1 . 
The individuation of the point that, in the TDR waveform, corresponds to section 
k = 1 may be complicated because of the masking effect caused by the transition 
coaxial cable/rods and the portion of the rods included in the rigid support (i.e., 
probe-head). 

On such bases, to discriminate the electric distance of the rods inserted in the 
SUT, L a pp, from the electric distance of the probe-head, Lq, a simple but effective 
procedure has to be performed. To evaluate Lq, it is necessary to perform preliminary 
TDR measurements on materials with well-known dielectric characteristics, such as 



air and de-ionized water 1 16]. Considering Fig. 13.41 this can be simply summarized 
in the following equations: 



^p. water — M) ~r -^/ £app, water ^ 
^p.air M) ' V ^air-^ 



(3.20) 
(3.21) 



where L is the physical length of the sensitive portion of the probe in air; £ app . water — 
80 is the apparent permittivity of water; and £ a ; r = 1 is the dielectric constant of air. 
Once Lq is individuated, it can be compensated for in successive measurements. 



3.6.2 Calibration Procedure in Frequency Domain 



Frequency domain reflectometry can rely on well-established calibration proce- 
dures. For single-port measurements, such as microwave reflectometry measure- 
ments, the short-open-load (SOL) calibration is largely used 11411 . The error model 
can be identified using a signal flow graph as shown in Fig. l3.15l Si i, sol is the re- 
sponse of the sample under test; Sn,M is the measured response (which includes the 
parasitic effects); and a, b and c are the parameters calculated from the open, short 
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and 50Q-load calibration measurements. The SOL-calibrated frequency response 
of the sample under test can be extracted from the signal flow graph according to 

S] 1 M — O- 

Sii,sol - . c ' r- (3.22) 

where a, b, and c include the contribution of Sh.mo, ■S'ii.ms. an d ^ii.ml, which are 
the reflection scattering parameters when an open, a short and a matched loads are 
connected: 

a = Su ML , (3.23) 

,_ Sll,MO + Sll,MS ~ ^iSj^ML C\JA\ 



■Sll.MO — 5n,MS 

2 (Si i ,ML ~ Sj i ,Mp) (Si 1 .MS ~ 25i i ,ML 
Sll,MO~ ^lKMS 



(3.25) 



The aforementioned error-compensation model is solved directly in FD; neverthe- 
less, it can implemented in the TD/FD transformation algorithm described in 11511 . 



thus extracting the calibrated Si i (/) from TDR measurements [17]. The aforemen- 
tioned calibration procedures were implemented in the transformation algorithm 
used for transforming TD data into FD data. An additional advantage of this strat- 
egy is that it overcomes the need of acquiring the input function, vrj. 
However, calibration standards are often not readily available for many commercial 
probes. In these cases, it is necessary to fabricate appropriate calibration fixtures or 
custom-made calibration standards. 

An innovative calibration procedure, namely triple short calibration (TSC), has 



been presented in 13511 : this procedure, which will be described in detail in Chap. [51 
relies on three short-circuit applied at three different sections along the probe. The 
major advantage is that realizing a short-circuit condition is easier than achieving 
an open-circuit condition or a matched condition, especially for probes as the one 
represented in Fig. 13.41 



3. 6.3 Time-gated Frequency Domain Approach 

Time gating is a calibration strategy that can be adopted for excluding systematic 
errors in the TD/FD combined approach |5J. Time gating implies the individuation 
and the removal of those portions of the waveform that are considered undesired 
effects (e.g., reflections caused by connection transitions): once these effects are 
removed from the TD data (via software), the data are transformed into the FD. As 
a result of the process, the gated frequency response will resemble more closely the 
response of the sole SUT. As a matter of fact, the time-gating option is often directly 
integrated within VNAs: in this case, the FD response is transformed in TD; then 
gating is applied, and, finally, the data are transformed back into the FD. 
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Fig. 3.16 Schematic dia- 
gram illustrating the imple- 
mented algorithm for the 
TD/FD combined approach 
and for extrapolating the 
calibrated reflection scat- 
tering parameter from TDR 
measurements 
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As an example, Fig. l3.17l shows the sequence of the gating operation for isolating 
the frequency-domain response of the outer connector of a cable. Considering the 
non-gated time domain response of the cable, the discontinuity on the left is due to 
the input connector and the discontinuity on the right is due to the output connector. 
By using the gating function, the effect of the mismatch can be removed in order to 
see the frequency response of just the output connector; in this way, the frequency 
response of the sole output connector is obtained. 

It must be pointed out, however, that gating is an extremely delicate strategy, and 
it should be adopted with particular care. 

Time-gated frequency-domain reflectometry can be particularly useful for removing 
undesired reflections from the environment while measuring the Si i (/) of antennas 
when the anechoic chamber is not available [19]; this aspect will be further dis- 
cussed in the following chapters. 



3. 6.4 Transmission Line Modeling and Inverse Modeling 



An additional strategy for compensating for systematic errors is to model the influ- 
ence of the parasitics through a transmission line (TL) model. More specifically, the 
SUT is schematized as a circuit model, and the parasitics are modeled as lumped cir- 
cuit elements. This approach can virtually compensate for all the systematic errors 
effects. 
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Fig. 3.17 Sequence of the gating operation, a frequency response of a cable, b time domain 
response, c gating applied to remove the first discontinuity, d frequency response with gating 
on and off [5] 



This approach is based on the fact that the Sn(f) at a specific section is to be 
seen as the resultant of infinite contributions, that can be mathematically repre- 
sented through the multiple reflection theory 02OI1 . The measurement cell can be 
considered as a multi-section transmission line, and M is the number of sections 
where impedance changes occur. Let k = 0, ...,M be the number of sections starting 
from the end of the probe (k = 0) up to the TDR output (k = M). For the sake of 
example, Fig. l3.18l shows the 'main' sections of a three rod probe. 



instrument 
connector 



probe-head 



reflectometer 




Fig. 3.18 Schematization of the 'main' sections of a typical three-rod probe 
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The scattering parameter at the gene ric section k, due to the sections up to the ]i h 
is given by the following equation J2Cll23L 13811 : 



c*,^ r s k {f)+s k n \f)e-^ 



Sll(/) ~i+w)sV(/)'- 2 * °' 26) 

where S]~[ (f) is the scattering parameter at the (k— l) rh section due to the first 
(k — 2) sections; r s k (f) is the reflection coefficient; y is the complex frequency- 



dependent propagation parameter [23]; and L k is the length of the U section. The 
losses due to relaxation and conductivity for a TEM mode propagation, y(/) L are 
given by 112511 

y(f)L = i2nfL^Lill ( 3.27) 

c 

However, an assumption must be made about the reflection at the end of the probe 
(i.e., section k = 0): for an open-ended probe, Sj j (/) = 1 ; for a short-circuited probe 
S° n (/) = - 1 ; and for a matched probe S° n (/) = 0. 

The reflection coefficient between the different sections k of the transmission line 
is 

Fs (/) = z k ^f) + z k (f) (3 ' 28) 

where Z k (f) is the impedance of the transmission line at section k and is given by 

Z k {f) = -?M= (3.29) 

Zo k being the characteristic impedance of the transmission line at section k. 
Equation (13.26b is recursive; therefore, to evaluate Si i (/) for the whole system (i.e., 
at the M th section), it is necessary to subsequently compute S k n (f) for k= 1, for 
k = 2, ... ,k = Af(ffl. 

The practical usefulness of the TL modeling will be discussed thoroughly in the 
following chapters. However, it is useful to anticipate one major application, which 
is the extraction of the Cole-Cole parameters of materials. The approach proposed 
for the accurate evaluation of the Cole-Cole parameters is the combination of mul- 
tiple strategies, and it comprises three major steps: 

• First, the calibrated Si [ (/) is evaluated (either through direct VNA measurements 
or through the TD/FD combined approach); 

• Secondly, the used measurement cell is accurately modeled as a TL, and the 
considered material is parameterized through the Cole-Cole model; 

• Finally, through the minimization of the deviations between the measured Si i (f) 
and the reflection scattering parameter of the implemented TL model, the Cole- 
Cole parameters of the material are obtained. 

These steps are summarized in the diagram in Fig. 13. 191 and will be described in 
detail in the following subsections. 
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Fig. 3.19 Schematization of the TL-based procedure for extracting the Cole-Cole parameters. 
The activities enclosed in the gray-colored area can be substituted by direct VNA measure- 
ments of the Si i (/) 
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Chapter 4 

Quantitative and Qualitative Characterization 

of Liquid Materials 



'Science exists, moreover, only as a journey toward troth. 
Stifle dissent and you end that journey.' 
John Charles Polanyi 



Abstract. In this chapter, broadband microwave reflectometry (BMR)-based meth- 
ods for measuring simultaneously the levels and the dielectric characteristics of liq- 
uid materials are presented. 

First, an approach based solely on time domain reflectometry (TDR) measure- 
ments and on the analysis of TDR waveforms is described. 

Secondly, a further enhancement of this method is accomplished by resorting 
to the combination of TDR with frequency domain reflectometry (FDR). Also this 
method provides, in one shot, the level and the dielectric permittivity of the con- 
sidered liquids (this time, the complete spectral signature is obtained in terms of 
Cole-Cole parameters). This approach is considered the optimal solution especially 
when dealing with dispersive media. In this second measurement method, the ac- 
curacy is enhanced also through the adoption of a transmission line (TL) modeling 
of the measurement cell and through the realization of a custom-made fixture that 
allows performing short-open-load (SOL) calibration measurements. 

Finally, BMR-based measurements of the dielectric characteristics are extended 
to edible liquids: in particular, vegetable oils are considered. 

Reported results show that all the proposed methods have strong potential for 
possible practical implementation in the field of industrial monitoring. 

It will be evident how, in all the considered approaches, the design of the probe 
is crucial for enhancing the final measurement accuracy. 



4.1 Introduction 

In recent years, dielectric spectroscopy measurements performed through broad- 
band microwave reflectometry (BMR) have become an effective tool for extracting 
information on the qualitative and quantitative status of materials, thus opening in- 
teresting perspectives for several monitoring applications. 
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With regards to the fluid processing-related industry, measurement and control 
of liquid parameters play an important role in industrial applications, for example, 
in the petrochemical industrial processing. In this field, additional stringent require- 
ments frequently deal with the spatial localization of non-miscible-layered liquids 
and with the depth estimation of different layers of liquid. As an example, during 
the refinery industrial processing, a portion of water is often present in the bottom 
of the tanks, as a non-miscible fraction, and its continuous monitoring is crucial. 
For this purpose, typically, various flow meters are used for fiscal and quality pur- 
poses; nevertheless, the a-priori knowledge of the permittivity of the involved fluid 
components, and the associated sample calibration are mandatory for the practical 
implementation of such methods. 

Starting from these considerations, in this chapter, first, a time domain-based 
approach for the simultaneous estimation of the dielectric characteristics of liquid 
materials and of their level is presented. It is shown that the analysis of time domain 
reflectometry (TDR) data allows, in one shot, the measurement of liquid levels, 
the determination of multiple interfaces in layered media, as well as the evalua- 
tion of dielectric properties. Experimental results demonstrate the robustness and 



reliability of the proposed method. As reported in 12 111 , which compares the perfor- 
mance of different liquid level sensors, TDR fares well with respect to other sensing 
technologies. 

Secondly, an enhancement of the simultaneous measurements of liquid levels and 
dielectric permittivity is achieved through the implementation of the TD/FD com- 
bined approach described in subsection l3.6.4l The proposed combined approach al- 
lows the evaluation of the Cole-Cole parameters of the considered samples. Indeed, 
accuracy enhancements are obtained also through some novel strategies that include 
i) an appropriate probe design; ii) the realization and the adoption of a custom-made 
fixture for allowing calibration measurements to be performed; and, most impor- 
tantly, iii) a targeted optimization routine that renders the method extremely robust. 
Results deriving from simple TDR technique are compared with those achieved, in 
frequency domain, both through the fast Fourier transform (FFT)-based processing 
of TDR data and through direct measurements with vector network analyzer (VNA). 
The comparison between the reference and the frequency-domain extrapolated val- 
ues of permittivity and levels shows a very good agreement. 

Finally, in this chapter, the estimation of the Cole-Cole parameters through BMR 
is applied to edible liquids: in particular, an interesting test-case for vegetable oils 
is presented. Also in this case, starting from TDR measurements, the proposed ap- 
proach allows the accurate estimation of the Cole-Cole parameters of the considered 
oils, thus providing a method useful for qualitative and anti-adulteration control pur- 
poses (which are regarded as crucial issues in the specific sector). 

The common important feature of the proposed approaches is the possibility 
of automating the necessary data-processing, thus providing results in real time. 
On such bases, the proposed procedures can be considered as excellent candi- 
dates for qualitative and quantitative liquid monitoring applications, with relevant 
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perspectives in the field of industrial monitoring, which not only does require accu- 
racy and reliability, but also low cost and portability of the involved devices. 

4.2 Estimation of Levels and Permittivities of Industrial 
Liquids Directly from TDR 

In this section, a TDR-based method for simultaneous quantitative and qualitative 
monitoring of liquids is described. As outlined in the previous chapters, the level of 
liquids in tanks can be inferred from the travel time of the electromagnetic signal 
along the probe. The basic idea concerns the use of a perfectly-matched coaxial 
probe immersed in a lossless non-dispersive liquid. Under these circumstances, the 
unknown dielectric characteristics can be evaluated from (13.14b . as 

where p is the reflection coefficient in time domain. 

The quantitative measurement of liquid level is based on the spatial localization 
of different interfaces. This task can be performed simply by considering the reflec- 
tion coefficient (or, even better, its derivative) versus the apparent distance. In fact, 
the clear peaks occurring at significant discontinuities are associated with a local 
maximum (or minimum) of the reflection coefficient derivative curve; hence, they 
spatially localize the points associated with different interfaces. For example, in the 
TDR waveform of the system schematized in Fig. 14. II a strong drop of the reflection 
coefficient would be clearly detectable at the air-liquid interface (due to the different 
dielectric constant of air and of the liquid). 

Once the dielectric characteristics of the liquid are known, the level of the liq- 
uid (or height, H) can be evaluated directly from ( 13.5b . which is rewritten here for 
convenience: 

h = -^ee- (4.2) 



where H is the physical length of the probe portion immersed in the liquid under 
test (LUT), and // app is the corresponding apparent distance. 

To summarize, the adopted measurement method comprises the following steps: 

1 . first, the apparent liquid level (// a pp) is evaluated from the derivative of the mea- 
sured TDR waveform; 

2. secondly, using d4.lt . the measured reflection coefficient value is used to calcu- 
late the apparent dielectric constant value of the considered medium; finally, 

3. equation ( 14.2b allows retrieving the liquid level, H. 

In spite of the simplicity of the listed steps, there are some important aspects that 
must be considered to obtain accurate results: these aspects are thoroughly discussed 
in subsection l4.2.2l 
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Fig. 4.1 Setup for TDR-based liquid level measurements: the probe is immersed in the liquid, 
and it is connected to the TDR instrument (which, in turn, is controlled through a computer) 



4.2.1 Probe Design and Realization 

As aforementioned, coaxial probes are the wisest choice for applications on liquids. 
This probe configuration can guarantee the insertion into the liquid samples, the me- 
chanical stability, the electromagnetic continuity with the connection to the coaxial 
cable, as well as an adequate impedance matching. This strongly enhances the over- 
all performance of the measurement method. The first aspect to take into account 
is the minimization of the impedance mismatches along the line (i.e., TDR unit, 
coaxial cable and probe in air), so that ( 14. H can be considered reasonably valid. 

The schematic configuration of the used probe is shown in Fig. 14.21 The designed 
coaxial probe is made up of a central cylindrical conductor and of a coaxial stain- 
less steel conductive shield. The values of the outer diameter of the inner conductor 
(a) is 10 mm, whereas the inner diameter of the hollow outer conductor (b) is 23 
mm: these dimensions allow a good 50 Q-matching, according to J3.14I ). The outer 
conductor is perforated, so as to allow fluid circulation when the coaxial probe is im- 
mersed in the liquid (four series of holes with a diameter of 8 mm and 3 cm- spaced, 
are made longitudinally along the external hollow conductor). The 1 cm-thick metal- 
lic flange facilitates the positioning of the probe in the tank/sample holder (hence, it 
can be customized in order to fit the required operating conditions). The longitudinal 
cross-section scheme, reported in Fig. 14.21 clearly shows the other design details. 
The probe-head and the outer conductor are 6 cm and 43.5 cm long, respectively. 
The inner conductor is centered on the coaxial probe-head, through a 2 cm-teflon 
ring, while the lower part of the coaxial probe (probe-end point) is open-ended. 

The adopted 50 Q-matched coaxial probe shows very good performance in terms 
of impedance mismatching minimization. This can be seen from Fig. 14.31 which 
shows that the measured reflection coefficient associated to the probe in air is ap- 
proximately zero. 
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Fig. 4.2 Coaxial probe picture, dimensional scheme of cross-section, and longitudinal geo- 
metric scheme 
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Fig. 4.3 TDR waveform and derivative of data for the coaxial probe in air 
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4.2.2 Sources of Errors and Compensation Strategies in TD 

The overall contribution to the uncertainty in TDR measurements can be associated 
to some error sources that are related to the to the specific measurement setup (and, 
consequently, to the electromagnetic characteristics of the signal propagating along 
the transmission line). Once the error components are characterized, the related ef- 
fects can be de-embedded from measurements. 

The overall uncertainty affecting TDR measurements can be mostly attributed to 
the following causes: 

1. the undesired impedance variations introduced by the non-ideal cable-to-probe 
connection, by the teflon ring inserted in the probe-head, etc.; 

2. the non-ideal 50 Q-matching of the coaxial cable and of the probe; 

3. the effects of losses and signal dissipation in the coaxial cables (which increase 
as the cable length increases); 

4. the assumptions of lossless and non-dispersive liquid; and, finally, 

5. the periodical presence of parasitic multiple reflections in TDR waveforms, 
caused by the signal traveling back and forth along the transmission line lit]. 

Although the errors sources listed above refer to the specific setup, they can be 
present in most BMR-based measurement systems. 

The aspects listed in points l)-3) can be attributed to the non-ideal characteristics 
of TDR equipment and the related effects mainly cause systematic errors, both in the 
estimation of the apparent distances and in the evaluation of the reflection coefficient 
values. On the other hand, the aspects listed in points 4) and 5), are mainly related 
to the intrinsic limitation of the method itself, which can circumvented through the 
approach described in Sect. 14.31 

4.2.2.1 Minimization of Errors in the Evaluation of Apparent Distance 

Error sources described in point 1) basically lead to an overestimation of the ap- 
parent distances; hence, they can be considered systematic. These effects can be 
minimized by calibration measurements on 'reference distances'. 

As seen in subsection l3.6.1l it is important to evaluate Lq (which, as reported in 
Fig.s l4.2l and l4.3l refers to the probe-head portion of the probe) through preliminary 
calibration measurements (in air and in de-ionized water), thus individuating the 
portion of the waveform L app that refers to the effective sensitive portion of the 
probe, L. 

4.2.2.2 Minimization of Errors in the Evaluation of the Amplitude of the 
Reflection Coefficient 

The parasitic effects described in points 2) and 3) are mainly related to an atten- 
uation and dissipation of the electromagnetic signal, which can be due both to the 
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cables and to the material under investigation itself. The related error contribution 
basically affects the dielectric constant estimation described by J4.lt . However, a 
separate calibration procedure can help discriminating the attenuations due to the 
sample under investigation from that due to the measurement setup. In fact, the 
coaxial cable loss attenuation constant, A, can be found through the following pro- 
cedure ]4j]: measuring the reflection coefficient when the distal end of the probe is 
open-circuited and short-circuited (i.e., p^ c eas and p^eas)' an ^ comparing the results 
with the theoretical values (i.e., Py c eal = +1 and Py eal = —1), the reflection coeffi- 
cient amplitudes can be corrected, according to 

A = exp(2aL cable ) (4.3) 

Pcorr = ^Pmeas (.4.4) 

where a [m~ ' ] is the attenuation coefficient of the coaxial cable, and L ca t,i e [m] is the 
length of the coaxial cable. The constant A accounts for signal losses over a path of 
length 2L ca bie, and equals 1 for ideal cables. For the specific 3 m-long coaxial cable 
used in the measurements reported in this section, the evaluated A is 1 .06. 

4.2.2.3 Intrinsic Limitation of the Proposed Method 

Finally, the implications of points 4) and 5) can be considered as intrinsic limi- 
tations of the proposed method itself and of the considered measurement setup, 
respectively. 

Generally, the presence of a steady curve portion in the TDR waveform indicates 
that the liquid can be considered as a lossless dielectric. Conversely, when dissipa- 
tive liquids are considered (such as electrolyte solutions or water-like liquids), the 
ohmic and polarization losses cause a strong variation of the reflection coefficient at 
the probe-portion filled with liquid. This problem can be circumvented through the 
estimation of the static electrical conductivity, or through the consideration of the 
FD approach (as detailed in the following section). 

As for the parasitic effects of multiple reflections described in point 5), they can 
be caused, for instance, by the connection between the TDR unit and the cable, by 
the cable-to-probe connection, and by the plastic ring insertion in the probe head. 
They cause a periodical presence of parasitic multiple reflections of TDR wave- 
forms. For example, in Fig. 14.31 it can be noticed that the effects of parasitic mul- 
tiple reflections are clearly detectable at longer distance after the probe-end point 
(approximately 4.9 m). Considering the distance where it occurs, this effect is at- 
tributable to the repetition of the reflection occurring at the probe head section; 
however, this does not play a dominant role in terms of error contribution for the 
considered application. 
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4.2.3 Validation of the Method on Dielectric Liquids 

The TDR-based method was tested on four liquids (with different levels), which 
were representative of 'high' (de-ionized water), 'medium' (acetone), and 'low' 
(diesel oil and fuel) permittivity values. These liquids have negligible ionic con- 
ductivity and relaxation effects in the considered frequency range (0-1.7 GHz). The 
liquids were placed inside a cylindrical sample holder, and their level (// re f) was 
measured directly from the reading of the lOOOml-graduated sample holder. 

TDR measurements were performed through the HL1500; the number of mea- 
surement averages for the acquisition was 10, and the number of acquired points 
was 2048. The TDR waveforms for the considered liquids are shown in Fig. 14.41 
The significant reflections, clearly shown by the derivative curves, correspond to 
coaxial cable-to-probe connection, to the air-liquid interface, and to the probe-end 
point. The distance differences between the probe-end peaks and the air-liquid in- 
terface peaks (located through the analysis of derivative data for each sample), give 
the apparent liquid lengths. 

It is important to point out that, for evaluating the dielectric constant from ( 14. It . 
the reflection coefficient value should be retrieved from the steady- state portion of 
the TDR waveform, as long as the low-loss non-dispersive liquid assumption can be 
considered valid. 
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Fig. 4.4 TDR waveform (black line with dots) and derivative of data (gray 
1 ine). a for a 41 .3 cm level of de-ionized water sample [probe completely full], b for 27.2 cm 
level of acetone sample, c for 32.1 cm level of fuel sample, d for 27.8 cm level of diesel oil 
sample 
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For each liquid the following data were evaluated: 

1 . the TDR-measured apparent liquid level, // app ; 

2. the measured reflection coefficient, p mea s; 

3. the corrected reflection coefficient according to ( 14.4b . p CO n-; 

4. the square root of the calculated dielectric constant y/e^', and, finally, 

5. the actual liquid level, H meas . 

To provide a simple performance assessment of the method, the relative percentage 
difference between reference and measured liquid levels was calculated: 

AH = \H m{ - # meas 1 100/// ref . (4.5) 



The results for level and dielectric constant estimation are summarized in Table l4.ll 
it is worth noting that the experimental results for the level measurements show a 
very good agreement with the reference levels; similarly, the measured dielectric 
constant values agree with those reported in IU8I1 . 



Table 4.1 Summarized results obtained from TDR measurements for the estimation of levels 
and dielectric constants of liquid samples reported in Fig. l4.4l (de-ionized water, acetone, fuel, 
and diesel oil) 



#ref 

(cm) 


probe end 

point 

(cm) 


probe-to- 
interface 
(cm) 


■liquid 
point 


^*app Pmeas 

(m) 


pcorr 


v/Eapp 


^meas 

(cm) 


AH 
(%) 


41.3 


776.5 


401.0 




deionized water 
375.5 -0.75 


-0.80 


8.97 


41.8 


1.2 


27.2 


546.5 


416.7 




acetone 
129.8 -0.61 


-0.65 


4.71 


27.6 


1.5 


32.1 


459.0 


411.0 




fuel 
48.0 -0.18 


-0.19 


1.47 


32.6 


1.5 


27.8 


457.2 


415.0 




diesel oil 
42.7 -0.19 


-0.20 


1.50 


28.1 


1.1 



4.2.3.1 Measurements on Stratified Liquids 

The proposed measurement system was also validated for measurements of the lev- 
els of stratified liquids (which can be particularly useful in the field of petrochem- 
ical, chemical, or agro-food, industrial processing control purposes). In this case, 
measurements were performed with the probe-end termination short-circuited. 

A 27.5 cm-high diesel oil layer was considered on a 9.5 cm-high tap water 
sample: the corresponding TDR waveform is reported in Fig. 14.51 The measured 
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5 6 7 

apparent distance (m) 

Fig. 4.5 TDR waveform and derivative of data for a diesel oil-tap water layered sample 
(short-circuited probe-termination) 



apparent height for tap water, //™ ter , was 81.1cm, and the corrected reflection 
coefficient was -0.8, corresponding to the expected value of dielectric constant 
(£water — 80). Also the third reflection related to diesel oil-water interface could be 
directly evaluated as difference between the total coaxial probe length, the coaxial 
probe in air portion and the diesel oil distance. 



4.3 Estimation of Levels and Permittivities of Industrial 
Liquids Using the TD/FD Approach 



The simple TD approach described in the previous section provides adequately ac- 
curate results as long as the investigated liquids can be considered non-dispersive 
and lossless. Conversely, when the dependence of permittivity on frequency or the 
effect of losses cannot be neglected, to enhance the measurement accuracy, specific 
dielectric models (such as the Cole-Cole formula) have to be adopted for describing 
permittivity. 

This can be achieved adopting a frequency-domain approach, thus considering 
the reflection scattering parameter, S\\(f). As mentioned in Sect. 13. 41 this parameter 
can be obtained either through direct frequency-domain measurements with a VNA 
or through an appropriate fast Fourier transform (FFT)-based processing of TDR 
data Jl MM, El- 

On such bases, the approach proposed herein for the simultaneous estimation 
of liquid levels and of permittivity (in terms of Cole-Cole parameters), relies on 
FD data and on the transmission line (TL) modeling of the measurement cell (i.e., 
probe + LUT). In particular, the measurement cell is modeled as a function of the 



4.3 Estimation of Levels and Permittivities of Industrial Liquids 



61 



'unknown' Cole-Cole parameters of the liquid filling the probe and of the 'unknown' 
liquid levels. Through the minimization (over the unknown variables) between the 
measured scattering parameter (Sh.meas) and the scattering parameter of the TL 
model (Si i. mod), the proposed procedure provides the accurate evaluation of the 
Cole-Cole parameters and of liquid levels. Furthermore, this procedure can be ex- 
tended to stratified liquids 12211 . 

As will be detailed later, an additional accuracy enhancement was obtained by 
realizing an ad-hoc fixture that would allow performing short-open-load (SOL) cal- 
ibration measurements, which are essential to reduce the effect of systematic errors 
(as seen in subsection l3.6.2l i. 

Experimental measurements were performed considering different levels of lay- 
ered samples (fuel/tap water). Since both fuel and tap water exhibit a strong dis- 
persive behavior (especially, the imaginary part of their complex permittivity), a 
complete spectral picture of the frequency-dependent dielectric properties can be 
obtained through the five-parameter Cole-Cole formula. 



4.3.1 Transmission Line Modeling of the Measurement Cell 

The approach proposed herein is based on the adoption of a TL model of the 
measurement cell. As aforementioned, the complex relative permittivity of fuel and 
water, £*, can be appropriately modeled through the Cole-Cole formula. Such a 
function, both for fuel and water, is incorporated in the S\ i (/) multiple-reflection 
model described by (13.26b . through (l3.27b -( [3T29l >. The TL model is parameterized 
also as a function of the unknown liquid levels. 

On such bases, as reported in in subsection 13.6.41 the minimization (over the 
unknown variables) of the squared difference between the measured, Sh.meas(/). 
and the modeled scattering parameters, Sii,mod(/)> leads to the estimation of 
the Cole-Cole parameters and of the corresponding liquid levels. For the adopted 
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measurement setup, the considered multi-section TL model comprises three sec- 
tions: section is the probe termination, section 1 corresponds to the probe portion 
filled with water, and section 2 corresponds to the portion filled with fuel. 

To estimate the input scattering parameter S\f 2 (f) of this multi-section TL 
model, the value of the reflection coefficient at the termination has to be assumed. 
The used probe has an open-ended configuration and, therefore, S\^° = 1 is as- 
sumed. Starting from this value and using (13.26b . firstly for section 1 and then for 
section 2, the input scattering parameter can be obtained. It is worth noting that 
T k=l = 1, due to the open-ended termination, while Zq^ = 50 Q for all the remain- 
ing sections (due to the specific probe design). 

The proposed model was implemented within the commercial microwave circuit 
simulator AWR Microwave Office (MWO), which also makes optimization routines 
directly available. In particular, the chosen objective function to be minimized is the 
sum of the squared magnitude of the differences between the measured and the 
MWO-modeled scattering parameters. 

This way, since the Cole-Cole formula is included in the modeled scattering pa- 
rameter, along with the liquid level, the minimization of the objective function suc- 
cessfully provides an accurate estimation of the five dielectric parameters (over the 
considered frequency range) and of the levels of liquids. 

It must be pointed out that the optimization is rather complex because of the 
large number of parameters involved: the fuel level (the water level is constrained 
to the difference between the known probe length and the fuel level), four Cole- 
Cole parameters for the fuel (the static conductivity is constrained to S m~ : ), and 
five Cole-Cole parameters for the water. The optimum Cole-Cole parameters and 
levels for fuel and water are evaluated through the following three- step optimization 
procedure: 

1. the Cole-Cole parameters for water and fuel are fixed at reasonable initial guess 
values (the dielectric properties of water are well known once the temperature is 
measured), thus the fuel level is optimized. This gives a reasonable initial guess 
for the fuel level; 

2. starting from the previous values, all parameters are optimized through a random 
optimization procedure that should take the objective function close to the global 
minimum and should guarantee that the subsequent optimization step does not 
fall into a local minimum; 

3. using the previous results as the new starting inputs, a simplex optimization 
procedure is applied, which forces the objective function to the final global 
minimum. 

The robustness and efficiency of the proposed procedure will be shown in the next 
subsections. 



4.3.2 Measurement Setup and Enhanced Probe Design 

As aforementioned, the probe configuration design plays a crucial role in terms of 
measurement accuracy. In fact, any impedance mismatch in the TL may affect the 
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evaluation of the liquid dielectric properties and levels. To circumvent this problem, 
the optimal strategy should simultaneously involve a suitable calibration procedure 
and a good impedance matching of each component of the measurement setup. 

In this regard, the most critical component to be designed is the probe-head sec- 
tion, which should guarantee a perfectly matched cable-to-probe connection, as well 
as the mechanical stability of the probe central conductor. Additionally, the probe- 
head structure should be included in the SOL calibration measurements, in order 
to de-embed its effect and to accurately derive the response related only to the 
probe section filled with the sample under test. As a matter of fact, all these re- 
quirements were considered in the design of the the enhanced configuration of the 
coaxial probe. Fig. 14.71 shows the 433 mm long coaxial probe, the 160 mm long 
probe-head portion as well as the calibration structure that exactly reproduces the 
internal probe-head portion: all these components are virtually 50 Q-matched, as 
explained in the following. For this probe, the inner and outer diameters of the probe 
conductors (a = 9 mm and b = 21 mm, respectively) were chosen, once again, to en- 
sure a 50 Q-matched characteristic impedance when the probe is in air, according 
to (El . 

The chosen dimensions limit the frequency bandwidth of the probe to about 
6 GHz, when the dielectric material considered is air. However, when the probe is in- 
serted into a different material, this limit decreases proportionally to the square root 
of the dielectric permittivity of the material. In particular, for the considered cases, 
the liquid with the highest permittivity is water, which limits the useful bandwidth 
of the probe to approximately 700 MHz. 




open*ended 
probe-end \ 



433 mm-long 
Cgagjgj probe 



50 ii calibration 
astrucwre willi 
double sided N- 
conneclors 




fc=21 mm 



perforation 



Fig. 4.7 Coaxial probe external configuration, details of calibration structure, internal probe 
head and inner probe conductor, external probe head shield and external probe shield, cross- 
section probe dimensions II22I1 
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The probe-head internal portion is made up of a 50 Q foam-dielectric coaxial 
cable (Cellflex LCF78-50JA), contained in the stainless steel outer conductor that 
is ground contacted. The 433 mm-long inner conductor of the probe is obtained 
removing the same length of the external shield and foam dielectric of the Cellflex 
cable. 

As aforementioned, to perform a SOL calibration at the exact section where the 
sensing element begins (section #2 of Fig. 14.61 ). a calibration structure reproducing 
the internal probe-head portion was specifically realized. This structure terminated 
with an N-type connector (thus ensuring the connection to N-type SOL standards). 

For the time-domain measurements, the HL1500 unit was used. Comparative fre- 
quency domain measurements were performed through the Agilent E8363B VNA, 
whose upper operating frequency is 40 GHz. The probe was connected to the TDR 
and to the VNA instrument through a 3 m-long RG213 coaxial cable. 

The HP85032B N-type SOL standard calibration kit was used for calibration 
measurements, both for VNA and TDR measurements. In particular, the SOL kit 
standards were connected at the end of the 'calibration structure' (which was con- 
nected to the instrument through the 3 m-long cable as well). In this way, all the 
parasitic effects caused by the probe-head (and, consequently, also of the cable) 
were adequately compensated for. 

4.3.3 Experimental Validation of the Method 

To assess the method, measurements were performed on stratified liquids, namely 
fuel/tap water layered samples. Three samples were considered. All measurements 
were performed at room temperature (26±1°C). The three samples differed in the 
relative level of fuel and water, as shown in Table 14.21 First, the apparent dielectric 

Table 4.2 Reference level of fuel (H^ ) and tap water (H^), for the three considered samples 





(mm) 


zjref 

(mm) 


Sample #1 


283.0±0.5 


150.0±0.5 


Sample #2 


219.0±0.5 


214.0±0.5 


Sample #3 


335.0±0.5 


98.0±0.5 



constant and the levels were estimated through the time-domain approach presented 
in the previous section. 

Secondly, the dielectric permittivity and the levels were evaluated by resorting 
to the proposed TL modeling and optimization procedure. For comparison, this 
TL-based approach was applied separately to two sets of data: to the scattering 
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parameter derived directly from VNA measurements (5^" a MEAS ) and to the scattering 
parameter derived from the FFT of TDR data (Sf[ MEA s)- 

In this way, it was demonstrated that the measurement accuracy in the evalua- 
tion of levels and dielectric permittivity is substantially improved with respect to 
the simple time-domain approach. Additionally, it was also shown that the TD/FD 
combined approach can provide as much accuracy as direct VNA measurements. 

Finally, it is important to highlight that, for all the TDR and VNA measurements, 
the SOL calibration procedure (with the calibration structure connected at the end of 
the feeding coaxial cable) was performed, according to the calibration error model 
reported in subsection l3.6.2| [51. Additionally, also the TDR measurements were cal- 
ibrated, according to the time-domain calibration procedure described in subsection 
WT2\ 

4.3.3.1 Time-Domain Approach 

As mentioned in the previous section, the apparent dielectric permittivity of fuel, 
£app,f , was extrapolated from the amplitude of the first reflection in the TDR wave- 
form, caused by the air-to-fuel impedance change. As an example, Fig. 14.81 shows 
the TDR waveform (together with its derivative), when the probe was immersed in 
sample #1. As highlighted in the figure, the significant reflection points are asso- 
ciated to the level estimation derived from the distance difference between peaks 
of the derivative curve. Based on the £ ap p.f evaluation and measuring the apparent 
distance associated to the fuel, the fuel level (Hf d ) was determined from (14. 2M 1 ! 

The water level, H^, was therefore simply derived from the difference between 
the overall probe length and the estimated actual fuel level. Table l4. 3 I summarizes the 
results obtained solely from the TDR waveforms analysis described in the previous 
section flaLZl]. The deviation between reference values and the TDR-estimated values 
is also reported. 

Table 4.3 Summarized results, obtained through the TD-based approach, for the three con- 
sidered samples. The apparent permittivity of fuel (£ apPi f), the estimated level of fuel (H]r ), 
and the estimated level of tap water (H$ ) are reported. Deviations of estimated levels from 
reference levels are also reported 





app,t 


(mm) 


(mm) 


deviation 
(mm) 


Sample #1 


2.23 


281 


152 


±2 


Sample #2 


2.22 


219 


214 





Sample #3 


2.23 


331 


102 


±4 



The apex td indicates that results were obtained from time-domain measurements only. 
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Fig. 4.8 Reflection coefficient (and corresponding first derivative) acquired with the probe 
immersed in sample #1 



4.3.3.2 Frequency-Domain Approach 

For a preliminary validation of the minimization algorithm, the reference parameters 
for the three samples of liquids were used to generate the ideal scattering parame- 
ters: the modeled Su(f) would correspond to those that would be measured under 
ideal conditions, that is without any uncertainty due to imperfect calibration, im- 
perfect equipment, noise, etc. (see Table 14.41 ). These ideal scattering functions were 
considered as the input for the subsequent optimization routine, thus obtaining the 
extrapolated results reported in Table l4.4l Analysis of the data reported in Table l4.4l 



Table 4.4 Fuel level (Hf), water level (// w ), and Cole-Cole parameters extrapolated from 
'simulated' frequency-domain measurements and compared to reference values 







Sample #1 




Sample #2 




Sample #3 




Ref. 


Extr. 


Ref. 


Extr. 


Ref. 


Extr. 


£ s ,f 


2.35 


2.35 


2.35 


2.35 


2.35 


2.35 


£~,f 


2.02 


2.03 


2.02 


2.13 


2.02 


1.99 


/ri(GHz) 


13.0 


12.1 


13.0 


6.6 


13.0 


14.9 


]8f 


0.25 


0.25 


0.25 


0.20 


0.25 


0.25 


iff (mm) 


283.0 


283.0 


219.0 


219.0 


335.0 


334.9 


£s,w 


79.8 


79.8 


79.8 


79.8 


79.8 


79.7 


£oo 5 w 


5.2 


6.1 


5.2 


6.1 


5.2 


3.9 


/i-.w (GHz) 


17.3 


17.1 


17.3 


17.1 


17.3 


18.5 


/3w 


0.00 


0.00 


0.00 


0.00 


0.00 


0.02 


rJo.w (mS itT 


L ) 2.58 


2.58 


2.58 


2.59 


2.58 


2.44 


fl w (lMl) 


150.0 


150.0 


214.0 


214.0 


98.0 


98.1 
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shows that, due to the limited measurement bandwidth (from 10 MHz to 700 MHz), 
the evaluation of &. and of f T is somewhat less accurate. On the other hand, both 
the e s and, most importantly, the liquid levels are estimated with high accuracy. 
This demonstrates that the proposed setup is definitely suitable to estimate the static 
permittivity and level of fuel. 

When £„<, and f r are of interest, the probe design should be modified, and a coaxial 
probe with smaller radius should be adopted in order to achieve larger working 
bandwidths, without any modification in the optimization procedure (this case will 
be considered in Sect. 14.41 ). 

After the proposed optimization technique was tested on the idealized data, 
demonstrating its reliability and accuracy, it was applied to 'real' measurements 
performed on the same samples. In particular, the procedure was implemented con- 
sidering as input first the scattering parameter directly measured by the VNA, and 
successively, the one obtained through the FFT-based algorithm of TD data (i.e., 
the TD/FD combined approach) [5]. The obtained results are reported in Table 14.51 
and Table 14.61 for VNA and FFT data, respectively. Analysis of the two tables 
demonstrates that the optimization procedure can suitably estimate fuel and wa- 
ter levels with a maximum deviation within the uncertainty in the evaluation of the 
reference levels. In particular, it is important to point out that the TD/FD combined 
approach does not cause substantial accuracy decrease with respect to the VNA 
measurements. This confirms that a low-cost TDR system, associated to the 

Table 4.5 Fuel static permittivity, fuel level, and water level extrapolated for the three sam- 
ples using frequency-domain minimization routine, applied on VNA measurements. Devia- 
tions from reference levels are also reported 





c vna 


rjvna 

(mm) 


(mm) 


deviation 
(mm) 


Sample #1 


2.37 


283.5 


149.5 


±0.5 


Sample #2 


2.24 


219.5 


213.5 


±0.5 


Sample #3 


2.17 


335.1 


97.9 


±0.1 



Table 4.6 Fuel static permittivity (£*), fuel level (WjF'), and water level (H^) extrapolated 
for the three samples using frequency-domain minimization routine, applied on FFT of TDR 
measurements. Deviations from reference levels are also reported 





c fft 

t s.f 


(mm) 


(mm) 


deviation 
(mm) 


Sample #1 


2.21 


282.7 


150.3 


±0.3 


Sample #2 


2.18 


219.3 


213.7 


±0.3 


Sample #3 


2.16 


334.7 


98.3 


±0.3 
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Fig. 4.9 Reflection scattering parameter for sample #1 as obtained from VNA measurements 
and evaluated through the optimization procedure, a magnitude, b phase 



proposed procedure, is the key issue for the practical implementation in industrial- 
related applications. On a side note, the decreasing trend of static permittivity from 
sample #1 to sample #3 can be explained by the steady increase of room tempera- 
ture during the experimental session. However, temperature increase was so small to 
cause only minimal changes in static permittivity, not detectable in the TDR traces, 
which only allow an approximate average permittivity (e ap p) to be estimated. As a 
matter of fact, this confirms the suitability of the frequency-domain approach also 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
/(GHz) 
(a) 




/(GHz) 
(b) 

Fig. 4.10 Reflection scattering parameter for sample #1 as obtained from FFT of TDR mea- 
surements and evaluated through the optimization procedure, a magnitude, b phase 



for sensing small deviations from the standard reference status, even when they are 
caused by the temperature variations. 

As an example, Fig. 14.91 and Fig. 14. 101 compare the optimized Sn(/) for sam- 
ple #1, with measurements obtained from VNA and from the TD/FD combined 
approach, respectively. 

It is worth noting that, although the S\ i (/) directly obtained from VNA mea- 
surements, as expected, appears to be less noisy than the Si i (/) obtained from the 
TD/FD combined approach, this does not have any relevant impact on the accuracy 
in the evaluation of quantitative and dielectric parameters. 
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4.3.3.3 Validation of the Robustness of the Procedure 

It is important to point out that the adopted optimization procedure converges to a 
final solution that depends on the initial estimates and on the intermediate conver- 
gence point of the random part of the optimization procedure itself. To assess the 
robustness of the overall procedure, it was repeatedly applied to the same sample, 
slightly varying the initial guesses for the Cole-Cole parameters of water and fuel. 
The obtained results, in terms of extrapolated fuel levels, show dispersion within 
0. 1 mm, thus demonstrating the robustness of the procedure. To summarize the ob- 
tained results and to compare the different methods, Table 14.71 compares the fuel 
levels extrapolated from the time-domain and from the two frequency-domain ap- 
proaches, with the reference values. 



Table 4.7 Reference and extrapolated fuel levels, in time-domain and frequency-domain ap- 
proaches, for the three samples 



Fuel level (mm) 




reference 


time 
domain 


frequency 

domain 

(VNA) 


frequency 

domain 

(TD/FD combined approach) 


Sample #1 


283.0 


281 


283.5 


282.7 


Sample #2 


219.0 


219 


219.5 


219.3 


Sample #3 


335.0 


331 


335.1 


334.7 



The total time required for the optimization (encompassing the steps outlined) 
was of the order of 10 minutes. 

Indeed, the time required for optimization might be dramatically reduced by us- 
ing a customized software routine rather than resorting to the commercial circuit 
simulator. 

It is worth mentioning that, in a practical monitoring application, after the first 
run of the optimizer has been completed and initial Cole-Cole parameters and liquid 
levels have been estimated, the continuous monitoring simply requires the optimizer 
to track the small changes in dielectric parameters and level that occur in time. This 
means that the further optimizations will only require a few iterations of the simplex 
algorithm, which however should not take more than a few seconds^. 



As a matter of fact, different minimization methods, such as the globalized bounded 
Nelder-Mead (GBNM), can be alternatively used. Nevertheless, this would probably im- 
prove the efficiency of the first run, but would probably have no relevant impact in the sub- 
sequent continuous monitoring, when the optimization routine should just track smooth 
temporal variations in liquid levels and permittivities. 
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4.4 Dielectric Spectroscopy of Edible Liquids: The Vegetable 
Oils Case-Study 

As previously discussed, the in situ evaluation of the Cole-Cole parameters is not 
an easy task, especially when reference data are missing. This is the case, for exam- 
ple, with vegetable oils, for whose dielectric characteristics only scarce reference 



data are available B19ll23ll . In fact, the complexity of the intrinsic characteristics of 
vegetable oils, such as density and viscosity, makes the investigation of this kind 
of material rather challenging, nonetheless interesting. Additionally, since quality 
control of vegetable oils is becoming more stringent (especially for avoiding adul- 
teration), the estimation of these dielectric parameters may be used as a reliable 
indicator for certifying the product quality. 

It is worth pointing out that although highly sophisticated methods for the anal- 
ysis of edible oils are available, many of them are not easily applicable for routine 
or continuous monitoring (e.g.. gas chromatography [20] and liquid chromatogra- 
phy 19D), whereas others (e.g., Fourier transform infrared (FTIR) I125tl . Raman spec- 



troscopy 11411 . and nuclear magnetic resonance (NMR) [24]) are rather complex and 
highly expensive. On the other hand, other simpler methodologies that are com- 
monly adopted for routine analysis (such as those based on chemical titrations), 
do not provide much information. Additionally, these methods of analysis require 
skilled operators, and, most importantly, they cannot provide any information about 
the possible presence of adulterants. As a matter of fact, all of aforementioned meth- 
ods are affected by some limitations: in particular, they cannot be performed on the 
process line and they are very laborious. 

On such bases, not only does the presented case study aim at providing a prelim- 
inary assessment of the use of microwave dielectric spectroscopy for the characteri- 
zation of vegetable oils, but it also investigates the possibility of adopting a suitable 
measurement procedure for a rapid and reliable evaluation of dielectric parameters. 
This is done in view of possible practical applications, which may be useful, for 
example, in the on-line monitoring of the characteristics of oils throughout a pro- 
duction process. 

The approach proposed for the accurate evaluation of the Cole-Cole parameters 
is based on a procedure similar to the one described in the previous section, and it 
comprises three major steps: 

• First, the TDR waveforms are acquired and processed in order to obtain the cor- 
responding reflection scattering parameter, Sn (/); 

• Secondly, the measurement cell (i.e., probe and filling liquid) is accurately mod- 
eled as a TL, in which the dielectric characteristics of the considered liquid are 
parameterized through the Cole-Cole formula; 

• Finally, the Cole-Cole parameters of the LUT are evaluated by minimizing 
the deviations between 1) the measured reflection scattering parameter, Sn(/), 
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obtained from the TD/FD approach; and 2) the modeled reflection scattering 
parameter obtained from the implemented TL model of the probe plus LUT, 
5h,mod(/)- 

It is worth pointing out that, for this specific monitoring purpose, the TL model 
includes has to include not only the Cole-Cole parameters of the LUT as minimiza- 
tion variables, but also some lumped circuit elements that account for additional 
parasitic effects that cannot be compensated for through traditional SOL calibration 
procedures. 

This procedure was first tested and validated on well-referenced liquids, and suc- 
cessively it was applied for the evaluation of the Cole-Cole parameters of oils. The 
main steps of the procedure are similar to those in the diagram in Fig. |3.19| and will 
be described in detail in the following subsections. 

To define the performance characteristics of this combined approach (particu- 
larly, in terms of repeatability and uncertainty), a rigorous metrological validation 
was also carried out through measurements on well-referenced materials. 

TDR measurements on vegetable oils were carried out through the TDR80E04. 
However, in view of possible practical implementation, it is worth underlining that 
the procedure described herein also apply for portable TDR instrumentation. 



4.4.1 Design of the Measurement Cell and Transmission Line 
Modeling 

Although the SOL calibration is the first strategy for reducing the systematic er- 
ror contributions, a further accuracy enhancement must take into account the par- 
asitic effects associated with the measurement cell. For this purpose, an accurate 
TL modeling of the specific measurement cell represents a successful key for a fur- 
ther enhancement of the final measurement accuracy. The implemented TL model 
must include all the undesired effects (e.g., undesired impedance mismatches) that 
inevitably influence the measurements. 

A picture of the probe designed for this specific application is shown in Fig. l4.1l| 
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The probe was short-circuited at the distal end; the schematic configuration of 
the probe is shown in Fig. 14.121 



lo SMA 3.5 mm 

COMIKX'tOl 



short circuiting 
term itiation 



Fig. 4.12 Schematic config- 
uration of the used coaxial 
probe. 




4,4 nun 



The TL model of the measurement cell was implemented in the MWO simulator. 
The basic TL model includes a length of short-circuited coaxial probe, with the 
physical dimensions of the actual probe, filled with a material characterized through 
its relative permittivity and loss tangent (whose frequency-dependent variation is 
described by (l238l ». 

Preliminary TDR measurements on well-referenced materials showed that two 
major parasitic effects were apparent in the acquired waveforms: one is caused by 
the transition probe/SMA connector, whereas the other is due to the non-ideality 
of the short circuit at the end of the probe. Therefore, to take into account these 
parasitics, the TL model was appropriately adjusted; in particular, the introduction 
of some lumped circuit elements appropriately resembled the contributions noticed 
in the actual measurements. 

In order to obtain a TL model where the only variable parameters would be the 
Cole-Cole parameters of the LUT, the first step was the optimization of the model 
(i.e., the individuation of the values of the added lumped elements). For this pur- 
pose, preliminary measurements were carried out on three well-referenced mate- 
rials, namely 1,1,1-trichloroethane, ethyl acetate, and air: for these materials, the 
Cole-Cole parameters were taken from the related literature |3[] and were consid- 
ered as known values. Through an optimization procedure detailed in Sect. 14.3.11 
the squared deviations between the measured S\i(f) and the modeled Sh.mod(/) 
were minimized, thus estimating the values of the lumped elements that appropri- 
ately described the parasitics of the measurement cell. 

Minimization was performed over a frequency rangqj that goes up to 4 GHz. As 
an example, Fig. 14.131 shows the comparison between the measured Su(f) and the 



Preliminary minimization routines were performed over different frequency ranges, and 
the 0-4 GHz frequency range proved the optimal frequency range both in terms of mini- 
mum of the objective function of the minimization and in terms of appropriateness of the 
optimized circuit elements. 



74 



4 Quantitative and Qualitative Characterization of Liquid Materials 




Ethyl acetate 

— ■ — TL model 

— • — TD/FD combined approach 



1 



2 
/(GHz) 

(a) 



180 



90 



r o 

s 



-90 



\ 






H \ 


\ 






T 


\ 






1 


\ 


-1 1 — ' 


1 — 1 ' — 


— i r- 51 — 



-180 

12 3 4 

/(GHz) 
(b) 

Fig. 4.13 Comparison between the modeled scattering parameter (black squares) and the 
measured Si i (/) (grey dots), for ethyl acetate, a magnitude, b phase. The rinse between the 
TL-modeled S\\ (/) and the measured S\\ (/) is 0.006, over the 0-4 GHz frequency range [80 



modeled Sn,MOto(f) (after the minimization procedure was carried out) for ethyl 
acetate: it can be seen that the reflection scattering parameter obtained from the TL 
model closely agrees with the measured scattering parameter, both in magnitude 
and phase. The root mean square error (rmse) between S\\(f) and S\\ mod(/) was 
evaluated as a figure of merit, and its value was 0.006. 
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The final TL model and the values of the added lumped (parasitic) elements are 
shown in Fig. 14.141 In the implemented model, TL\ represents the SMA 3.5mm 
adapter; TLi represents the probe filled with the LUT; C\, L\ represent the non- 
ideal transition adapter/probe; and R\, R2, L2 represent the non-ideal short circuit at 
the end of the probe. C\ is proportional to the real part of the dielectric permittivity 
of the considered liquid, Real(e^ quid ). 



r> ^ > - 



L ''■■- TL 



Cj-Real^*) 



Reference plane for the \ 
calibrated S n (f) measurements 



j_ 



SMA 3.5 mm 
;er transition 




C, =0.155 pF 
£, = 0.524 nH 



R, = 0.150 n 
R 2 = 2.464 £3 
L, = 0.085 nH 



Non-ideal short 
circuit 



TL 2 andg*j. ^ depend on the specific liquid under test 
Fig. 4.14 Transmission line model (as implemented in MWO) of the used measurement cell 



4.4.1.1 Validation of the Method through Measurements on 
Well-Referenced Materials 

The appropriateness of the TL model in estimating the Cole-Cole parameters of 
the LUTs was verified through further measurements on 1,1,1 -trichloroethane and 
ethyl acetate. In this case, starting from the implemented TL model (now includ- 
ing all the evaluated parasitic components), the Cole-Cole parameters of the liquid 
were considered as the minimization variables. TDR measurements (along with the 
corresponding frequency domain processing and minimization) were repeated ten 
times for each considered material (each time, by emptying and refilling the probe). 
In this way, it was possible to assess the repeatability of measurements and to 
reduce random errors contributions. The four evaluated Cole-Cole parameters are 
summarized in Table l4.8l and Table |4.9l (static conductivity is negligible for both liq- 
uids): each set of parameters was evaluated as the average of the results of ten min- 
imizations over the ten measurements. In order to identify a figure of merit for the 
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metrological performance of the entire measurement procedure, the type A evalua- 
tion of the uncertainty (ma) was conducted for each Cole-Cole parameter, according 
to the following expression jig] : 



u A 



»(n-l)S ( 



(4.6) 



where n is the number of measurements, xi is the result of each evaluated parameter, 
and x is the average of the n repetitions. In the measurements reported herein, the 
number of repetitions n is 10. 

Table 4.8 Comparison between the reference and the evaluated Cole-Cole parameters (at 
20.0°C) for ethyl acetate. The experimental uncertainties for the evaluated values are also 
reported 



ethyl acetate 


parameter 


reference 
value 


averaged 
evaluated 
value 


experimental 
uncertainty 


£s 


6. 04 ±0.02 


6.00 


0.01 


6*> 


2.48 ±0.09 


2.71 


0.10 


/r(GHz) 


36.2±0.8 


29.6 


1.0 


]8 


0.06±0.01 


0.00 


0.00 



Table 4.9 Comparison between the reference and the evaluated Cole-Cole parameters (at 
20.0°C) for 1,1,1-trichloroethane. The experimental uncertainties for the evaluated values are 
also reported 



1,1,1- 
trichloroethane 


parameter 


reference 
value 


averaged 
evaluated 
value 


experimental 
uncertainty 


£s 


7.20±0.10 


7.25 


0.02 


&o 


2.08±0.02 


2.12 


0.15 


A (GHz) 


28.9±0.6 


26.7 


0.9 


]8 


0.030 ±0.002 


0.00 


0.00 



From Table l4.8l and Table |4.9l the calculated values of the static permittivity (e s ) 
closely agree with those reported in the literature; additionally, the corresponding 
experimental uncertainty (ma) is remarkably low. Similarly, the values obtained for 
Eoo and for f T are in good agreement with the literature values (which were obtained 
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through VNA-based measurements), despite being affected by a higher uncertainty. 
The higher uncertainty is most probably due to the fact that the relaxation frequen- 
cies of 1,1,1-trichloroethane and of ethyl acetate are much higher than the upper 
limit individuated for the optimization of the model (i.e., 4 GHz) and even higher 
than the frequency limit of the used TDR instrument (i.e., 15 GHz). In fact, addi- 



tional measurements on ethanol (whose f r is in the order of 1 GHz) 1 12], provided 
accurate results for all the four Cole-Cole parameters, thus validating the appropri- 
ateness of the proposed procedure. The four Cole-Cole parameters for ethanol and 
the related experimental uncertainties are reported in Table 14.101 



Table 4.10 Comparison between the reference and the evaluated Cole-Cole parameters (at 
20.8°C) for ethanol. The experimental uncertainties for the evaluated values are also reported 



ethanol 


parameter 


reference 
value 

CI 


averaged 
evaluated 
value 


experimental 
uncertainty 


£s 


25.04 ±0.04 


24.57 


0.02 


£*, 


4.481 ±0.027 


4.15 


0.05 


/i-(GHz) 


0.856 ±0.004 


0.92 


0.01 


]8 


0.00 


0.00 


0.00 



4.4.2 Experimental Results for Vegetable Oils 

When the proposed procedure was first applied to vegetable oils, the resulting val- 
ues of the dielectric parameters were affected by a relatively high uncertainty. Ad- 
ditionally, the obtained results were rather different from the data available in the 
literature. In particular, the values of the relaxation frequencies were systematically 
overestimated. 

The reason for this discrepancy could be mostly attributable to some intrinsic 
properties of the oils: these characteristics make the minimization over the four 
Cole-Cole parameter^] rather troublesome. In fact, oils are known to exhibit a re- 
laxation frequencjo lower than 1 GHz; additionally, their dielectric properties antic- 
ipate a low value of the dielectric increment (e s — £„<,). 

On such bases, to reduce the number of variables to optimize, the values of the 
static permittivity were assessed through an alternative method, and the minimization 
routine was performed over the remaining three Cole-Cole parameters (i.e., /3 , £*,, and 
f T ). The evaluation of the e s of the oils was performed through independent capacitive 
measurements carried out with an LCR meter. This strategy was deemed as the best 
solution to obtain accurate results over the whole considered frequency range. 



For vegetable oils, static electrical conductivity is practically zero. 
5 As reported in 1 19], vegetable oils exhibit a first dielectric relaxation at approximately 
1 MHz; however, this first relaxation was not considered in the analysis reported herein. 
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4.4.2.1 Evaluation of the Static Permittivity through LCR Measurements 

LCR measurements were performed using the same coaxial probe, but removing the 
short circuit at the distal end. The probe (filled with the LUT) was considered as a 
cylindrical capacitor whose capacitance (Cliquy) ideally depends on the liquid filling 
the space between the cylindrical plateqj, according to the following equation: 

C ^ = 2 «^WIa) (4 - 7) 

where £o is the permittivity of free space; e s lqul is the static permittivity of the 
LUT; h is the length of cylindrical plates (i.e., the length of the probe); a is the 
external diameter of the inner conductor; and b is the inner diameter of the external 
conductor. 

The measured capacitance (C meas ) can be regarded as the parallel equivalent of 
three separate contributions: 

Lmeas = Miq U id + Cpar t Cf r (4.o) 

where C par is the parasitic capacitance (mostly due to the SMA adapter used to 
connect the probe), and Cf r is the capacitance due to the fringing field at the probe 



end I13H - LCR-based capacitive measurements with the probe in air at increasing 
frequencies showed that, after a slight initial decrease, above 1 kHz the measured 
capacitance was constant, thus indicating that the effect of parasitic and fringing 
capacitances was stable; therefore, measurements at 1 kHz were used to estimate 
the value of C pa r + Cf r . This way, capacitance measurements on the materials under 
test could be corrected, thus obtaining Cii qu id and the e s of the liquid sample. 

Successively, the LCR-based approach for the evaluation of £ s was validated 
through preliminary measurements on reference materials (i.e., ethyl acetate, 1,1,1- 
trichloroethane, chloroform); measurements were repeated ten times for each con- 
sidered material, thus compensating for random errors contributions. Results for the 
evaluation of the values of e s are summarized in Table 14.111 and are compared to 
the literature values. The LCR-based method was used for the evaluation of the e s 
values of oils; also in this case, measurements were repeated ten times for each type 
of oil. 



4.4.2.2 Evaluation of the Remaining Cole-Cole Parameters for Vegetable 
Oils 

The values of the static permittivity were considered as known (fixed) values in the 
subsequent minimization routines carried out on FFT-transformed TDR data, em- 
ploying the previously described approach: as expected, since these minimizations 



LCR impedance measurements at low frequencies (i.e., at 1 kHz), confirmed the capacitor- 
like behavior of the probe; in fact, the phase angle of the impedance was -90°. 
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Table 4.11 Comparison between the reference and the LCR-evaluated values of e s , for the 
reference materials (at 20.0°C). Experimental uncertainty is also reported 



material 



reference 
value 



bibliographic 
reference 



averaged 
evaluated 
value 



experimental 
uncertainty 



chloroform 



4.82 ±0.02 



Id 



jui 



4.81 



0.02 



ethyl acetate 



6.04 ±0.02 



IHf 



6.04 



0.02 



1,1,1-trichloroethane 7.20±0.10 



7.42 



0.04 



were performed only on three Cole-Cole parameters, the accuracy of the final results 
was enhanced. 

Also in this case, TDR measurements were repeated ten times for each type of 
oil; for each measurement, the minimization routine was performed and the corre- 
sponding values of £«,, f T and /3 were evaluated. The minimization procedure over 
the Cole-Cole parameters was performed on the 0-1 GHz frequency range; in fact, 
this frequency range was found to adequately include the major dielectric response 
of oils. 

Ten different types of vegetable oils were considered: in particular, four olive 
oils with different acidity values (ranging from 0.3% to 4.0%) and six other kinds 
of vegetable oils (i.e., peanut, sunflower, corn, castor, various seed, and soybean). 

Table 14.121 summarizes the Cole-Cole parameters for the ten oils (each set of 
parameters was calculated as the average of ten repetitions), and the correspond- 
ing standard deviation (a). The obtained results for the Cole-Cole parameters 
of the considered vegetable oils are in good agreement with the few data available 
in the literature Ill0lll9ll23ll . It is worth noting that, for the results obtained for oils, 
the standard deviation was reported (rather than the experimental uncertainty), so as 
to provide an indicator of the metrological performance of the method, when used 
without performing a series of repeated measurements. 

As expected, from Table l4~T2l it is apparent that /3, £«,, and e s change slightly 
among the different oils. Conversely, / r seems to change significantly between one 
type of oil and another. To better evidence this trend, for each oil, Fig. 14.151 shows 
the expanded uncertainty range ±2.26 a; where a is the standard deviation cor- 
responding to the Gaussian probability distribution (which, in turn, was verified 
through the ^-square test). The expanded uncertainty was evaluated considering the 
r-distribution of Student, and 2.26 is the f-score that corresponds to a confidence 
level of 95% and to nine degrees of freedom llql . On the basis of the obtained 
results, the proposed technique can be useful for quality monitoring of oils: for ex- 
ample, by measuring the relaxation frequency, it might be possible to discriminate 
different kinds of oil. Additionally, using the averaged relaxation frequency as a 
reference value (i.e., the expected value that the considered oil should exhibit), any 
deviation may be associated with the possible presence of adulterants. 



so 



4 Quantitative and Qualitative Characterization of Liquid Materials 



Table 4.12 Averaged Cole-Cole parameters for the ten types of vegetable oil (at 20.0°C). 
Standard deviation values are also reported 



type of oil 


fs 


Oe, 


£„ 


o&. 


ft 

(MHz) 


(MHz) 


P 


<Jp 


olive (ac.=0.3%) 


3.08 


0.01 


2.39 


0.01 


315 


11 


0.33 


0.01 


olive (ac.=1.2%) 


3.14 


0.01 


2.38 


0.01 


288 


9 


0.36 


0.02 


olive (ac.=1.6%) 


3.19 


0.03 


2.36 


0.01 


259 


4 


0.40 


0.01 


olive (ac.=4.0%) 


3.19 


0.03 


2.34 


0.02 


249 


5 


0.42 


0.02 


peanut 


3.05 


0.01 


2.40 


0.01 


334 


5 


0.28 


0.01 


sunflower 


3.12 


0.04 


2.40 


0.01 


292 


5 


0.31 


0.01 


corn 


3.11 


0.02 


2.41 


0.01 


309 


6 


0.29 


0.02 


castor 


4.69 


0.01 


2.56 


0.01 


122 


2 


0.42 


0.01 


various seeds 


3.10 


0.02 


2.43 


0.01 


371 


5 


0.27 


0.01 


soybean 


3.09 


0.01 


2.41 


0.02 


390 


7 


0.30 


0.02 
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Fig. 4.15 Summarized comparative results of the averaged values of f r for the considered 
oils. The corresponding expanded uncertainty bars (with a confidence level of 95%) are also 
reported 



On a side note, it is important to point out a few aspects related to the used probe. 
Generally, a short probe provides several benefits: a lower attenuation of the elec- 
tromagnetic signal is achieved, small volumes of liquid are required, and, finally, 
a higher mechanical stability between inner and outer conductors can be ensured. 
On the other hand, short probes do not provide enough electrical length for the 
electromagnetic signal in the LUT at low frequencies: this may represent an issue 
when determining the value of the static permittivity, e s . In this regard, the choice 
of the optimal length of the probe should start from preliminary considerations of 
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the intended purposes. As a matter of fact, the choice of the most appropriate probe 
length should guarantee an optimal trade-off of contrasting effects, such as an ad- 
equate performance at low frequency and the need to limit the quantity of oil re- 
quired for performing the tests. This alternative solution may as well be regarded as 
a practical alternative for avoiding the use of the LCR meter for evaluating the static 
permittivity. 

4.4.2.3 Results on Mixtures of Oils 

After the validation of the proposed procedure for quality-control of pure olive 
oils, the successive phase was to extend the methodology to more complex systems 
(such as mixtures of oils), so as to test the method also for anti-adulteration control 
purposes. 

On such bases, the proposed procedure was applied for the evaluation of the 
Cole-Cole parameters of some mixtures of oils (Ml, ..., M4). In particular, to re- 
produce some typical adulteration cases, the mixtures were prepared by mixing a 
new olive oil with sunflower oil, in different volumetric percentages. Measurements 
were performed one single time for each mixture; this was done to mimic a typical 
practical case in which the personnel in charge of control have to rely on single-shot 
measurements. Summarized results are reported in Table l4.13l Also in this case, it is 
apparent that the parameter that is more sensitive to the change in the composition 
of the mixture is the relaxation frequency. More specifically, it can be seen that start- 
ing from pure olive oil, the relaxation frequency value increases proportionally to 
the percentage amount of sunflower oil. On the other hand, the remaining Cole-Cole 
parameters do not exhibit significant changes. 

On a final note, it is worth noting that results show that the difference in / r for 
the considered oil mixtures is of the order of a few MHz (see Table 14.131 ). This 
figure is comparable to the achievable uncertainty for measurements on pure oils, 
when single-shot measurements are performed (see Table l4. 12t . Overall, it appears 
that a further minimization of the uncertainty on / r evaluation would be necessary 
to make the method more effective. This could be achieved, for example, exploit- 
ing enhanced probe geometries, as currently under development. Moreover, suitable 

Table 4.13 Relaxation frequencies for the olive-sunflower oil mixtures 



Mixture 


Sunflower oil percentage (%) 


ft 

(MHz) 


olive oil 





245 


Ml 


10 


250 


M2 


20 


257 


M3 


30 


259 


M4 


40 


263 


Sunflower oil 


100 


295 
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data analysis techniques (such as Principal Component Analysis, PCA) might be 
employed to distinguish adulterated oils from pure olive oils on the basis of the 
whole combined set of Cole-Cole parameters. 

4.4.2.4 Considerations on Possible Practical Implementation 

Results show that different oils exhibit different dielectric characteristics, particu- 
larly in terms of relaxation frequency. The uncertainty value associated with the 
evaluation of the relaxation frequency provides a self-consistent parameter for test- 
ing the performance and the suitability of dielectric spectroscopy for quality control 
of vegetable oils. 

For this purpose, results may be collected in a database (together with the as- 
sociated confidence intervals) and used as reference values for successive 'offline' 
measurements on unknown oils. 

Most importantly, results on mixture of oils demonstrate that dielectric spec- 
troscopy is sensitive to the presence of adulterants in oils. In practical applications, 
provided that some reference dielectric data are given for specific unadulterated oil, 
the possible presence of adulterants may be inferred, for example, as a deviation 
between reference and measured dielectric data. The major advantage of employing 
this method would be that the personnel in charge of control may perform the check 
directly in situ, without overwhelming the few available laboratories with unnec- 
essary requests. In other words, if the dielectric spectroscopy-based analysis gives 
a 'positive' result (i.e., adulterant might be present), then further specific analyses 
may be required from dedicated laboratories; on the other hand, if the result is neg- 
ative, then there is no need to carry out additional expensive and time-consuming 
analyses. 
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Chapter 5 

Qualitative Characterization of Granular 

Materials and Moisture Measurements 



'There are two possible outcomes: if the result confirms the 
hypothesis, then you 've made a measurement. If the result is 
contrary to the hypothesis, then you 've made a discovery. ' 
Enrico Fermi 



Abstract. This chapter focuses on broadband microwave refiectometry (BMR) ap- 
plications for monitoring water content (0) and static electrical conductivity (Oq) of 
granular materials, with particular emphasis on applications for soil measurements. 
First, a TDR-based method for inferring from measurements of the apparent di- 
electric permittivity is presented. This approach, which relies on the individuation 
of so-called calibration curves, is discussed in detail and the related metrological as- 
sessment is provided. Successively, a more accurate method that takes into account 
the frequency-dependence of the dielectric permittivity of the moistened granular 
material (considering the permittivity of each single constituent) is presented. This 
application is also used as test-case for validating an innovative calibration proce- 
dure that is particualrly useful when the traditional short-open-load (SOL) calibra- 
tion cannot be performed. Furthermore, the adoption of antennas in place of the 
traditional probes is discussed (thus assessing the possibility of guaranteeing a non- 
invasive approach). Finally, two innovative strategies for enhancing and simplifying 
TDR-based measurements of static electrical conductivity Co (typically used in soil 
science) are presented. 



5.1 Introduction 

Measuring water content of granular materials represents one of the major practi- 
cal applications of broadband microwave refiectometry (BMR), especially for soil 
moisture measurements. 

Indeed, refiectometric techniques have become increasingly attractive solutions 
for soil moisture sensing, with applications that focus on environmental monitoring 
purposes Il60ll and agricultural water management Il8ll5lll . In fact, mostly thanks to 



the versatility of the technique, a wide variety of monitoring tasks can be effectively 
performed through BMR I14I1 . 
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Although other moisture sensing techniques could arguably provide as much 
measurement accuracy, their use is often hindered by practical limitation^]. Con- 
versely, BMR holds considerable potential for practical implementation; in particu- 
lar, time domain reflectometry (TDR)-based systems have become competitive with 
respect to other techniques, such as capacitive methods and ultrasonic sensors. 

TDR-based systems can be customized according to the operating conditions, 
thus providing high adaptability to the specific application. Such a feature, com- 
bined with low-cost instrumentation and with an adequate measurement accuracy, 
has contributed to the spread of TDR in soil science, not only for moisture content 



evaluation but also for determining the moisture content profile I23l 12911 . On such 



bases, relevant research effort has been dedicated to developing innovative TDR 



probes, enhanced measurement procedures, and sound theoretical modeling Q8H . 

The typical TDR-based instrumental setup for in situ moisture sensing includes a 
portable reflectometer and a three-rod probe that can be easily inserted in the soil (or 
granular material in general) under investigation, thus acting as measurement cell. 
As mentioned in the previous chapters, despite having some shortcomings with re- 
gards to calibration, this probe configuration is widely used for TDR-based moisture 
content measurements. 

The present chapter is structured as follows. First, the evaluation of 9 through 
TDR-measurements of the apparent dielectric constant, £ app , is considered. In par- 
ticular, the adoption of the so-called calibration curves, which relate e app to 9, is 
discussed: the metrological analysis of this method demonstrates that it can ensure 
accurate results, provided that for each different material a specific calibration curve 
is individuated. 

Secondly, the adoption of an enhanced approach that also takes into account 
the frequency-dependence of the dielectric permittivity of the investigated material 
is fully described. Similarly to the previously reported applications, the proposed 
method relies on the combination of TDR, frequency domain reflectometry (FDR), 
and a transmission line (TL) modeling of the measurement setup. Most importantly, 
this application was also used as a case-study for testing an innovative FD-based 
calibration technique that relies on three short-circuits applied at different probe 
sections. This calibration procedure is particularly important when multi-rod probes 
are used; in fact, in this case, traditional short-open-load (SOL) calibration is often 
practically unfeasible, due to the impossibility of connecting standard loads to the 
electrodes. 

Successively, a noninvasive approach for moisture monitoring is also assessed: 
this method focuses on the possibility of adopting antennas as sensing elements. 
More specifically, the antenna is placed on the top surface of the material under test 
(MUT), and the shift in the resonant frequency of the antenna is associated to the 
corresponding moisture content of the sample. 



1 As an example, the gamma ray attenuation method is accurate but highly costly, and the 
presence of a radioactive source requires competent handling thus making the method 
unappealing for routine applications 14511 . 
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Finally, the last section of this chapter describes an enhanced TDR-based method 
for measuring static electrical conductivity (Ob), especially useful for granular ma- 
terials, such as soil. In particular, two innovative approaches for simplifying and 
'speeding-up' the time-consuming preliminary procedures that have to be carried 
out for TDR-based estimation of Co are presented. The first method relies on the 
combination of the TDR measurements with a transmission-line modeling (TLM) 
of the measurement cell; the second method simply relies on a couple of indepen- 
dent capacitance measurements (ICM) performed through an LCR meter. 

5.2 Dielectric Models for the Estimation of Water Content 

Basically, in BMR, the value of the volumetric water content, 6 , is inferred from per- 
mittivity measurement^. In fact, the relative permittivity of water (approximately 
80) is considerably higher than the typical relative permittivity of many granular ma- 
terials of interest (e.g., soil); therefore, the presence of water changes considerably 
the overall dielectric permittivity of the MUT 

Different methods can be adopted for deriving from TDR measurements. The 
simplest approach is based on the evaluation of 6 from measurement of the apparent 
dielectric constant (£ app ) of the considered material llul . 

A simple approach for TDR-based soil moisture monitoring relies on the use of 
empirical material-specific calibration curves that provide the functional relation- 
ship 6 — £ ap p llul . These curves must be determined through preliminary measure- 
ments on samples moistened at pre-fixed known moisture levels (0 re f): a known 
volume of the considered material is moistened with fixed amounts of water and, 
for each level of moisture, the corresponding e app is retrieved from measurement 
of the apparent length. The points (0 re fj£a PP ) are fitted through a polynomial, and 
the obtained curve is used as reference for the successive measurements on samples 
of the same material, with unknown moisture levels lid] . As aforementioned, this 
approach is simple and adequately accurate. 

Another widespread approach relies on the adoption of empirical models, which 
fit empirical relationship to experimental data. These models have some shortcom- 
ings (such as the limited possibility of extrapolation outside the moisture range of 
the original set of experiments 1140); nevertheless, they are routinely used especially 
for soil-like materials. 

In particular, for soils, one of the most widely used relationship between e app and 



9 is the Topp's equation H46l 16311 : 



E app = 3.03 + 9.3 0+ 146 2 - 76.7 3 . (5.1) 



Volumetric water content is expressed as 8 = V W at/Vtot, where K W at is the volume of water, 
and V tot is the total volume of the mixture. The corresponding gravimetric water content 
can be calculated as (pwat/Pb) 9, where p wat is the density of water, and p\, is the bulk 
density of the moist granular material. 
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Topp's equation^ is sufficiently general (in fact, it app ears to be independent of 
soil texture and structure, salinity and temperature) [52]. However, although ( 15.1b 
has proved suitable for a wide range of homogeneous soils, it is inappropriate for 
clayey soils J19I1 and for soils with high content of organic materials 15611 . In these 
latter cases, a third degree polynomial may still be used, but the coefficients should 



be extrpolated through a fitting procedure 16211 . 



Malicki et al. suggested that the influence of the soil matrix was remarkable and, 



for this reason, they suggested including either dry bulk density or porosity B42I1 . 
However, the empirical functions used in current research for TDR data conversion 
still remain anything but universal models 11411 . 

A more accurate 6 — £ app relationship can be obtained through volumetric mixing 
models, which take into account the volume fraction and the dielectric constant of 
each single constituent that is present in the MUT. The dielectric mixing models are, 
therefore, particularly useful for the characterization of heterogeneous materials, 
such as soils with low dry bulk density, a large amount of bound water, or a relatively 
large permittivity of the solid phase. 

In this approach, soils are usually considered as three-phased systems (made up 
of soil particles, water and air) or as four-phased systems (considering that water 
close to the soil particles behaves differently) 14311 . 

The typical expression of a three-phased mixing model is as follows 

ei = ee^ + (l-<j))e^ + (<j)-9)ei (5.2) 

where the subscripts m, s, a, and w stand for measured, dry soil, air, and water, 
respectively; is the porosity; /3 is a parameter that depends on the spatial structure 
of the mixture and its orientation with the respect to the outer electric field [49]. 

The four-phased systems take into account that the water that is close to the 
soil particles has lower mobility than free water; hence, the thin layer of water that 
surrounds soil particles has a lower dielectric permittivity than free water 12011 . The 
expression of a typical four-phase mixing model is the following |68]: 

£ = (1 - <t>) ef + (<j> - 8) £ a « + (0 - ft,) e t l + 0b W C (5.3) 

where the subscripts m, s, a, fw, and bw stand for measured, dry soil, air, free water 
and bound water, respectively; a is a parameter that takes into account the geometry 
of the media particles with respect to the applied electromagnetic signal; and 0b w is 
the fraction of bound water, which can be obtained from 

9 bvj = n m 8p b S e (5.4) 



Practically, for estimating water content from TDR measurements of the apparent dielec- 
tric constant, the reciprocal of d5.lt is used: 

6=4.3x l(T 6 £ a 3 pp -5.5x l(T 4 £ a 2 pp + 2.92x l(T 2 £app -5.3 x 1(T 2 . 
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where n m is the number of monolayers of water molecules of thickness 8 bounded 
to the soil particles; S e is the soil specific area; and p\, = p s (l — (j>) is the soil bulk 
density, p s being the soil dry density 11911 . For complex media, such as soil, a « 0.5; 
nevertheless, in general, a is a fitting parameter J43ll . 

A further enhancement of these models can be achieved by taking into account 



the frequency-dependence of the permittivity for each of the soil constituents [68]. 
In fact, when dealing with complex materials or when a high percentage of water is 
present in the medium, the dielectric mixing model should also consider the depen- 
dence of the dielectric characteristics on frequency and the effect of Ob 1271 . Taking 
into account the frequency dependence of permittivity is usually accomplished with 
the so-called inverse modeling, which may be developed either directly in time do- 
main J32I1 or in frequency domain I9L I33I1 . 



Dielectric mixing models have proved to be suitable in many practical applica- 
tions and their perspectives for further use seem to be very promising, not only for 
soil-like materials, but for a variety of materials (e.g., materials for construction 



1141). 
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In this section, the estimation of moisture content from TDR-measurement of the 
apparent dielectric constant is addressed. In particular, the adoption of £ app — 6 cal- 
ibration curves is described in detail, and the strategies for deriving accurate and 
reliable curves are discussed. 

As aforementioned, for a given material, the e app — calibration curve can be 
assessed as follows. Samples of the considered material are moistened at prefixed 
values of moisture (<9 re f), and the corresponding apparent dielectric constant (e app ) 
is measured through the TDR method. The (e app , re f) points are fitted through, thus 
obtaining a regression curve. 

In this way, in successive moisture content measurements on the same type of 
material, it is enough to measure e app , and the corresponding moisture level is simply 
retrieved from the calibration curve, with a certain confidence level. However, it 
must be pointed out that calibration curves are specific for each kind of material; 
hence, different curves must be derived for different materials (also, for example, 
for different types of soil). 

As a matter of fact, different procedures are used for the extrapolation of empir- 
ical calibration curves, which, as a result, are not readily comparable. In fact, the 
individuation of an optimal calibration curve for each of the materials under inves- 
tigation is a key point to improve the accuracy of results. Despite the deep interest 
surrounding TDR, neither an assessment of standard measurement procedures, nor 
the corresponding metrological characterization has been thoroughly investigated. 

To fill this gap, the strategies for obtaining accurate calibration curves were 
specifically investigated. Furthermore, the metrological performance of the method 
was assessed through a statistical analysis on the results obtained through the two 
TDR instruments considered in this book (namely, the HL1500 and the TDR80E04). 
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Measurements were carried out on the same samples through each instrument: 
starting from measurements on well-referenced siliceous sand samples, the analysis 
was extended to agrofood granular materials (i.e., bran and corn flour). The com- 
parison of results and related uncertainty limitation led to a consistent assessment 
of the method from a metrological point of view. 

As already mentioned on several occasions, a metrological characterization must 
involve the consideration of both systematic and random errors. To minimize and 
to compensate for systematic error effects, a preliminary calibration of the experi- 
mental setup is mandatory (as detailed in Sect. l3.6.Tb . Nevertheless, further residual 
uncertainty contributions are mainly related to random effects. Therefore, to eval- 
uate the final uncertainty, a statistical analysis on repeated TDR-measurements of 
dielectric constant was approached. To carry out the analysis, the uncertainty prop- 
agation theory and the evaluation of the expanded uncertainty associated to the t- 
Student distribution were considered for the estimation of the final measurement 
uncertainty. 

Once the uncertainty, associated to the e app evaluation, was quantified for each 
measurement point, results deriving from two sets of measurements (repeated with 
the same instrument in different times) were unified in a single 6 — £ app fitting curve. 
Subsequently, a non-linear regression method was considered in order to obtain a 
specific e app — 6 regression (calibration) curve and the related uncertainty range 
corresponding to the confidence level of 95%. In this way, a robust method for the 
evaluation of the uncertainty limitation was provided, along with some practical 
hints for minimizing error contributions. 



5.3.1 Details on the Experimental Procedure 

The three materials considered for the comparative analysis were sand, corn flour, 
and bran; their mass densities were 1.53, 0.66 and 0.67 gem -3 , respectively. 

Each material was placed inside a 'large box' (approximately lmxlmx0.3m) 
and was moistened at pre-established levels, according to the volumetric method 
reported in 1471 . To achieve the pre-established moisture levels, an appropriate vol- 
ume of water was weighed through an electronic balance (with an uncertainty of 
0. 1 g), and added to the MUT Material and water were then mixed, so as to ensure 
a homogeneous moisture level: this was verified through preliminary TDR mea- 
surements of the dielectric constant in different points of the material contained in 
the 'large box'. Successively, a sample was removed from the box, placed inside a 
cylindrical sample holder (with a diameter of 7 cm, graduated up to 2000 ml) and 
characterized. Preliminary measurements confirmed that boundary effects due to 
the used holder were negligible. As a matter of fact, this is one of the advantages of 
TDR method over other well-known moisture content measurement methods, such 
as capacitive method or ultrasonic sensor. In fact, in the last two cases, the involved 
low-frequency range requires the use of large boxes intended to prevent interfer- 
ences caused by boundary effects. 



5.3 Evaluation of Moisture Content Directly from TDR 91 

TDR measurements were performed through a 15 cm-long commercial three-rod 
probe (CS630, Campbell Scientific). 

Before proceeding with the detailed discussion on the proposed experimental and 
metrological approaches, let us briefly summarize the key-points of the procedure. 
After a specific probe calibration (as reported in subsection l3.6.1l l. for each of the 
pre-established moisture levels (0 re f) at which the samples were moistened, the cor- 
responding e app value was evaluated. This way, the functional relationships between 
the apparent dielectric constant and the corresponding moisture level were identified 
through an appropriate fitting of the data. The overall uncertainty of e app was evalu- 
ated so as to draw a preliminary comparison between the different performances of 
the used TDR instruments and to understand to which extent results depend either 
on the instrument or on other experimental conditions. 

Successively, the fitted data were used to extrapolate specific calibration curves, 
for each considered material and for each instrument. Simultaneously, an evaluation 
of the overall uncertainty in the estimation of 6 from the calibration curves was 
done through the non-linear regression method. This was done to provide a compre- 
hensive panorama of the measurement system as a whole and to validate a general 
method which, in turn, can be applied also for other types of materials. 

As reported in |3] and according to the calibration procedure described in Sect. 
13.6.11 a ten-measurement set was preliminarily carried out with the three-rod probe- 
footnotemark immersed first in air and then in demineralized water. 



5.3.2 Uncertainty Evaluation for Apparent Dielectric Permittivity 
Measurement and Individuation of the Calibration Curves 

The evaluation of the apparent dielectric constant from TDR measurements relies 
on ( 13.5b . which is reported here for convenience: 

Lapp ^ (5.5) 

where L app is the apparent probe length in the MUT and L is the length of the probe 
in air. To take into account the uncertainty sources, a ten-measurement set of the 
apparent length L app was acquired for each sample. For each set of measurements, 
both the mean apparent length L app and the standard deviation o> were calcu- 
latecQ, then the expanded uncertainty TJi~. was evaluated through the f-Student dis- 
tribution, according to the relation: 

U LiaB = 2.26a, (5.6) 



A coaxial geometry for the probe would guarantee a reliable matched transmission line. 
Nevertheless, three-rod probes are definitely more practically usable when investigating 
granular materials, thanks to the easier insertion in the MUT. 
4 <Ti. t is the standard deviation corresponding to the Gaussian distribution, which was ver- 
ified through the ^-square test. 
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where 2.26 is the /-score associated with a confidence level of 95% and nine degrees 
of freedom. Once L app and E/z™, were both known for each of the ten-measurement 
sets, referring to Ul as the expanded uncertainty of the probe electrical length in 
air, the measurement uncertainty was evaluated using the uncertainty propagation 
theory B3411 . according to the following equation: 



U. 




\ 



-app(0 



Ui, 



[2 ^app ( ,. 



2^U L 



(5.7) 



In the above equation, U Ej is the propagated uncertainty associated to the j'-th mea- 
surement of the e ap p value. Finally, data with corresponding uncertainty bars were 
fitted in a 6 — e app curve. 

To determine the £ app — 6 relationship, the experimental data were fitted, thus 
obtaining the optimal calibration curve. As also reported in many related papers, 
third order polynomial curves prove to be a good fitting of experimental data: 



0=B O +Bie app +B2e : 



-app 



-B 3 e. 



app 



(5.8) 



where Bo, B\, B2 and B3 are the regression coefficients. This way, ( 15.81 ) represents 
the calibration curve: in correspondence of each measured dielectric constant value 
£ app , the curve provides the corresponding 6. 

To characterize the extrapolated value of 9 from a metrological point of view, 
the associated uncertainty was evaluated through the non-linear regression theory 
1 3411 . Hence, the variance analysis for the single values expected from the previous 
equation was conducted according to 



var [9] 



dd_ 

dB t 



var[Bi}+2^ 



do de 

~dBj~dB~j 



cov[Bi,Bj] (5.9) 



where var [6] is the variance of the moisture level; a 2 is the variance between ex- 
perimental and fitted data; n is the number of experimental points; i, j = 0,1 ,2,3 ; and 
cov [Bj,Bj] is the covariance between Bj and Bj parameters. 

For a confidence level of 95%, the equations associated to the lower and upper 
confidence limits for the regression curve are given by the following equations: 
alpha 



^low = {Bo + B\ e app + B2£ app + S3£ app ) — t n 



-4,1 



a v var [01 

T 



(5.10) 



5.3 Evaluation of Moisture Content Directly from TDR 93 



£u P = {B + Bie aw + B 2 e^ p + B i ej pp ) +t n _ 4A _ayJvar[d] (5.11) 

where a is the significance level (5%) [34]. 



5.3.3 Measurement Results and Instrumental Performance 
Comparison 

Preliminarily, to validate the practical implementation of the proposed methodology, 
repeatability tests were carried out on the considered material samples (sand, bran, 
and corn flour), with both the HL1500 and the TDR80E04. For each instrument, 
two series of measurements were performed, in different times, on each MUT: as 
expected, both for the HL1500 and for the TDR80E04, a good agreement between 
the first and the second series of measurements was observed. 

The repeated measurements were grouped to extrapolate specific fitting curves 
for each instrument. The resulting fitting curves were unified in a single 'global' 
curve for each material and for each instrument (Fig. 15.11) . and the corresponding 
uncertainty values were calculated for each measured e app . Particularly, to verify the 
adequacy of the normal function distribution in describing the uncertainty distribu- 
tion, the r-Student distribution was considered for the whole series of measurement 
data: the hypothesis was confirmed assuming a confidence level of 95% [34]. This 
way, in Fig. 15.11 to each measured £ app value the corresponding absolute uncertainty 
bar (U e ) is associated, exhibiting a maximum percentage relative standard uncer- 
tainty lower than 10% for the HL1500, and lower than 5% for the TDR80E04. As 
expected, the uncertainty affecting e app measured through the HL1500 is higher than 
the uncertainty affecting the e app obtained through the TDR80E04. The metrological 
characterization of the proposed method involves the analysis of regression curves 
on each material sample and for both HL1500 and TDR80E04. The associated con- 
fidence intervals are reported in Fig. s !5. 2145. 41 in terms of lower and upper limits of 
global calibration curves (Li ow and L uv ), as calculated from ( 15.101 ) and ( 15. 111 ). 

The graphs show the good agreement between measurements conducted through 
the HL1500 and through the TDR80E04; the corresponding value of the confidence 
band allows deducing the final uncertainty in moisture level evaluation. Results re- 
lated with regression curves for lower and upper limits of confidence bands are 
summarized in Table 15.11 It is worth pointing out that even though in some cases 

Table 5.1 Absolute uncertainty range of regression curve for each considered material, asso- 
ciated to a confidence level of 95%, calculated for the two TDR instruments 
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Fig. 5.1 Comparison between global fitting curves obtained using HL1500 and TDR80E04. 
a for sand, b bran, c corn flour. Absolute uncertainty bars are also reported 
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(a) 




(b) 



Fig. 5.2 Global regression curve for sand, a obtained using HL1500. b obtained using 
TDR80E04. Related confidence limits, associated to the 95% of probability occurrence, are 
also reported 

it appears that the distribution of the HL1500 measurement points is narrower than 
the TDR80E04 measurement points distribution around their respective extrapolat- 
ing polynomial, this is merely due to the fact that upper and lower limits mainly 
indicate how well the regression curve fits the experimental points, and should not 
be attributed to a better performance of the HL1500. 



5.3.3.1 Consideration on Practical Implementation 

As expected, TDR80E04 exhibits a better performance in terms of accuracy than its 
less-expensive portable counterpart, the HL1500. Nevertheless, the adoption of the 
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Fig. 5.3 Global regression curve for bran, a obtained using HL1500. b obtained using the 
TDR80E04. Related confidence limits, associated to the 95% of probability occurrence, are 
also reported 



low-cost portable HL1500, in conjunction with the proposed measurement proce- 
dure, guarantees adequate performance, thus allowing a reliable monitoring of the 
water content. These results are achieved thanks to the strategies described in Sect. 
13.6.11 such as a preliminary calibration procedure making the approach extremely 
robust, as well as techniques of waveform interpretation and compensation of resid- 
ual error contributions. 

As already demonstrated in the cases of liquid materials, this is particularly im- 
portant in view of possible industrial implementation, which not only does require 
accuracy and reliability, but also low cost and portability. 
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Fig. 5.4 Global regression curve for corn flour, a obtained using the HL1500. b obtained 
using the TDR80E04. Related confidence limits, associated to the 95% of probability occur- 
rence, are also reported 



5.4 Moisture Content Measurements through TD/FD 
Combined Approach and TL Modeling 

5.4.1 Moisture Content Measurements through Dielectric Mixing 
Model 



As mentioned in Sect. 15.21 a further enhancement of the accuracy in the estimation 
of moisture content can be achieved by considering the dielectric permittivity of 
each of the soil constituents (i.e., air, water, and soil particles) and their dependence 
on frequency. 
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In the proposed method, the three-rod probe (i.e., the 15 cm-long Campbell 
Scientific CS630 130, connected to the HL1500) inserted in the moistened soil 
(considered as a dispersive material) is modeled as a transmission-line. In partic- 
ular, considering a TEM mode propagating along the open-terminated probe, the 
frequency-dependent reflection coefficient Si 1 .mod (/I at the beginning of the rods 
(section A in Fig. 15.5b ) can be simply written as 15011 

Sn,MOD(/)=^ 2ir(/) ' (5.12) 

where / is the rod length, i 2 = — 1, and k*(f) = 27r/(/io£o£*c(/)) is the complex 
propagation constant, associated to jAq and £o, and to the complex relative permit- 
tivity of the soil £* c (/) . 

The key point for the implementation of Eq. (15.12b is the choice of an appropriate 
mixing model for £ s * c (/), thus accurately representing the dielectric characteristics 
of the soil-medium. 

As well known, the frequency-dependent dielectric permittivity of a homoge- 
neous material can be suitably described by the Debye equation 12611 . Conversely, 
the heterogeneity of soil dictates the use of more complex models that could pro- 
vide a more truthful representation of the dielectric behavior. Various models have 
been used to describe mixtures; nevertheless, no complete model is available that 
can describe the dielectric properties of soil B31I1 . A detailed review of the advances 



in soil mixing models can be found in 155I1 . 



The experiments reported in this section focused on sandy soils, for which the 
effect of bound water is negligible; therefore, a three-phased dielectric mixing model 
is enough to adequately describe the dielectric properties 16811 . More specifically, for 
sand moistened at a volumetric water content, the corresponding complex relative 



permittivity £ s * c (/) can be suitably described through the following equation B68I1 : 



< c (/) = [(l-<?)e S p + 0(4c(/)) a + (0-C] 1/a -i^ (5-13) 

where (j) is the porosity of sand; £ sp is the relative permittivity of the sand particles 
(assumed constant with frequency); £$, c (f) is the frequency-dependent complex rel- 
ative permittivity of water; £ a is the relative permittivity of air (assumed equal to 1); 
Ob is the static electrical conductivity of the moistened sand; and a is a fitting pa- 
rameter that takes into account the geometry of the sand particles. 



5.4.2 Triple-Short Calibration Procedure 

The adoption of dielectric mixing models for moisture content estimation has sev- 
eral major advantages; nevertheless, its measurement accuracy is strongly related to 
the adoption of a frequency-domain calibration (FDC) procedure, which typically 
involves the use of coaxial electrical standards, such as short-open-load (SOL). Gen- 
erally, for performing a calibration, the electrical standards should be connected to 
the section that corresponds to the beginning of the material under test (i.e., the 
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Fig. 5.5 Schematic representation of the experimental setup. The rectangle highlights all the 
effects that are included in the error network (not to scale) 14811 



beginning of the probe rods, section A in Fig. 15.51 ). Unfortunately, the structure of 
commercial three-rod probes does not allow such a calibration because the open 
and matched loads insertion would require detaching the rods from the probe head 
and the exact knowledge of the value of the matching impedance. As a matter of 
fact, SOL standards for three-rod probes are not commercially available, and even 
in-house solutions for calibration fixtures mimicking open and load conditions are 
difficult to be reproduced. 

To overcome such a limitation, an innovative calibration procedure that does not 
require either the open or the load conditions was specifically developed. In fact, 
the proposed calibration method relies only on the short-circuit conditions obtained 
at three different sections along the probe (triple short calibration - TSC) [48]. 

The TSC procedure was tested and validated through preliminary measurements 
on well-referenced materials. Successively, it was applied to sandy soil brought 
at different levels of moisture. Results show that the combination of a frequency- 
dependent dielectric mixing modeling, of a minimization routine, and of the pro- 
posed TSC procedure is the key strategy for obtaining accurate TDR moisture 
content evaluations while reducing the complexity of the calibration phase. Another 
possible alternative solution for allowing traditional SOL calibration could be to de- 
sign and realize a custom-made probe, similarly to the case presented in Sect. 14.31 
However, the TSC procedure presented herein can be readily applied to commer- 
cially available multi-rod probes. 

As aforementioned, to effectively enhance the accuracy in the estimation of 
soil moisture content, the reflection coefficient Sh,meas(/) at section A of the 
probe (see Fig. 15.5b must be evaluated. In fact, it is worth mentioning that, simi- 
larly to the previously-discussed cases, TDR intrinsically measures the reflection 
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coefficient along the whole transmission chain (i.e., starting from the instrument 
output connector, up to the probe end). As a direct consequence, the resulting re- 
sponse is inevitably affected by all the parasitics and losses introduced by the cable, 
connections, section transitions, and probe-head portion. A calibration procedure 
allows reducing the influence of these parasitics. 

As well known, the aim of the calibration procedure is to assess the three scat- 
tering parameters that characterize the 'error network', depicted in Fig. 15.51 which 
models all systematic errors effects between the instrument output connector and 
the section corresponding to the beginning of the MUT yj]. Three short circuits 
were applied at different sections along the probe; in fact, these loads can be easily 
applied at whichever distance from the probe head. 

For this specific configuration, the first short was directly applied to section A, 
whereas the other two shorts were applied at distances l\ and 1% from the probe head, 
respectively. When a short circuit condition is realized at a generic distance / from 
the probe head, the reflection coefficient at section A (Sn scA ) can be expressed as 
follows: 

Sn x , A = -e- 2ipi (5.14) 

where /3 = 27i:/(/Xo £ o) i s me propagation constant in air. The 'error network', 
which is characterized by its scattering parameters Sn e! £12 e! an d S22,e, transforms 
such reflection coefficient at section A in a new one, namely Si i B , measured at the 
TDR connector (section B): 

c2 -2ijil 

Jncc 

(5.15) 



Therefore, if the values of Sn sc B corresponding to three different positions of the 
short circuit (namely, /= 0, I = h, and / = l£) are measured, a system of three 
equations in three unknowns is obtained, from which the three error parameters 
Sue; S\2,e', and S22.E can be retrieved. In particular, after appropriate simplifica- 
tions within the Mathematica™ software, the following expressions are obtained 
for the error network parameters: 



5i 



IE 



-^iisco.B ( s nsci,B-Sn sc2 , B )-Sn sc2 , B (Sn sc0B -Sn SC | B )cot(/3li) 

A 

Snscifi (Sii 8o0)B -5 n )cot(^/ 2 ) 
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(5.17) 



5.4 Moisture Content Measurements through TD/FD Combined Approach 101 

-i( 5 Hsci.B -Su x2£ ) + {Sn KOfi -Sn scl , B )cot(j3h) 



She 



( s n sc0 .n-Sn sc2 , B )cot(pl2) 



(5.18) 



where A is given by 



A = - i ( 5 ii scl . B -5ii sc2B )-(5ii sc0B -5ii sclB )cot(/3/i) + 

+ ( 5 ll S cO,B-"Vsc2.B) COt (/^2) 

Su no , Sit ,„, and Su , „ are the reflection coefficients measured at section B for 
the short circuits located at /= 0, / = l\ , and 1 — 1%, respectively. 

Once the error network is characterized by means of the TSC procedure, the 
measurements performed at section B (Sub) can be corrected so as to obtain the 
corresponding reflection coefficient at the beginning of the rods (5ha): 

s Sll.B-Sll.E (519) 

Sll,BJ22,E + S nE - Si 13S223 

A crucial point of the TSC procedure is the optimal choice of the distances l\ and 

h. 

The impedance along a line presents a periodic behavior, with a period equal to 
half wavelength. This means that at the frequency for which fa is equal to half wave- 
length, the distal short appears again as a short at section A, causing a singularity in 
the error network parameters (as a matter of fact, only two of the three loads remain 
independent). On the other hand, choosing very short distances for the three shorts 
would compromise measurement accuracy at low frequencies, because the phase 
differences between the shorts tend to vanish. 

Since the TDR unit used for these measurements has an operating bandwidth of 
approximately 1.5 GHz, a distance fe equal to 10 cm was chosen, yielding a singu- 
larity exactly at the upper bound of the bandwidth. For symmetry, 1\ was set equal 
to 5 cm. 

5.4.2.1 Evaluation of Moisture Content 

Soil moisture content was evaluated through the minimization of the squared mag- 
nitude of the difference between the 5h.mod(/) modeled through ( 15.121 1 and ( 15.131 ). 
and the Sh.measC/") that was derived from TDR measurements and referred to the 
corresponding calibration section (i.e., section A in Fig. 15.51 ). The measurement 
cell is modeled as the portion of a transmission line, and the soil is characterized 
by means of an appropriate dielectric model. Subsequently, the reflection scatter- 
ing parameter of the implemented model, Sh.mod(/)> 1S evaluated as a function 
of the unknown value 6. Finally, through the aforementioned minimization routine 
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(in which Su mbas(/) 1S derived from the TD/FD combined approach), it is possible 
to retrieve the unknown value of moisture content. To carry out this minimization 
procedure, all the parameters in ( 15.13b . with the exception of 6, should be known. 

For this purpose, the sand porosity was preliminarily assessed through the so- 
called water saturation methoqj. 

As for the complex permittivity of water, its frequency-dependent behavior was 
obtained through a Debye model; in turn, the dependence of the Debye parameters 
on temperature was taken into account through empirical expressions reported in 
the literature 03 511 . Static electrical conductivity was measured through the well- 
known Castiglione-Shouse method [8], which allows the estimation of Oq from TDR 
measurements (in particular, the same probe was used). 

Finally, the two remaining unknowns (namely, e sp and a) were estimated as fol- 
lows: two measurements were performed, on dry sand (6 = 0) and on sand moist- 
ened at a known 6 value (namely, 25%); successively, e sp and a were extracted 
through a preliminary minimization routine (which was performed over e sp and a, 
and considering 6 as a known value). On a side note, it is worth mentioning that the 
value of e sp might also be directly measured through the TDR technique using, for 
example, the procedure reported in B5411 . 

5.4.2.2 Preliminary Validation of the TSC Procedure 

The proposed TSC procedure was preliminarily tested on bi-distilled water and 
ethanol. To this end, the three TDR measurements on the short circuits were used 
to derive the error network parameters, according to ( 15. 16b - (f5TT8l ). Successively, the 
error correction procedure, based on ( 15.191 1, was applied to Sn mbas(/) measured 
at section B for the two liquids. Fig. |5.6(a)| shows, for the case of bi-distilled wa- 
ter, a comparison among the absolute values of the Sn(f) at section B, obtained 



through the TD/FD combined approach 111 111 : the corrected Sn at section A; and 



the theoretical S\ \ behavior obtained through TL modeling of the rods immersed 
in water (characterized through its Debye model). Similarly, Fig. |5.6(b)| shows the 
comparison for the case of ethanol. 

In the TL modeling, for the probe rods, an effective length of 152 mm (including 
the fringing effect at the rods distal end) was always considered. 

The analyzed frequency band was limited to 500 MHz, because the three-rod 
probe configuration and dimensions do not allow accurate and repeatable measure- 
ments above such frequency. The figures clearly show the effect of the adopted 
calibration and of the error correction procedures, yielding corrected reflection co- 
efficients curves that are in close agreement with the ones expected from theory. 



A known volume of water (v w ) was added to a known volume of dry sand (v s ); the volume 
of the sand included the pores. After the sand was saturated with water, the total volume of 
the mixture (v m i x ) was measured. Eventually, the volume fraction of the pores was obtained 
as (v w + v s ) — v m j x , from which the sand porosity could be easily derived. 
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Fig. 5.6 Results for the preliminary test of the TSC procedure: comparison of the magnitudes 
of the measured, of the corrected, and of the modeled S\ i . a for bi-distilled water at 25.6°C. 
bforethanolat25.4°c(I! 



5.4.2.3 TSC Application Test-Case: Soil Moisture Estimation 

Construction sand was used as test material; the sample was placed inside a cylin- 
drical holder (with a diameter of 7.8 cm and a height of 25 cm). Preliminary mea- 
surements confirmed that boundary effects due to the used holder were negligible. 

To verify the suitability of the proposed method for the estimation of 0, the sand 
sample was moistened at five prefixed water contents levels (0), corresponding to 
5%, 10%, 15%, 20%, and 25%. The probe was inserted in the sand and was not 
removed until all the experiments were carried out. Moistening was performed by 
adding the fixed amount of water from the top-surface of the sample; when the dis- 
tribution of water appeared homogeneous (which took approximately 30 hours for 
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each moisture level) the TDR measurements could be carried out. To avoid varia- 
tions in moisture profile along the probe, due to gravity effects, the cylindrical holder 
was kept horizontal during all measurements. The sample temperature was kept at 
27 ± \°C during all the experimental session. 

For each moisture level, two TDR waveforms were acquired: one with a very long 
acquisition window and the other with a shorter one. The former was used for the 
evaluation of Oq through the Castiglione-Shouse method Q8[| ; whereas the latter was 
used both for the estimation of the dielectric constant through the traditional TDR 
method [10] and for obtaining ShmeasC/O to be used in the subsequent minimiza- 
tion procedure. For the sake of example, Fig. 15.71 shows a typical TDR waveform 
(corresponding to the dry sand sample), acquired with a short acquisition window. 
Finally, to perform the aforementioned TSC, the TDR waveforms corresponding 
to three short-circuit conditions obtained at different sections of the probe were 
acquired. The short-circuit conditions were achieved by short-circuiting the probe 
rods (at the desired section) by means of a perforated copper slab with a 'guillotine 
configuration' (see Fig. l5.8t . All the TDR waveforms measured at section B were 
preliminarily FFT-processed and corrected by means of ( 15.191 ), thus obtaining the 
Sh.meas(/) at section A to be used in the minimization procedure. 
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As aforementioned, the minimization procedure was initially applied to measure- 
ments on dry sand and on sand with = 25%, in order to estimate the values of e sp 
and a. The optimized values for such parameters, retrieved through a simplex opti- 
mization procedure, are £ sp = 4.70 and a = 0.624. Moreover, sand porosity (j) was 
measured to be 0.39. 

Successively, the minimization procedure was applied to Siimeas(/) f° r eacn 
moistened soil sample, in order to retrieve the optimized value for the volumetric 
water content 9. 

For each moistening step, the obtained results are summarized in Table l5.2l which 
reports the nominal moisture level, the moisture level estimated through the mini- 
mization procedure, and the evaluated static conductivity. Results in Table l5.2l show 
an excellent agreement between estimated and nominal moisture levels for high wa- 
ter content, with a slightly decreasing performance for low-humidity conditions. 
This last effect is attributable to the fact that low-humidity samples exhibit low per- 
mittivity values and would require longer probes to provide more electrical length 
for the TDR signal. 



Table 5.2 Reference volumetric water content (0 re f) and estimated volumetric water content 
(0est), for five moisture levels. The corresponding values of the static electrical conductivity 
(<Tn), evaluated through the Castiglione-Shouse method, are also reported M48I1 
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Fig. |5.9| shows a comparison between ISi^measC/OI an d |Sh.mod(/)|> after opti- 
mization, for 6 = 5% [Fig. |5.9(a)[ and for 6 = 20% [Fig. |5. 9(b)) . A good agreement 
can be observed between the experimental and analytical curves over all the consid- 
ered frequency range, with the exception of the low frequency band. Indeed, a fine 
tuning of the acquisition window length should suitably lead to a balance of some 
contrasting effects, thus preserving the accuracy of the frequency-domain data. In 
particular, as will be detailed in the following chapter, the optimal time windowing 
should balance 1) the frequency resolution; 2) the signal to noise ratio (SNR); and 
3) the inclusion of a sufficient number of TDR multiple reflections so as to obtain a 
complete spectral picture. 

Fig. 15.91 also clearly shows the higher noise in the experimental curves and 
the consequent worsening in the model fitting for lower humidity levels. Finally, 
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Fig. 5.9 Results for the described moisture evaluation method: comparison between 
|5'll,MEAs(/lJ_ an d |Sii,mod(/)|, after optimization, for the sand sample, a at 9 = 5%. b 
at = 20% \4$\ 



Fig. 15.101 shows the permittivity spectra of the sand sample, for = 20%, as ob- 
tained from (15.13b after optimization. It can be seen that the imaginary part of the 
permittivity of the soil is not negligible in the low-frequency range and that it shows 
a markedly dispersive behavior. 

On a side note, it is worth noting that for 6 = 20%, the apparent dielectric con- 
stant (e ap p) directly evaluated from the time-domain TDR traces is 1 1 .83. This result 
is due to the fact that e app incorporates the effect of both the real and the imag- 
inary parts of the dielectric permittivity (e* c ) in the frequency range of the test 
signal Q. 
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Fig. 5.10 Permittivity spectra, as obtained from the adopted dielectric mixing model after op- 
timization, for the sand sample moistened at 6 = 20%. a real part of permittivity, b imaginary 
part of permittivity 14811 



5.5 Enhancements of TDR-Based Static Electrical Conductivity 
Measurement 



Static electrical conductivity (Gq) is a fundamental physical parameter that is ex- 
tremely useful in many research and industrial areas. Two major fields in which 
the estimation of Gq is particularly important are soil science and the agro-food in- 
dustry: in the former, Co can be associated with several different properties such as 
salinity, solute transport, etc.; whereas in the latter, Gq may be used as an indicator 
of the quality status of the products. 
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However, the evaluation of the Go of soils has proved to be a difficult task; and 
this is the case of granular materials in general 11511 . The well-established four- 
point technique for electrical conductivity/resistivity measurement^] imposes strin- 
gent operating requirements that are not always easily achievable, especially when 
dealing with granular moistened material^]. 

For this reason, there has been a growing interest toward alternative measurement 
techniques, such as the Sigma probe 12511 . the electromagnetic induction method 
H, and TDR fig]. 

Indeed, the most attractive trait of the TDR-based approach is that the test-signal 
is a broadband electromagnetic signal: the value of the static electrical conductivity 
is extrapolated through appropriate analysis and processing of the TDR waveforms; 
hence, the DC-related limitations are overcome. It goes without saying that the addi- 
tional advantages of using TDR relate to the possibility to simultaneously estimate 
several other material properties, such as dielectric permittivity llal . water content 
I17I1 . salinity 15 311 . and solute transport 13611 . Finally, the wide availability of portable 
low-cost TDR instruments and the capability of obtaining accurate real-time, in situ 
monitoring make the TDR technique definitely attractive also in view of practical 
applications. 

However, the accuracy of traditional TDR-based electrical conductivity measure- 
ments is strongly influenced by the preliminary calibration technique. In this regard, 
two innovative approaches for speeding up the preliminary calibration procedures 
required for the TDR-based measurement of Gq were assessed. The first method 
relies on the combination of the TDR measurement with a transmission-line mod- 
eling (TLM); the second method simply relies on a couple of independent capac- 
itance measurements (ICM) performed through an LCR meter. The metrological 
performance of this last method is compared with that of the traditional one, thus 
validating its applicability. Experimental results and related uncertainty analysis on 
various samples demonstrate that the proposed method is suitable for a simpler and 
accurate estimation of the ao, also when dealing with moistened soils. 

5.5.1 Traditional TDR-Based Static Electrical Conductivity 
Measurement 

Traditionally, TDR-based measurements of Gq of granular materials are carried out 
employing multi-rod probes. 



In the four-point technique, a current is forced through the two amperometric electrodes, 
placed at the ends of the sample, and the corresponding voltage value is measured at the 
two voltmetric electrodes. 
7 For example, good and electrically-stable contacts between the four electrodes and the 
sample are essential 16711 . Nevertheless, unless very expensive materials are used for the 
electrodes (e.g., titanium, gold), electrochemical processes due, for example, to the pres- 
ence of water in the soil, occur at the electrodes interfaces, thus rapidly degrading the 
electric contacts. Indeed, this kind of measurements should be performed in static condi- 
tions, therefore, as well known, the electrode polarization effect is an unavoidable limit 
that practically impairs the measurement feasibility. 
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TDR-based conductivity measurements rely on the Giese-Tiemann (G-T) method, 
which is based on the linear relation between Gq and the static conductance of the 



load, G s (22D: 

cto = K p G s (5.20) 

where K v is the probe constant. The value of G s is evaluated from measurements 
of the reflection coefficient at longer times (i.e., when the TDR has achieved the 
steady state, approximately corresponding to the zero-frequency response), of the 
probe inserted in the MUT, p^: 

G,g-t= 1 )~ P °° (5.21) 

Ztdr 1 + P~ 

where Ztdr is the output impedance of the used TDR instrument. The value of the 
probe constant can be estimated from the probe geometry: 

^P = ^p (5.22) 

where L is the length of the used probe, Z p is the characteristic impedance of the 
probe, c is the velocity of light in vacuum, and £o is the permittivity of vacuum. In 
( 15.22I I. the evaluation of Z p is not always an easy task. For this reason, to obtain 
more accurate estimation of the electrical conductivity, the proportionality constant 
K v is usually determined empirically, through multiple preliminary calibration mea- 
surements on electrolyte solutions of known conductivity. 

The G-T approach neglects the series resistance of the cable, connectors, and ca- 
ble tester; as a result, for higher electrical conductivit y va lues (i.e., Go > 0.2 Sm -1 ), 



this method underestimates Gq. For this reason, in I28H . the series resistor model 
was introduced, thus considering the coaxial cable and the sample as two resistors 
in series. The G-T method was modified as follows: 

(1/G S )— Rcable 

where -R C abie 1S the resistance due to the cable. Castiglione-Shouse (C-S) proposed an 
alternative approach for taking into account cable resistance; their method required 
the correction of the measured reflection coefficient by scaling it with respect to the 
reflection coefficients measured with the probe in air, Poo.air, and short-circuited at 
the distal end, p^.sc, respectively [8[]: 

~ o P ~ P°°,air . , ,~ 0/1 n 

P-,scaled = 2 hi. (5.24) 

Poo. air Poo,sc 

However, Lin et al. argued that, despite providing accurate results, the C-S method 



is theoretically incorrect and that the series resistor model is indeed accurate 114 111 



In particular, Lin et al. suggested that the C-S scaling method actually accounts for 
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the errors due to imperfect amplitude calibration (and not for cable resistance) and 
that the C-S method overestimates CTq by a constant rate fl40L kill : this error can 



be compensated for if the probe constant is evaluated every time the measurement 
setup changes (even when the probes are identically manufactured). On the other 
hand, the probe constant in the series resistors model has to be determined only 
once for identical probes, also regardless of the cable length. 

Finally, it is worth mentioning that for highly conductive media, because of the 
strong attenuation of the signal, it becomes difficult to measure simultaneously the 



apparent dielectric permittivity and the electrical conductivity. In B4411 . it was shown 
that for 0.4 Sm -1 < <Jq < 0.75 Sm -1 , a possible solution could be to insulate the 
rod-probes with a high-dielectric coating. 

5.5.2 Innovative Calibration Strategies: The TLM and the ICM 
Methods 

It is apparent from ( 15.20b and ( 15. 2U that, once K p is determined for the specific 
used experimental setup, Go can be easily evaluated through a single TDR mea- 
surement on the considered sample. Therefore, the crucial step of this measurement 
method, and hence the related metrological performance, is strictly related to the 
accurate evaluation of K p . Since the estimation of K v is associated to a preliminary 
calibration, it goes without saying that the minimization of the systematic errors 
contributions is strictly related to this procedure; consequently, the evaluation of K v 
strongly influences the resulting accuracy in the estimation of Ob- 

Starting from these considerations, two alternative and innovative approaches for 



an accurate and simple evaluation of K p were assessed IU3I1 . Both these methods 
require only a couple of preliminary measurements: one in air and the other in a 
reference liquid material, such as water, at known temperature (thus avoiding the 
otherwise necessary calibration with several reference samples of known conduc- 
tivity, as required by the traditional C-S method). 

The first proposed method allows accurately estimating K p through the suitable 
combination of TDR measurements with a specific TLM of the measurement setup. 

On the other hand, the second proposed method leads to the accurate estimation 
of K v through two ICMs performed through an LCR meter. As will be detailed later, 
the ICM method has the great advantage of being very rapid, without compromising 
the final accuracy. 

The proposed approaches were validated through repeated Co measurements: 1) 
on reference electrolyte solutions; 2) on sand samples progressively moistened with 
de-ionized water; and 3) on sand samples progressively moistened with salted water. 
The obtained results are compared with those achieved through the traditional C-S 
method. 

Additionally, the performance comparison between C-S and ICM methods is 
evaluated through a rigorous metrological analysis, which, in turn, provides the mea- 
surement uncertainty in the final estimations of Ob- 
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5.5.2.1 TLM Method 

In the TLM method, the TDR probe surrounded by the material is considered as a 



portion of transmission line in which a TEM mode propagates 1221,1641 . 16911 . Consid- 
ering the proportionality between the static capacitance and the conductance for 
such a line, and taking into account the relation between the static capacitance 
per unit length and the characteristic impedance, the following expression can be 
derived: 

<ro = ^f^. (5.25) 

To apply ( 15.25b , Z p and L have to be accurately determined. 

The probe length might be derived from a separate measurement or from the 
probe datasheet. Nevertheless, the fringing effect at the open end of the probe would 
be neglected. As for the evaluation of Zp, for the adopted three-rod probe, an ana- 
lytical formulation relates the impedance in air to the geometrical parameters of the 
probe [5]. However, this formulation can only give an approximate value of Zp, es- 
sentially because the probe rods are not perfectly parallel, and their distance slightly 
varies along the probe. Therefore, these first approximations would limit the ac- 
curacy of Go measurements. In order to overcome this limitation, a specific TLM 
of the experimental setup was developed [21]. This way, by using a suitably imple- 
mented optimization procedure it is possible to extrapolate the accurate values for Zp 
andL. 

More in detail, the experimental measurement system (TDR instrument, cable, 
and three-rod probe) was modeled through a three-section transmission line, plus 
some lumped elements in order to accurately take into account the parasitic effects 
arising from the TDR instrument-to-cable connection and from the probe-head (see 
Fig. 15.11b . The analysis of measured TDR waveforms exhibited a purely induc- 
tive effect at the TDR instrument-to-cable connection and a more complex behavior 
(approximately resembled by a C-L cell) at the probe-head. The TLM was imple- 
mented within the microwave circuit simulator AWR Microwave Office (MWO), 
which allows computing the frequency dependent reflection coefficient at the input 
section. Subsequently, an inversion processing, based on inverse fast Fourier trans- 
form (IFFT), provides the corresponding TDR response. 
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Fig. 5.11 TLM of the experimental setup yj] 
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The first step of the TLM method involves a couple of TDR measurements, which 
are performed on the probe in air and immersed in de-ionized water at a known tem- 
perature, respectively. Secondly, the same experimental conditions are simulated 
through the TL model described in Fig. 15.111 Therefore, by employing the sim- 
plex optimization procedure (available within MWO), the optimal values of probe 
length, characteristic impedance, and each parasitic element are evaluated, through 
the minimization between measured and simulated TDR waveforms. 

In these simulations, the frequency dependent dielectric parameters of water are 
obtained through the well-known Debye modeling [35]. The optimized waveforms 
(together with the corresponding measured waveforms) are reported in Fig. |5. 12(a)! 
for probe in air, and in Fig. |5.12(b)| for probe in water, respectively. An excellent 
fitting quality is observed. The optimized values obtained for Zp and L are 172 Q. and 
15.2 cm, respectively. It is worth noting that the nominal probe length is 15 cm and 
that the approximate characteristic impedance is 180 Q: this value was calculated 
from the nominal geometrical parameters. The slightly larger length obtained from 
simulations is attributed to the aforementioned fringing effect of the open-ended 
termination. As for the optimized value obtained for the characteristic impedance, 
in the simulation, the non-uniform geometry of the probe is appropriately considered 
150. The obtained values of Zp and L led to a value of the probe constant, K p , equal 
to 3.1m" 1 . 

Once K v is estimated, <7o can be evaluated through (15.25b . measuring the static 
conductance G s of the probe immersed in the MUT In turn, G s is related to the 
measured static reflection coefficient: in the TLM method, the ohmic losses of the 



connecting cable are compensated for through the series resistor-model B28H . Fi- 
nally, non-idealities in the TDR instrument are corrected by dividing all the values 
of the measured static reflection coefficients by the one obtained with the output 
connector left open, whose corresponding ideal value is +1. The overall adopted 
procedure can be summarized in the following steps: 

1. the static reflection coefficient, p oc , with the TDR instrument output connector 
left open is measured. 

2. The parasitic resistance, Rp, of the connecting cable is evaluated from the static 
reflection coefficient, p sc (appropriately corrected through p oc ) measured with the 
probe in air short-circuited at its distal end, i.e.: 

R p =Zo—^. (5.26) 

Poc 

The resistance of the probe rods does not appear in d5.26l l since it is negligible 
compared to the resistance of the cable. 

3. The static reflection coefficient with the probe inserted in the MUT (p sam pie)> is 
measured, and the corrected reflection coefficient p CO rr,TLM is derived, i.e.: 

Psample , co m 

Pcorr.TLM — • (•>•*') 
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Fig. 5.12 Comparison between measured and simulated TDR waveforms for the probe, a in 
air. b in de-ionized water 11311 



4. The measured static resistance, .R m eas, is evaluated by 

1 + Pcorr,TLM 



R» 



Zo 



1-P. 



(5.28) 



corr.TLM 



5. The probe static resistance is obtained subtracting R p from /? mea s 

6. Once the static resistance and, hence, the conductance, are evaluated, the static 
conductivity is finally obtained from (15.25b . 
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It is important to emphasize that the proposed TLM method can virtually take into 
account every single parasitic contribution associated to a specific experimental 
setup. As a result, the influence of systematic error contributions is dramatically 
reduced. 

5.5.2.2 ICM Method 

Considering a TEM mode propagating through the three-rod probe, the characteris- 
tic impedance of the probe in air may be written as follows: 



1 
cCo 



Z P = — (5.29) 



where Co is the capacitance per unit length of the probe in air. Combining ( 15.25b 
and d5.29l l, the following expression is obtained: 



G s £q 



<7o = tF = G s K p (5.30) 



where C s o = CqL is the capacitance of the probe in air. 

When the probe is inserted in a sample material, the quasi-static capacitance of 
the probe will include two contributions: 1) a parasitic capacitance, C p , that is due to 
the connecting coaxial cable and to the probe-head; 2) the static capacitance of the 
probe filled with the MUT, C s . Therefore, the total capacitance measured through an 
LCR meter, C m , will be given by the following expression: 

C m = C p + C s = C p + C s o£s (5.3 1) 

where £ s is the static permittivity of the material filling the probe. 

Consequently, from ICMs on two materials with well-known static dielectric per- 
mittivity (e s i and e S 2, respectively), it is possible to retrieve the value of C p and, most 
importantly, of C s o, according to 

Cso = c m2 -c ml (532) 

£ S 2 — £sl 

where C m \ and C m 2 are the LCR-measured capacitances in the two reference mate- 
rials, respectively. Finally, from J5.30I ). it is possible to estimate K p . 

For the ICM method, the Agilent 4263B LCR meter was used directly connected 
to the TDR probe; and air [66] and pure water B24I1 were chosen as reference materi- 
als. LCR measurements were performed at a frequency of 100 kHz, with a measure- 
ment average of 256 and temperature kept constant at 25°C. LCR measurements 
showed that, at this frequency, the probe behaves just like a capacitor; in fact, the 
phase angle of the corresponding impedance resulted approximately -90°, thus con- 
firming the validity of the approximations introduced using ( 15.3 ll rl. 



8 It is worth mentioning that, at 100 kHz, the dielectric behavior of air and water is practi- 
cally static, thus dissipative and dispersive effects are avoided. 
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The advantage of the ICM method is that, thanks to the simplicity of the used 
equations, it leads very rapidly and without much effort to the estimation of K p . 



5.5.3 Validation of the Methods 

The experimental setup for the refiectometry measurements includes the HL1500 
TDR unit 0], the 15 cm-long three-rod probe (Campbell Scientific CS630) H, and 
a 3.5 m-long 50 Q-matched coaxial cable. 

The first step was the evaluation of the probe constant Xp,c-s to be used for sub- 
sequent measurements based on the C-S method. After preparing nine electrolyte 
solutions, the corresponding values of o"b. re f were measured through a conductiv- 
ity meter and considered as reference values. The value of K p c-s was evaluated 
from the linear fitting of the measurement points in the G s — o"b plane, as shown 
in Fig. 15.131 The obtained value for K p cs is 2.9m" 1 , whereas the correspond- 
ing standard uncertainty evaluated according to the formulation described in 01311 is 
UKp = 0.2m . 
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Fig. 5.13 G s — Oo.ref regression curve obtained from the measurements on nine electrolyte 
solutions differing in reference electrical conductivities. The curve slope gives the A" p to be 
used for the estimation of o"b through the C-S method n 1 3M 



At the same time, TDR waveforms at longer distances were acquired on the same 
electrolyte solutions and were used to calculate G s for each electrolyte solution, 
as the average value of TDR measurements in the considered long-distance range. 
Fig. 15.141 shows the typical TDR waveforms measured for the nine electrolyte so- 
lutions, which differed in reference values of electrical conductivity. It can be seen 
that the waveforms are quite stable in the distance range 750-1250m; therefore, 



116 



5 Qualitative Characterization of Granular Materials and Moisture Measurements 



l.o 



■g 

3 0.4 



S 0.2 

S 

Q 

a 
8 o.o 



-0.2 
-0.4 

















\ OmS/m 




begitniLQ 


gof the probe ; !-- 


/ ; 12 mS/m 




;^/-/-fKrtta 








^/ / ; 34 mS/m 






/ -, - _ 














J 




_, 'y 54-mSVm 








V 

1 


end of theprobe ml/ ~*—S-*< 






_^/ i 81 mS/m 












_^/' ; ft" mS/m 












^^^^■LiLL*^— y ■•; 





4.5 



5.0 



5.5 6.0 6.5 7.0 

apparent distance [m] 



7.5 



8.0 







o.s 


\ niS m 


1 12 mS/m 


0.6 




■ 25 mS,'m 


0.4 




:34mS/m 


0.2 








: 54 mS/m 


n r. 










i»lmS/m 


0.2 




T 


0.4 


;99mS/m 
stable steady state 



700 1200 1700 

apparent distance Im] 



Fig. 5.14 TDR waveforms acquired at longer distances on nine electrolyte solutions differing 
in reference electrical conductivities (left). The behavior of the same waveforms at the steady- 
state condition is also reported (right) 11311 



the static reflection coefficients and the corresponding evaluated uncertainties were 
referred to this range. 

Successively, the values of the electrical conductivities (Co.tlm) of the electrolyte 
solutions were also evaluated through the TLM method using, however, the same 
TDR waveforms. To this purpose, preliminary TDR measurements with the probe 
in air and in pure water were performed in order to obtain the optimized values for 
the TLM model parameters, as discussed in the previous sections. 

Furthermore, considering the TDR-probe in air and in pure water, and performing 
capacitance measurements, the probe constant deriving from the ICM method was 
estimated, as detailed in the previous sections. Results are summarized in Table 15.31 
in which the measured C m , and the estimated K p icm are reported along with the 
corresponding absolute standard uncertainties. 

The electrical conductivities deriving from the ICM method, (Jo.icm, were evalu- 
ated for the electrolyte solutions. 

Table l5.4l shows the comparison between the experimental results for Oo.tlm an( ^ 
<7o,icm> obtained on the reference electrolyte solutions. The good agreement among 
results provides a useful preliminary validation of the two considered approaches. 
After the preliminary validation of the method on reference electrolyte solutions, 
the described methods were used to measure the electrical conductivity of moist- 
ened sand. To this purpose, sand samples were watered at ten different moisture 
levels (ten samples were moistened with de-ionized water and ten with salted wa- 
ter), thus obtaining a wide range of experimental conditions. TDR waveforms at 
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Table 5.3 Results for LCR capacitance measurements on pure water and air. The correspond- 
ing literature values of static permittivities and the estimated K p value are also reported, with 
associated standard uncertainties 



material 


t- m 

(pF) 


£s 


^"p.ICM 

(m- 1 ) 


air 


294.07 ±2.85 


1.0005 ±0.0002 


3 1 ±0 1 


water 


517.13±4.27 


78.5 ±0.1 





Table 5.4 Comparison between static electrical conductivities of electrolyte solutions (at 
different CTo.ref) evaluated through the TLM (Ootlm) an d the ICM (Oqicm) methods. The 
absolute standard uncertainty (only for the ICM method) is also reported 



C0,TLM 

(mSm" 1 ) 


0h,ICM 

(mSm -1 ) 


1.2 


1.2±0.1 


12.6 


12.3±0.3 


23.8 


23.3±0.6 


35.0 


34.2±0.8 


45.2 


44.3±1.1 


56.7 


55.5±1.3 


68.8 


67.4±1.6 


89.3 


87.5±2.1 


105.8 


103.6±2.5 



longer distances were acquired on all the prepared samples, and the corresponding 
electrical conductivity values (Ob.c-S, Oo,tlm> an d Ou.icm ) were evaluated. 

Table 15.51 and Table 15 . 6l summarize the results obtained for the o"b measurements 
(through the C-S, the TLM, and the ICM methods), on sand moistened with de- 
ionized water and on sand moistened with salted water, respectively. The obtained 
results show an overall good agreement among the values obtained through the three 
methods. In particular, considering the sand samples, it is worth noting that the 
uncertainty related to the ICM method is always lower than the uncertainty related 
to the C-S method: this is particularly evident for higher moisture levels. 

It is important to note that the final uncertainty on <7o,c-s an d Co.icm 1S mostly due 
to the uncertainty in the evaluation of the probe constant through the C-S method 
and through the ICM method, respectively. In fact, the percentage standard uncer- 
tainties are 7.5% and 2.3%, for K p cs and K p jqm > respectively. Conversely, the 
contributions of the uncertainty related to the other quantities involved (i.e., p sam pie, 
Poc, Psc) hardly influences the final uncertainty of the calculated value of o"b- in fact, 
the maximum values of the percentage standard uncertainty for static reflection co- 
efficients, not reported here for the sake of brevity, are in the order of 0. 15%. 
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Table 5.5 Results obtained for measurements on sand samples mixed with deionized water, 
at ten different moisture levels. The estimated values for Oq c-Si Co.tlm; ar| d °0,ICM an d for 
the corresponding standard uncertainties are compared 



e 

(%) 


Ob.C-S 
(mSm -1 ) 


0U,TLM 

(mSm -1 ) 


0h,ICM 
(mSm -1 ) 


0.0 


- 


0.0 


- 


2.0 


0.6±0.1 


0.8 


0.7±0.1 


4.0 


1.4±0.1 


1.6 


1.6±0.1 


6.0 


2.6±0.2 


2.9 


2.9±0.1 


9.1 


3.6±0.3 


4.1 


4.0±0.1 


12.1 


4.2±0.3 


4.7 


4.6±0.1 


15.1 


5.7±0.4 


6.4 


6.3±0.1 


18.1 


7.3±0.6 


8.1 


8.0±0.2 


21.1 


11.5±0.9 


12.7 


12.4±0.3 


24.4 


15.7± 1.2 


17.3 


16.9±0.4 



Table 5.6 Results obtained for measurements on sand samples mixed with salted water, at 
ten different moisture levels. The estimated values for <7o,C-S> Oo.tlm ar, d Oo.icm an d f° r the 
corresponding standard uncertainties are compared 



d 

(%) 


0b.c-s 
(mSnr 1 ) 


OiD,TLM 

(mSm" 1 ) 


0h,ICM 

(mSm" 1 ) 


0.0 


- 


0.0 


- 


2.0 


0.8±0.1 


1.0 


1.0±0.1 


4.0 


2.5±0.2 


2.8 


2.7±0.1 


6.0 


4.1±0.3 


4.6 


4.5±0.1 


9.1 


6.3±0.5 


7.0 


6.9±0.2 


12.1 


9.9±0.8 


11.0 


10.8±0.3 


15.1 


14.3±1.1 


15.8 


15.6±0.4 


18.1 


20.2±1.5 


22.3 


21.8±0.5 


21.1 


30.7±2.3 


33.8 


33.0±0.8 


24.5 


41.5±3.1 


45.6 


44.6±1.1 



Finally, it must be pointed out that the K^cs was evaluated through nine pre- 
liminary measurements on the electrolyte solutions, and it was found that MKp.c-s 
increases dramatically when less than nine points were considered for the linear re- 
gression. Conversely, only two preliminary measurements were enough to obtain an 
accurate estimation of the probe constant through the ICM method. This is one of 
the major advantages of the described methods. 
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5.5.4 Practical Considerations 

Preliminary results obtained comparing the traditional and the described methods 
highlighted similar performance in terms of estimation of the electrical static con- 
ductivity, both for the electrolyte solutions and for the moistened sand samples. 
However, it is important to underline that the two proposed approaches are advanta- 
geous in two major ways. In fact, the TLM method clearly discovers and compen- 
sates for the major parasitic effects that influence the final measurement accuracy. 
As for the ICM method, the simple requirement of two capacitance measurements 
simplifies and speeds up the entire procedure. Additionally, the addressed metro- 
logical validation clearly demonstrates that the uncertainty in the estimation of elec- 
trical conductivity through the ICM method is always smaller than (or, at worst, 
equal to) the uncertainty related to the traditional C-S method. On such basis, the 
described approaches represent a viable way toward a cost-effective and time-saving 
TDR monitoring of conductivity for quality control of both liquid and granular 
materials. 

5.6 Noninvasive Moisture Content Measurements 

In some monitoring or diagnostic applications, the noninvasive approach may rep- 
resent a stringent requirement. Therefore, with regards to the BMR capabilities to 
satisfy such needs, much research effort has been paid to designing and realizing 
probes that could guarantee a noninvasive approach. For example, metallic strips 
have been investigated as possible surface-TDR probes [58]. Similarly, serpentine- 
shaped wires (attached to a planar acrylic pad) have also been considered as a non- 
invasive solution for BMR-based moisture monitoring [59]. 

Starting from these considerations, in this section, the possibility of noninva- 
sive moisture monitoring through a suitable combination of BMR with the use of 
a microstrip antenna as sensing element is addressed. This method may become 
particularly suitable for monitoring applications in which the noninvasive approach 
is mandatory (e.g., construction materials, walls, architectonic structures, frescoes, 
etc.). 

As will be detailed in the following subsections, the dielectric properties of the 
considered material are associated to the changes in the reflection scattering param- 
eter of the antenna, Su (/), occurring when the antenna is placed in contact with the 
surface of the MUT 1570. More specifically, the changes in the resonant frequency 
(/res) of the antenna can be associated to the moisture content of the MUT. 

The main steps of the adopted method can be summarized as follows: 

• a microstrip antenna (which can be specifically designed and realized) is placed 
on the surface of the MUT, and the TDR waveform of the system (antenna and 
MUT) is acquired [see Fig. 15.15ft ; 

• the Si i (/) of the system is evaluated through the TD/FD combined approach 
described in Chapt. [3] 0] ; 
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Fig. 5.15 Schematic diagram of the experimental setup used for the noninvasive measure- 
ments 11211 



• the S\ i (/) is used to evaluate the resonant frequency of the system (occurring at 
the minimum of the Su(f) curve); and finally, 

• the evaluated / res is associated to the moisture level of the sample. 

As well known, changes in the resonant frequency of the antenna could be directly 
determined through a vector network analyzer (VNA). However, as already men- 
tioned, using a TDR instrument in place of a VNA can reduce the costs related 
to the experimental setup. In this regard, as also mentioned in previous cases, al- 
though the reported results were obtained through the TDR80E04 16 ill , a benchtop 
and expensive instrument, the proposed method can be equally implemented using 
affordable TDR instruments. 



5.6.1 Basic Theory of Microstrip Antennas 

Before proceeding with the detailed description of the experimental procedure, 
in this subsection the basic theoretical aspects behind the proposed method are 
addressed. 

A microstrip antenna consists of a dielectric substrate with a radiating patch on 
one side and a ground plane on the other side [see Fig.s |5.16(a)| and r5.16(b)| . This 
kind of antennas is inexpensive to be fabricated, and a high level of reproducibil- 
ity can be achieved. As well known, an antenna is characterized by its resonant 
frequency, which is defined by the following equation 13711 : 



/re 



2i e ff\/^eff 



(5.33) 



where c = 3 x 10 8 ms _1 is the speed of light in free space, L^g is the effective length 
of the radiating patch, and e e ff is the effective relative permittivity of the antenna. 
L e ff and §.g can be evaluated considering the following equations: 



L eff = L + 2AL = L 



2/z 



(5.34) 
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Fig. 5.16 a geometric configuration of the microstrip antenna, b photograph of the used 
antenna I12fl 



£eff : 



£r.sub ' 



£r,sub 



1 



10- 



-1/2 



(5.35) 



where e rsu b is the relative permittivity of the antenna substrate; h is its thickness; 
w and L are the width and the length of the radiating patch, respectively; and £ m is 
the relative permittivity of the medium in which the antenna radiates [see Fig. 15.161 



1651 II 



Considering the previous formulas, when the antenna radiates in free space, then 
£ m equals 1 (which is the relative permittivity of air), and the antenna will resonate 
at a specific frequency, f®£. On the other hand, when a different material is placed in 
front of the antenna, £ m changes, and / res will change accordingly. The direct impli- 
cation is that the resonant frequency is closely related to the dielectric characteristics 
of the medium in which the antenna radiates. 

As for the practical evaluation of / res , it can be measured from the minimum of 
the magnitude of Sn (/), which is given by 
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n ~~ Zant 



|Sn 



ri+Z 1 



ANT 



(5.36) 



where r\ — y/no/e m and Zant is the antenna impedance. Equation ( 15.361 ) shows that 
Sn(f) of the antenna is intrinsically associated to £ m . 

It is important to underline that in the previous equations the thickness of the 
material sample (d) in front of the antenna is assumed as semi-infinite. This as- 
sumption remains valid as long as the power fed into the antenna is sufficiently 
low. In this way, the antenna is not affected by what is beyond the back of the 
MUT (i.e., only the reflections at the front-face of the MUT have to be considered) 



[see Fig. 15. 151 . 

It is worth pointing out that Eqs. ( 15.33b - (15.35b represent a simplified model 
for describing the dependence of the resonant frequency on the relative dielectric 
permittivity. Indeed, a more accurate model should also take into account the effect 
of dielectric losses (by considering the complex relative permittivity of the material 
in which the antenna radiates) and the thickness of the MUT |6J]. Nevertheless, as 
confirmed by the experiments reported in the following subsection, ( 15.33b - ( 15.35b 
provide sound qualitative information on the dependence of / les on e m . 



5.6.2 Experimental Validation of the Method: Measurements on 
Moistened Sand Samples 

The used antenna is a simple and inexpensive A/2-resonant microstrip antenna re- 
alized on an RF35 substrate with a coaxial/transmission-line excitation [see Fig. 



15.1611 ; its working frequency (in air) is approximately 2 GHz. 

The antenna was fed through an SMA connector, thus allowing performing the 
SOL calibration measurements as required by the TD/FD combined approach II ill . 

5.6.2.1 Validation of the Method on Reference Materials 

Preliminary measurements with the antenna used as probe were performed on seven 
well-referenced materials: air, butan-1-ol, n-hexane, chlorobenzene, 2-heptanol, 2- 
propanol, and iso-propoxyethanol. This was done for validating the theoretical 
model described by Eqs. d5.33b - d5.35b . thus verifying the direct relation between 
£ m and the corresponding / res . 

The reference liquids were placed inside the box; the antenna was placed in con- 
tact with each of them, and the corresponding resonant frequency (/ res -ref) was mea- 
sured. The frequency-dependent permittivities of the reference liquids were mea- 
sured according to the procedure described in Chapt.[3]and compared with reference 
values reported in Hand [24|]. A good agreement between measured and theoretical 
data was observed 01211 . 
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5.6.2.2 Measurements on Sand Samples 

After the validation of the reference-model, the proposed method was tested on three 
different types of construction sand, differing in granulometry (ranging from approx- 
imately 1 mm to 4 mm) and composition (sand #1 and sand #2 were siliceous sands, 
whereas sand #3 was a calcareous sand). Moisture measurements were performed 
on sand samples brought at increasing (pre-established) levels of moisture, through 
the gravimetric method. Each gravimetric water content progressively achieved was 
then associated to the corresponding volumetric water content (0) through the mass 
density of the sand. It is important to underline that samples were moistened al- 
most up to water saturation level. In this way, for each type of sand, a wide range 
of reference moisture levels was obtained (referred to as sl,s2, ...,s6, respectively). 

The TDR80E04 was also used for performing 'traditional' comparative TDR 
measurement on the same samples. For this purpose, the apparent relative permit- 
tivity values of the reference samples, £ app , were evaluated (directly in TD) through 
separate measurements of the apparent length (L app ) of the 15 cm-long three-rod 
probe (Campbell CS630) inserted into the MUT. The e app values were evaluated 
through (15.51 [2]. Successively, on the same sample, measurements were carried out 
using the antenna as a probe. For the measurements with the antenna, the sample 
was placed in a top-opened plastic box with dimensions 20 cm x 20 cm x 4 cm. Full- 
wave electromagnetic simulations of the radiation pattern of the antenna, confirmed 
that boundary effects caused by the box were negligible. 

For all the samples (*1,*2, ..., s6) of each sand, the Si\(f) of the antenna in 
contact with the MUT was evaluated. 

As an example, the changes in the S\ \ (/) magnitude curves for increasing values of 
moisture of the MUT for sand #3 are shown in Fig. 15.171 starting from the antenna 
in air, each curve refers to a different moisture level (6) ranging from 0% (si) to the 
maximum considered moisture level (s6). 

To have a tentative value of the relative permittivity of the MUT, traditional com- 
parative TDR measurement of the e app were performed, so as to associate these data 
to the / res measured for the same sample. Results are summarized in Table 15.71 
As expected, as water content is increased, the value of e app progressively differs 
from the value of the theoretical e m . This behavior is mainly attributable to two 
main effects. The first is the fact that dielectric losses become more relevant as the 
amount of water is increased. In fact, the measured e app carries a strong contribu- 
tion of the imaginary part of the relative permittivity (e"), and cannot be merely 
assumed equal to e'. In this case, the theoretical model should include a correction 
factor accounting for e" |6J]. Secondly and more importantly, as discussed in ||6[] 
and as confirmed by measurements, the proposed method is best suited for monitor- 
ing materials whose corresponding relative permittivity does not exceed 10-12. For 
higher permittivity values, due to the intrinsic degradation of the resonant behav- 
ior of microstrip antennas, a strong attenuation of the resonance effect is observed. 
Therefore, the theoretical model described by Eqs. (15.33l l-( l5~3~5T > becomes less ac- 
curate in describing the changes in the resonant frequency as a function of the per- 
mittivity. Deviation between e m and e app are highlighted, for sand #3, in Table [ 
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Table 5.7 Resonant frequencies (/ res ) associated to the pre-established water content level 
(6) corresponding to the apparent permittivity values (e ap p). for each type of sand (#1, #2, 
#3)d 



sand #1 ( granulometry < 1 mm) 


MUT 


6(cm J 


/enr*) 


Eapp 


/ res (GHz) 


si 


0.000 




2.75 


1.958 


s2 


0.031 




3.22 


1.955 


s3 


0.062 




3.94 


1.949 


s4 


0.093 




4.72 


1.944 


s5 


0.124 




5.39 


1.909 


s6 


0.156 




6.55 


1.885 






sand #1 


! (granulometry 2-^4 mm 


) 


MUT 


e(cm J 


/cm 3 ) 


Eapp 


/res (GHz) 


si 


0.000 




2.80 


1.965 


s2 


0.034 




3.13 


1.958 


s3 


0.067 




3.54 


1.949 


s4 


0.101 




4.85 


1.941 


s5 


0.137 




6.29 


1.906 


s6 


0.171 




7.30 


1.906 




sand #3 (granulometry 2 mm) 


MUT 


8 (cm a 


W 


£app 


/res (GHz) 


si 


0.000 




3.91 


1.949 


s2 


0.038 




4.13 


1.947 


s3 


0.076 




4.66 


1.930 


s4 


0.114 




5.69 


1.912 


s5 


0.153 




8.33 


1.891 


s6 


0.191 




10.14 


1.866 



Results reported herein show that the proposed method provides good results for 
volumetric moisture content lower than 20% (e.g., in the considered cases £ app = 
10). Fig. 15.181 shows the final calibration curves relating the resonant frequency 
of each moistened sample to the corresponding moisture level. The curve of sand 
#3 is evidently more linear than the curves referred to sand #1 and sand #2: this 
is attributed to the different granulometric and porosimetric characteristics of the 
considered sands. In fact, the grains of sand #3 have rather uniform dimensions 
(2 mm). Conversely, the grains of sand #1 and of sand #2 have different dimensions 
(in the ranges reported in Table 15 .71 . For this reason, for lower moisture levels sand 
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Fig. 5.17 Magnitude of S\\ (/), for increasing moisture levels (.si being dry sand) for sand 
#3fH 



Table 5.8 Deviation between theoretical values relative permittivity (e m ) and measured ap- 
parent permittivities (e ap p) corresponding to each moisture content level, for sand #3 IU2I1 



moistened 
MUT 


e 

(cm 3 /cm 3 ) 


^m 


A£ — \£m ^appj 


si 


0.000 


5.18 


1.27 


s2 


0.038 


5.39 


1.26 


s3 


0.076 


7.19 


2.53 


s4 


0.114 


9.33 


3.64 


s5 


0.153 


11.8 


3.47 


s6 


0.191 


14.8 


4.66 



#1 and sand #2 show a similar behavior; on the other hand, as the moisture content 
increases, the two curves show a different trend, since granulometry and porosity 
start to affect differently the water absorption behavior. 



5.6.2.3 Effect of Electrical Conductivity 



Considering Fig. 15.171 as the water content is increased, two combined effects can 
be noticed: a shift of / res towards lower frequencies, and a change in the value of 
the minimum of |5i i (/) | . The first is mainly attributable to a change of the real part 
of the permittivity, whereas the variation in the minimum of the Si \ (/) magnitude 
is influenced both by the real and by the imaginary parts of the relative permittivity 
of the watered material. 
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Fig. 5.18 Calibration curves relating the resonant frequency of the antenna (/ res ) to the cor- 
responding moisture content of the MUT (0) for the three types of sand II 1211 



As a further validation of these results and to investigate how the imaginary per- 
mittivity influences the Su(f), full-wave simulations were performed through a 
commercial simulator, CST Microwave Studio. More specifically, the time-domain 
solver, based on the well-known Finite Integration Technique, was used. Simula- 
tions were performed by varying the value of the electrical conductivity of the MUT 
(<7mut), for several fixed values of the real part of the dielectric permittivity (£mut)- 

Fig. 15.191 shows how the minimum of the Su (/) magnitude changes as a result 
of variations of cJmut- It can be seen that when the permittivity is higher than 3, 
an increase of the electric conductivity causes |Sn(/)| to decrease. Simulation re- 
sults in Fig. 15.191 are also compared to measurements performed on sand samples 
moistened with approximately the same amount of bi-distilled water (6 =8%), with 
a conductivity ranging from to 107mSm _1 . The experimental trend is in good 
agreement with the simulations. 



5.6.2.4 Sensing Volume 

Another aspect that is worth taking into account is the sensing volume, which rep- 
resents an issue also for traditional moisture content measurements. For a given an- 
tenna, the sensing volume mainly depends on the relative permittivity of the MUT. 
To assess how the relative permittivity of the MUT affects the sensing volume, full- 
wave simulations were performed through CST Microwave Studio. For different 
fixed values of dielectric permittivity, the thickness of the MUT was progressively 
increased, and the corresponding changes in the resonant frequency of the antenna 
were observed. Results of these simulations are shown in Fig. 15.201 
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Fig. 5.19 Comparison between simulation results (for different fixed permittivity values) and 
experimental measurements performed on 8% -moistened sand #3 with a water conductivity 
ranging from to 107 mSm -1 11211 
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Fig. 5.20 Full- wave simulation results showing the change in the resonant frequency (/ res ) 
as a function of the thickness of the MUT, for several fixed values of £mut 11211 



As a result, for an assumed dielectric permittivity, the predicted sensing 
volume can be associated to the maximum sample thickness beyond which no 
relevant changes in resonant frequency are present. It must be pointed out that 
simulations were performed considering a MUT with fixed transverse dimensions 
(20cmx20cm) and varying thickness; indeed, a more rigorous analysis of the sens- 
ing volume should consider infinite transverse dimensions. However, considering 
the directivity of the used (and simulated) antenna, the sensing volume analysis re- 
ported herein can be considered adequately accurate. 
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Chapter 6 

BMR Characterization of Antennas through the 

Combined TD/FD Approach 



'Science progresses best when obser\'ations force us to alter 
our preconceptions. ' 
Vera Rubin 

Abstract. Antenna characterization measurements are traditionally performed in 
the frequency domain (FD) through vector network analyzers (VNAs) in anechoic 
chambers. Nevertheless, the expensiveness of these facilities limits the possibility of 
using this approach for routine measurements. In this chapter, the strategies for the 
accurate evaluation of the reflection scattering parameter (Si i (/)) of antennas, start- 
ing from simple time domain reflectometry (TDR) measurements, are described. As 
a matter of fact, not only are instruments operating in time domain (TD) usually less 
expensive than VNAs; but they allow performing time windowing, which is the key 
for excluding the unwanted spurious reflection from the environment, thus avoiding 
the use of anechoic chamber. 



6.1 Introduction 

Traditionally, measurements on antennas are carried out in the frequency domain 
(FD) through vector network analyzers (VNAs) (6j], in dedicated facilities (i.e., 
anechoic chamber). Although this procedure undoubtedly provides highly accurate 
results, the high costs involved make such procedure virtually impracticable for rou- 
tine measurements. For this reason, there is an abiding interest for assessing alterna- 
tive methods (less expensive and yet accurate) for antenna characterization. Starting 
from these considerations, this chapter focuses on the possibility of measuring the 
reflection scattering parameter, Sn(/), of antennas through the TD/FD combined 
approach, particularly useful for evaluating the impedance matching performance 
and the resonance characteristics 11711 . 

TDR has been largely used for the characterization of electronic devices, espe- 
cially for characterization at low frequencies where VNA have some intrinsic limi- 
tations 111811 ; yet the characterization of antenna through a TD-based approach is still 
an open issue. 

As a matter of fact, some of the advantages of the VNA (e.g., its high dynamic 
range, frequency resolution and bandwidth) can be achieved with instruments op- 
erating in TD, by appropriately tuning some measurement parameters, such as the 
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time window (T w ), the number of averages, and the sampling period (T c ). More- 
over, when reference standards are not readily available, de-embedding the effect 
of fixtures is easier when using TDR rather than when using a VNA (for example, 
through a suitable time domain gating) |16]. 

It is worth mentioning that many existing works deal with the possibility of com- 
pensating errors in time domain antenna measurements through time gating tech- 
niques j2ll7U8ll9lllOl.llll.li3l. li5n . Nevertheless, no work describes a simple pro- 
cedure to be practically implemented even when basic TD instrumentation (i.e., 
without any specific additional tool) is used. On such basis, an in-depth analysis 
that pinpoints the crucial steps for an accurate TDR-based antenna characterization 
is hereby provided. 

Similarly to the applications discussed previously, the overall procedure includes 
two main steps, namely the TDR waveforms acquisition and the correspondent FD 
signal evaluation. In the following subsections, both steps are discussed in detail. 
Starting from theoretical considerations, it is demonstrated that the optimal choice 
of a specific time window can suitably balance all the effects that limit the final 
measurement accuracy. The optimal trade-off is reached by 

1. optimizing the frequency resolution; 

2. maximizing the signal to noise ratio (SNR); 

3. including as many TDR multiple reflections so as to obtain a complete spectral 
picture; and 

4. excluding the spurious reflections coming from the surrounding objects. 

To assess the proposed method from an experimental point of view, the S\\(f) 
of two antennas was evaluated from TDR measurements performed under differ- 
ent experimental conditions (i.e., choosing different acquisition windows, applying 
digital filters on TD data, placing reflecting objects near the antenna). The crucial 
parameters related to the achievement of a measurement accuracy comparable to 
that provided by the use of a VNA in an anechoic chamber, were individuated and 
optimized. The obtained data were compared to reference- VNA measurements per- 
formed in an anechoic chamber (5 , n ire f(/)), thus definitely validating the accuracy 
enhancement yfl, Q5Q . 

On the basis of the aforementioned discussion, the ultimate goal is to demonstrate 
that TDR-based measurements, in conjunction with a specific data processing, can 
be regarded as a robust and cost-effective procedure for accurate characterization of 
antennas. 

This paves the way for the use of inexpensive TDR devices in practical appli- 
cations (without compromising measurement accuracy and reliability), and for the 
adoption of such a method as the next best thing to VNA measurements in anechoic 
chamber. 



6.2 Measurement Setup for the Validation of the Method 

The first considered device is Alien ALR-8610-AC antenna: a patch antenna in- 
tended for a Radio Frequency Identification (RFId) reader (Fig. |6.1(a)| ). The second 
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Fig. 6.1 a RFId-reader antenna, b biconical antenna |3|] 



is the Clampco Sistemi AP3000 biconical antenna: a wideband biconical antenna 
(Fig. lQM - 

These antennas have very different characteristics in terms of operating frequency 
band; in fact, the former is a narrow-band antenna and operates in the 865 MHz- 
940 MHz frequency range, whereas the latter is a wide-band antenna designed to 
work in a larger frequency range (80MHz-3GHz). The different performance of 
these two antennas allow assessing the proposed approach for different possible 
practical conditions. 

The TDR measurements were performed through the TDR80E04 module. 



6.3 Acquisition in Time Domain 

To obtain the corresponding FD information, TDR waveforms are suitably win- 
dowed and processed through the FFT-based algorithm described in subsections 
I3.4.1l and l3.6.2l It is worth noting that for the TD-acquisition, a rectangular window- 
ing is a wise choice, since it provides good balance between amplitude accuracy and 
frequency resolution for the subsequent FD-transformation. 

Averaged TDR measurements were performed setting the maximum number of 
acquisition points (4000 points), thus forcing the time domain resolution (or sam- 
pling period, T c ) to the minimum allowed by the instrument for each considered 
time window. 

When acquiring a TD waveform for a fixed maximum number of measurement 
points, the only parameter that can be modified is the time window duration. As will 
be clarified, the choice of this parameter has direct consequences in terms of mea- 
surement performance. In order to understand how the time windowing affects the 
final results, let us briefly recall how the TDR instrument performs measurements. 
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As already mentioned, the reflected waveform usually includes the so-called mul- 
tiple reflections due to the signal traveling back and forth between the measurement 
instrument and the device under test. Multiple reflections contain information on 
the spectral behavior of the whole system, including those contributions that in- 
trinsically 'belong' to the device. Therefore, as a general rule, a time window long 
enough to include all the multiple reflections occurring before the steady-state con- 
dition is achieved, should be set. In this way, the subsequent FD-transformation 
would provide a complete spectral picture. 

Additionally, longer time windows guarantee a better frequency resolution, A/, 
of the FD-transformed data, according to the following equation: 

A/=^ (6.1) 

where T w is the duration of the time window. 

However, when characterizing antennas, other issues come along. In fact, when 
the electromagnetic signal generated by the TDR unit enters the antenna connected 
to it, part of the signal is reflected back due to the impedance mismatch of the phys- 
ical device (i.e., at the antenna port); nevertheless, another portion of the signal is 
radiated according to the radiation pattern of the antenna. As soon as the radiated 
signal reaches an object, it is reflected and the antenna receives the signal that is 
bounced back. As a consequence, the measured TDR waveform might include un- 
wanted contributions due to some spurious reflections. Such effects can be excluded 
by choosing an appropriate time window that leaves them out. This is a relevant 
advantage of TD-based measurements over the direct FD-measurements, and can 
be considered as the practical way for avoiding the use of the anechoic chamber. 
Fig. |6.2(a)| shows a typical schematization of a TDR antenna waveform: in par- 
ticular, the grey-colored area corresponds to the portion of the signal that should 
be included in the acquisition window. Fig. |6.2(b)| shows a schematization of the 
measurement setup in which an antenna is connected to the TDR unit through an 
L c Tong cable and the nearest reflecting object is placed at a distance d from the an- 
tenna. Therefore, the time window that excludes possible unwanted reflections can 
be calculated through the following equation: 

7w = 2^ + ^J (6.2) 

where c = 3 x 10 8 ms _1 is the speed of light in vacuum, L is the fraction of L c that 
corresponds to the portion of the cable included in the time window^, and v is ve- 
locity of the signal propagated in the cable. The incident signal portion is windowed 
out, as required for the successive FD-processing [16]. 



1 In practice, the L-long cable portion inclusion is necessary for effectively performing the 
subsequent FD-transformation, as reported in 141, 1 1 411 ■ 
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Fig. 6.2 a Schematization of a typical TDR waveform, b schematization of possible interfer- 
ences from the environment ]5|] 



6.4 RFId Antenna Results 



The configuration of the Alien ALR-8610-AC antenna antenna is shown in Fig. 
|6.1(a)| The truncated edge-geometry of this antenna (highlighted with circles in the 
figure) is realized to guarantee circular polarization with only one feed point and 
to generate two closely-spaced resonant frequencies, between which the antenna 
should operate 11211 . 

The antenna was connected to the TDR module through a 28 cm-long LMR195 
cable with a nominal impedance of 50 Q. It is worth underlining an important as- 
pect on the feed of this antenna under test (AUT). The terminal part of the central 
conductor of the coaxial cable that feeds the antenna is not surrounded either by 
dielectric insulator or by external jacket: this introduces an additional impedance 
mismatch (between the 50 Q-cable and the antenna) that inevitably affects the TDR 
response of the antenna yj] . 



138 6 BMR Characterization of Antennas through the Combined TD/FD Approach 

Before proceeding with the individuation of the most appropriate time window, it 
is useful to analyze the TDR waveform associated to the antenna. Fig. |6.3(aj| shows 
the TDR waveform of the RFId antenna for a lOOns-long time window; Fig. |6.3(b)| 
zooms the 0-15 ns range. 

The first portion of the waveform at approximately 0.25 V corresponds to the 
50 Q-matched cable portion. The measured time between the beginning of the ac- 
quisition window and the end of the cable is 2.8 ns, corresponding to the length L 
in ( 16.2b . After this portion of cable, there is a brief but abrupt change in the ampli- 
tude of the waveform caused by the considerable impedance mismatch introduced 
by the 'naked' termination of the feeding cable (feed pin). The following portion of 
the waveform carries the antenna 'imprint' and, subsequently, the reflected signal 
approaches the steady-state condition (around the value of 0.5 V that corresponds 
to the open circuit) through several multiple reflections, whose shape is related to 
the resonant behavior of the antenna. The attenuating response of the multiple re- 
flections can be clearly distinguished until their peak-to-peak signal excursion is 
overwhelmed by noise. 



6.4.1 Practical Guidelines for Retrieving Accurate Measurements 

In this subsection, the results achieved through the FD-transformation of TDR data 
are reported. 

To have a definitive reference to which TDR-based measurements could be com- 
pared, the Sn re f (/) of the AUT was measured through an HP8753C VNA, equipped 
with the 85047A S-Parameter Test Set, inside an anechoic chamber (Fig. 16.41 ). After 
the short-open-load (SOL) calibration procedure, VNA measurement was carried 
out in the frequency range between 8 10 MHz and 960 MHz (with a frequency reso- 
lution of 1 MHz), since this range corresponds to the operating band of the antenna. 
The S\ i iIB f (/) curves (magnitude and phase) are shown in Fig. |6.5(a")] and Fig. |6.5(b)| 

As for the TD measurements, to avoid undesired spurious reflections, measure- 
ments were performed outdoor, making sure that no reflecting object was in the 
nearby of the antenna. A first set of measurement was carried out choosing differ- 
ent time windows, thus evidencing the consequences in the FD-transformed data. 
Therefore, different TDR waveforms of the AUT were acquired, varying the time 
windows (T w ), from 10 ns to 100 ns, and the S\ i (/) data were evaluated. Results are 
reported in Fig. |6.5(a)| and Fig. |6.5(b)[ for magnitude and phase, respectively, and 
they are compared to the Sn. re f(/). 

The root mean square error (rmse) between the magnitude of 5 , n jre f(/) and the 
magnitude of each of the S\ i (/) evaluated from TDR measurement was calculated 
for qualitatively identifying the most appropriate time window duration: results are 
summarized in Table 16.11 The rmse analysis anticipates that the most appropriate 
windowing seems to be T w — 32 ns. Starting from this experimental validation, in 
the following paragraphs an in-depth analysis is conducted, thus demonstrating how 
the optimal time windowing can be predicted, thus avoiding the use of reference 
VNA measurements. 
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Fig. 6.3 a TDR waveform for the RFId andtenna for a 100 ns- window, b zoom of the begin- 
ning of the waveform, c zoom of the settling of the waveform 
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6 BMR Characterization of Antennas through the Combined TD/FD Approach 



Fig. 6.4 Antenna mea- 
surements performed in 
anechoic chamber [5] 




Table 6.1 rmse values, over the 8 10-960 MHz frequency range, between |Sn(/)| measured 
(for different acquisition windows) through the TD/FD combined approach and \S\\ re f(/)l 



7w 
(ns) 


rmse 


100 


0.040 


50 


0.037 


40 


0.019 


35 


0.019 


32 


0.010 


20 


0.021 


10 


0.068 



6.4.1.1 Lower Limit for the Acquisition Time Window 



As can be seen from Fig. |6. 5(a)] and Fig. |6. 5(b)} the S\\{f) curve corresponding to 
a 10 ns-long time window is far from representing the actual antenna behavior. This 
result can be easily explained considering that according to J6.ll ). the frequency res- 
olution corresponding to T„ = 10 ns is approximately 100 MHz, which is not even 
adequate to resolve the two closely-spaced resonant frequencies that characterize 
the AUT. In fact, looking at the specifications, it is known (and it is confirmed by 
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Fig. 6.5 Comparison between the S\i(f) measured through the TD/FD combined approach 
(for different time windows) and the Sn. re f(/). a magnitude, b phase 



the 5 , ii, re f(/)), that the AUT presents two resonant frequencies in its operating fre- 
quency range. The analysis of the FD-transformed data in the 8 10-960 MHz range, 
confirmed that for T w = 10 ns, T w = 20 ns and 7 W = 30 ns, the corresponding fre- 
quency resolution ensures the presence of one, three, and four spectral components, 
respectively. As a matter of fact, for T w > 20 ns the frequency resolution seems to be 
adequate for discriminating the two resonant peaks. Additionally, Fig. [63(b)] shows 
that at low times (< 20 ns), multiple reflections effects are still considerable and 
thus they cannot be windowed out, since the related spectral content is relevant. The 
frequency resolution requirements, along with the aforementioned considerations, 
suggest that the lowest limit for the windowing is T w low > 30 ns. 
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6.4.1.2 Upper Limit for the Acquisition Time Window and Interferences 
from Spurious Reflections 

As discussed in the previous section, a longer time window provides both a better 
frequency resolution and a more complete spectral representation of the response of 
the device under test. 

Nevertheless, considering Fig. 16.51 a longer time window does not necessarily 
guarantee better results, as can be seen, for example, observing the Su(f) curve 
corresponding to the lOOns-long time window. In fact, when time window is in- 
creased, noise contribution becomes more relevant; this is the intrinsic measurement 
limitation for the used instrument. 

To verify that the corruption of the measured waveform occurring at longer times 
(Fig. |6.3(c)[ ) is mostly due to noise contributions, the TDR response of the AUT 
was simulated through a full-wave simulation software, CST Microwave Studio. 
The simulation results are reported in Fig. l6.6l for the longer times, and they show a 
good agreement with the measurement data. As expected, the only difference is the 
noise absence, since the ideal condition is reproduced. 
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Fig. 6.6 Simulated TDR waveform of the RFId antenna 



As mentioned in the previous subsections, another critical issue related to the 
upper limit for the acquisition window is related to the possible presence of objects 
in the nearby of the AUT. 

To verify the influence of the spurious reflection, a second set of measurement 
was performed outdoor. Another set of S\ i (/) measurements was performed using a 
fixed time window and placing a reflecting object (a steel plate) at a varying distance 
from the AUT: this allowed making considerations on the interferences occurring, 
for example, when measurements are performed inside a laboratory. 
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After arbitrarily setting the time window duration at 27 ns, different waveforms 
were acquired placing a reflecting steel-made plate in front of the AUT, at a progres- 
sively increasing distance (from 1 m to 4 m, with step of 1 m). The steel plate was 
placed in the direction of maximum radiation of the AUT (i.e., in front of the AUT). 
Then the scattering parameter was evaluated: results are shown in Fig. 16.71 
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Fig. 6.7 Scattering parameter measurement with a steel plate placed at distance d in front of 
the antenna 



As expected, as the steel plate is taken far from the antenna, the S\ i (/) curves 
shift towards the curve of the S\ \ (f) when the steel plate is not present. In particular, 
the Sn (/) curve corresponding to a distance of 4 m from the antenna is practically 
superimposed to the Sn(/) of the antenna in free space: the steel plate seems to 
have no effect at all. Indeed, with reference to Fig. |6.2(b)| and according to ( I6.2I ). 
the distance d corresponding to a time window of 27 ns (which includes the 2.8 ns- 
portion that corresponds to the cable) is approximately 3.6 m. 

The interference of surrounding objects can actually be seen directly in the TD. 
In this regard, Fig. l6.8l shows the comparison between the TDR waveforms with the 
steel plate placed at approximately 1 m from the AUT and without. 

The two curves overlap up to approximately 10.4 ns (including 2.8 ns of the ca- 
ble): at this point some differences are clearly detectable. Indeed, according to ( 16.21 ). 
10 ns correspond to a distance of approximately 1 m. This validates the fact that in- 
terferences due to the presence of the reflecting objects in an indoor condition may 
actually be left out through an appropriate choice of the time window. 

Finally, additional measurements were performed comparing TDR waveforms 
acquired outdoor (making sure that no reflecting object was around) and indoor. As 
expected, both waveforms appear noisy at longer times, but the indoor-waveform 
appears corrupted by some additional spurious signals that are definitely attributable 
to the reflections from the environment: this effect starts to become relevant for a 
corresponding distance of 5 m, which is approximately the distances of the wall in 
the room. 
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Fig. 6.8 Comparison between TDR waveforms of the antenna, with and without the steel 
plate placed at a distance of 1 m from the AUT 



Therefore, referring to the TDR waveform reported in Fig. |6. 3(c)] it appears that 
the upper limit for the time windowing in order to ensure a SNR> 1 is approximately 
Tw.up <35ns. SNR measurements, directly available on the used instrument, con- 
firmed these theoretical considerations. 

Indeed, longer time windows can be used, provided that the influence of noise 
is reduced through an appropriate post-processing of the acquired TDR waveform. 
In particular, a time domain filter can be applied before the FD-transformation. The 
feasibility of this approach was confirmed by the results achieved applying a 5th 
order-moving average filter (available within the MATLAB filtering toolbox) to the 
TD data of a 100 ns-long time window. Fig. l6.9l shows the comparison between the 
>Sii,ref(/) curve and the S\\ (/) curves obtained with and without the application of 
the filter. 

As expected, the S\ \ (/) curve obtained from the filtered TDR waveforms pro- 
vides a more truthful picture of the AUT response. This is confirmed by the evalua- 
tion of the rmse between the magnitudes of Si i (/) (with and without the application 
of the filter) and of S\\ re f(/) : the rmse is 0.01 for the filtered data and 0.04 for the 
non-filtered data. 

These results suggest that the choice of a long time window, such as the consid- 
ered 100 ns-long time window, in conjunction with an appropriate time-domain fil- 
tering, can efficiently circumvent the noise limitation, and can simultaneously guar- 
antee a better frequency resolution. Additional experimental results (not reported 
here for the sake of brevity), confirmed that a long-time windowing, together with 
the aforementioned filtering procedure, is definitely the most appropriate choice 
when the frequency response of the device must be retrieved also at the lowest 
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Fig. 6.9 Comparison between the |5ii(/)| measured through the TD/FD combined approach 
(with and without the application of the digital filter) and the iSn re f(/)l 

frequencies. On the other hand, when exploring a limited frequency bandwidth, as 
in the considered test-case, the optimal time windowing (for the considered AUT) 
is 30ns< r w <35ns. 



6.5 Biconical Antenna Results 

As a further experimental validation of the method, the second considered antenna 
was a biconical antenna (Fig. |6.1(b)| i, generally used for Electromagnetic Compati- 
bility (EMC) measurements: the wideband characteristics of this antenna anticipate 
different criteria for the choice of the optimal time window d 1 511 . 

Similarly to the previous case, the TDR measurements on the antenna were per- 
formed by choosing several time windows, ranging from 15 ns to 150 ns; once again, 
the corresponding Su(f) was extrapolated and compared to the 5n^ e f(/) m terms 
of rinse (Table [672b . 



Table 6.2 rmse values between |Sn(/)| measured (for different acquisition windows) 
through the TD/FD combined approach and \S\\ re f(/)|> for the biconical antenna 



7w 
(ns) 


rmse 


150 


0.010 


100 


0.006 


80 


0.007 


40 


0.010 


35 


0.013 


15 


0.019 
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Fig. 6.10 TDR waveform of the biconical antenna for a 100 ns-long time window (left), zoom 
of the initial portion of the waveform (right) [3] 




Fig. 6.11 Comparison between |5 , ii re f(/)| and the |Sn(/)| curves obtained from FD- 
transformed data for different time windows (r w ), for the biconical antenna yO 



The analysis of the TDR waveform of the biconical antenna (reported in 
Fig. 16.10b clearly shows that, differently from the previous case, the multiple re- 
flections quickly die out and the signal reaches a steady-state condition in about 
20 ns. Nevertheless, such a short T w would not accurately represent the antenna re- 
sponse at lower frequencies. In fact, Fig. 16. Ill shows that even a 30 ns-long time 
window, although providing overall accurate results, fails in accurately representing 
the antenna response at low frequencies. On the other hand, if the time window is 
too long, this will translate in inaccurate results at high frequencies. In fact, since the 
maximum number of acquisition points is fixed, then the sampling period (T c ) may 
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become too high, and hence the sampling frequency too low: this represents the up- 
per limit for the maximum frequency /m- In order to retrieve the frequency response 
up to /m = 3 GHz, then f$ should be at least 6 GHz, so as to satisfy the sampling 
theorem. As a result, the corresponding T c should be lower than 0.166 ns. Consider- 
ing that, for a fixed time window, the used instrument provides 4000 measurement 
points in TD, the maximum time window must be shorter than 666 ns. 

Indeed, the maximum time window should be definitely shorter. In fact, the used 
TDR instrument performs a real-time sampling only for time windows shorter than 
120 ns: beyond this limit, the instrument samples in equivalent- time. As a matter of 
fact, equivalent-time sampling is inappropriate when such one-shot signals are con- 
sidered. Therefore, it is advisable to use windows shorter than 100 ns. Additionally, 
it is worth pointing out that, in this case, the intrinsic noise limitation is not critical, 
since the steady-state condition is not reached through an oscillating and attenuating 
transient. 

The aforementioned considerations may be summarized as follows: if T w is too 
short, then results will be less accurate at low frequencies, conversely, if T w is too 
long, this will result in inaccurate results at high frequency. 

On such basis, Fig. 16. Ill shows that a lOOns-long time window provides over- 
all accurate results over the entire considered frequency range. The evaluation of 
the rmse values (reported in Table [672b confirmed that a lOOns-long time window 
(for which real-time sampling is assured), provides the best results in terms of ac- 
curacy. The rmse value corresponding to a 150ns-long time window confirms the 
performance degradation for higher time windows. 
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